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HIGHLIGHT  SUMMARY 

Information  about  the  soil  permeability  is  vital  to  the  sucess 
of  many  geotechnical  works.   Various  predictive  equations  based  on 
volumetric  and  grain  size  parameters  have  been  proposed  to  estimate 
the  permeability.   These  predictive  equations  have  had  limited  success 
for  coarser  soils,  but  have  been  completely  unsatisfactory  for  fine- 
grained soils.  Knowledge  of  the  soil  structure  is  essential  to 
understanding  and  predicting  the  engineering  behavior  of  fine-grained 
soils.   The  volumetric  and  grain  size  parameters  are  poor  indicators 
of  soil  structure  and  cannot  be  expected  to  accurately  reflect  changes 
in  permeability.  This  investigation  utilizes  pore  size  distribution 
measurements  to  examine  the  relationship  between  permeability  and  soil 
fabric  of  laboratory  compacted  clayey  silts. 

A  method  of  approximating  the  differential  pore  size  distribution 
is  described  which  simplifies  the  interpretation  of  pore  size 
distributions.   The  most  frequent  or  modal  pore  diameters  served  as 
useful  comparative  parameters  when  contrasting  pore  size  distribution 
curves  for  various  soils. 

Three  blends  of  silt  and  kaolin  were  compacted  by  kneading 
compaction  at  two  or  three  levels  of  effort.   For  each  soil  type  and 
compactive  effort,  samples  were  prepared  at  three  water  contents: 


optimum,  dry  of  optimum,  and  wet  of  optimum.   Closed  system  falling 
head  permeability  tests  under  back  pressure  were  performed  on  each  of 
the  compacted  samples.   Freeze  drying  was  successfully  used  to 
dehydrate  specimens  prior  to  pore  size  measurements.   The  mercury 
intrusion  technique  was  used  to  determine  the  pore  size  distribution. 

The  pore  size  distributions  of  the  soils  tested  were  bimodal  with 
a  large  pore  mode  occurring  between  10  and  1  pm  and  a  small  pore  mode 
occurring  at  0.1  ym.   Varying  the  compaction  variables  produced  changes 
in  the  position  and  frequency  of  the  large  pore  mode,  but  caused  no 
change  in  the  pore  size  distribution  about  the  small  pore  mode  for  a 
given  soil  type.   Permeability  generally  increased  with  increases  in 
the  frequency  and  position  of  the  large  pore  mode. 

Three  theoretical  permeability  models  which  relate  pore  size 
distribution  parameters  to  permeability  were  examined.   Although  the 
models  were  not  completely  successful  in  predicting  permeability,  the 
pore  size  parameters  from  two  of  the  models  did  reflect,  with  several 
modifications,  the  influence  of  pore  size  distribution  on  permeability. 
These  parameters  were  successfully  used  to  determine  empirical 
permeability  prediction  equations  for  the  soils  tested. 

This  study  demonstrates  the  usefulness  of  pore  size  distribution 
measurements  to  interpret  changes  in  soil  fabric  and  to  predict  how 
these  changes  affect  an  engineering  property. 


INTRODUCTION 

Permeability  is  one  of  the  most  important  yet  variable  properties 
in  Geotechnical  Engineering.   Conventionally,  permeability  is  evaluated 
either  by  in  situ  or  laboratory  tests.   In  situ  tests  are  costly  and 
difficult  to  interpret  because  of  the  complex  boundary  conditions 
involved.   Laboratory  tests  also  present  formidable  problems,  the  most 
troublesome  being  obtaining  representative  samples.   To  circumvent 
the  testing  problems,  methods  of  estimating  permeability  from  grain 
size  and  volumetric  characteristics  have  previously  been  attempted. 
These  predictive  equations  have  had  limited  success  for  coarser  soils, 
and  have  proved  to  be  completely  unsatisfactory  for  fine-grained  soils. 
The  fundamental  problem  with  these  equations  is  that  permeability  is 
extremely  sensitive  to  subtle  changes  in  soil  structure,  a  property 
which  bulk  volumetric  parameters  do  not  accurately  reflect. 

This  investigation  takes  a  more  fundamental  approach  by  examining 
the  relationship  between  permeability  and  pore  size  distribution  of 
compacted  fine-grained  soils.   Three  models  of  porous  media,  which 
utilize  pore  size  distribution  parameters  to  calculate  permeability, 
are  examined.   Simple  linear  regression  is  used  to  determine  permea- 
bility predicting  equations  from  these  pore  size  parameters  and 
permeability  measurements. 


The  fabric  of  compacted  soils  is  also  examined.   The  influence  of 
varying  water  content  and  compaction  effort  on  the  pore  size  distri- 
bution of  three  soils  is  considered  in  light  of  existing  models  for  the 
structure  of  compacted  clay. 

Pore  size  distributions  were  determined  by  the  method  of  mercury 
intrusion.  Dehydration  of  the  soils  prior  to  pore  size  measurements 
was  accomplished  by  freeze  drying.   Both  differential  and  cumulative 
pore  size  distributions  are  presented. 

A  closed  system  falling  head  test  was  used  to  measure  the 
permeability  of  the  compacted  samples.   Back  pressure  was  applied  to 
the  samples  during  permeation  to  insure  saturation,  and  the  permeability 
tests  were  conducted  at  low  hydraulic  gradients. 


1  -  LITERATURE  REVIEW  AND  GENERAL  DISCUSSION 

1-1  Permeability  and  Fabric  of  Compacted  Fine-Grained  Soils 

Permeability  is  one  of  the  most  fundamental  yet  variable  properties 
of  soils  in  Geotechnical  Engineering.   Permeability  must  be  considered 
in  problems  of  seepage,  drainage,  dewatering,  rate  of  settlement  and 
rate  of  strength  increase  with  pore  pressure  dissipation.  Values  of 
soil  permeability  may  range  across  twelve  orders  of  magnitude  and  are 
extremely  sensitive  to  changes  in  soil  composition  and  structure. 

The  measurement  of  soil  permeability  in  the  laboratory  presents  a 
formidable  problem,  but  when  these  laboratory  results  are  extrapolated 
to  field  permeability  values  for  analysis  purposes,  enormous  if  not 
intolerable  variability  may  result.   The  primary  reason  for  this  is 
that  there  exists  no  reliable  way  of  determining  the  field  soil 
structure  and  thus,  no  way  to  replicate  it  in  the  laboratory.   The 
macrofabric  of  the  in  situ  soil  must  also  be  considered  (Rowe,  1972), 
but  for  compacted  soils,  these  features,  such  as  fissures,  joints,  and 
stratification,  are  not  as  prevalent. 

The  objectives  of  this  study  are  to  use  pore  size  distribution 
measurements  as  an  indicator  of  soil  fabric: 

1)  to  improve  the  understanding  of  the  fabric  of  compacted 
fine-grained  soils,  and 

2)  to  determine  a  relation  between  pore  size  distribution  and 
permeability  measurements  for  these  same  soils. 


This  section  reviews  the  results  of  previous  studies  concerned  with 
the  permeability  of  compacted  clays,  discusses  the  importance  of  soil 
fabric  of  compacted  clays,  and  examines  several  permeability  prediction 
models  which  incorporate  pore  size  distribution  parameters. 

1-1.1  Permeability 
Larabe  (1954)  was  one  of  the  first  investigators  to  examine  the 
permeability  of  compacted  fine-grained  soils.   He  noted  the  pronounced 
affect  that  molding  water  content  had  on  the  permeability  for  a  given 
soil  compacted  at  a  given  effort.   Lambe  concluded  that  since  these 
changes  in  permeability  could  not  be  accounted  for  by  changes  in  dry 
density  or  void  ratio,  they  must  be  caused  by  changes  in  the  "arrange- 
ment of  soil  particles  or  'structure'."  Subsequent  studies  of  the 
permeability  of  compacted  clays  and  silts  to  water  by  Bjerrum  and  Huder 
(1957),  and  Mitchell,  Hooper  and  Campanella  (1965)  found  results 
similar  to  those  obtained  by  Lambe.   Figure  1,  taken  from  Mitchell  et  al. 
(1965),  shows  the  relationship  between  compaction  water  content  and 
permeability  of  a  silty  clay  prepared  at  three  compaction  efforts.  The 
relation  shown  by  this  Figure  may  be  considered  typical  for  compacted 
fine-grained  soils.   Summarizing  the  results  of  the  three  investigations 
listed  above  (Lambe,  L;  Bjerrum  et  al. ,  B;  and  Mitchell  et  al.  ,  M)  : 

1.   The  permeability  gradually  decreases  with  increasing  water 
content  on  the  dry  side  of  optimum  for  a  given  compaction 
effort  (L,  B,  M) ;  though,  Mitchell  et  aL  reported  one  soil 
for  which  this  was  not  the  case. 
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FIGURE    I       PERMEABILITY    OF  A  SILTY  CLAY  AT  THREE  COMPACTIVE 
EFFORTS    (FROM   MITCHELL  etal  1965) 


2.  Near  optimum,  the  permeability  undergoes  a  significant 
decrease  of  two  to  three  orders  of  magnitude  less  than 
the  dry  side  permeabilities  with  increasing  water  content 
for  a  given  compaction  effort  (L,  B,  M) . 

3.  Wet  of  optimum,  the  permeability  may  increase  or  decrease 
slightly  with  increasing  water  content  for  a  given 
compaction  effort,  but  generally  it  remains  within  the 
same  order  of  magnitude  as  the  permeability  at  optimum 
(L,  B,  M). 

4.  Increasing  the  degree  of  saturation  by  back-pressuring 
during  permeation  slightly  increases  the  permeability, 
all  other  factors  being  constant  (B,  M) . -  This  effect 
is  more  pronounced  on  dry  side  samples  (B) .   The 
permeability  increases  approximately  with  the  cube  of 
the  degree  of  saturation,  all  other  factors  being 
constant  (M) . 

5.  Varying  the  compaction  effort  shifts  the  position  of  the 
permeability-molding  water  content  curve,  but  the  general 
trends  (1  to  3  above)  remain  the  same,  as  shown  in  Figure 
1  (M). 

6.  For  a  given  water  content  wet  of  optimum  and  a  given  dry 
density  and  soil  type,  static  compaction  produces  a 
sample  with  a  higher  permeability  than  one  compacted  by 
kneading  compaction  (M) . 

Langfelder,  Chen  and  Justice  (1968)  measured  air  permeabilities  of 
compacted  fine-grained  soils  and  concurred  with  results  1  and  2  listed 


above.   They  also  determined  that  for  water  contents  at  or  above 
optimum,  the  air  permeability  was  essentially  zero.   For  most  of  the 
soils  which  were  tested  at  varying  gradients,  an  "apparent  threshold 
gradient"  was  obtained,  whose  magnitude  increased  as  the  air 
permeability  decreased. 

The  Kozeny-Carman  equation  (Carman,  1956)  has  frequently  been  used 
to  predict  the  permeability  of  other  types  of  porous  media  and  is  given 
as: 
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where  k  is  the  permeability  (L/t) 

Y  is  the  unit  weight  of  the  permeant 

U  is  the  absolute  viscosity  of  the  permeant 

e  is  the  void  ratio 

S  is  the  internal  surface  area  per  volume  of  solids 

k  is'  the  pore  shape  factor 

k°  is  the  tortuosity 

This  expression  has  not  been  found  to  be  suitable  for  predicting  the 
permeability  of  fine-grained  soils  as  shown  by  Lambe  (1954),  Michaels 
and  Lin  (1954),  and  Olsen  (1962). 

Michaels  and  Lin  (1954)  performed  permeability  tests  on  kaolin  beds 
prepared  and  permeated  with  fluids  of  varying  polarity  and  compressed 
to  decreasing  void  ratios.   They  found  that  the  rate  of  change  of  per- 
meability was  greater  than  that  predicted  by  the  Kozeny-Carman  equation. 
Because  these  deviations  were  of  the  same  order  for  the  various  fluids 
tested,  they  concluded  that  the  effect  was  caused  "...mainly  by  changes 
in  aggregate  size  and  particle  orientation  during  mechanical  compaction", 
rather  than  by  interactions  between  the  fluids  and  the  particle  • 
surface.   The  permeability  of  kaolin  decreased  markedly  with  increasing 
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polarity  of  the  fluid.   They  attributed  this  to  the  existence  of  a 
more  dispersed  structure  resulting  from  an  increase  in  the  polarity  of 
the  fluid  in  which  the  clay  beds  were  prepared.   Finally,  they 
concluded  that  interfacial  effects  such  as  adsorption  and  counter- 
electroosmosis  are  responsible  for  only  minor  changes  in  the 
permeability  of  kaolin. 

Olsen  (1962)  conducted  one  dimensional  consolidation-permeability 
tests  on  kaolinite,  illite  and  Boston  blue  clay  and  also  found 
discrepancies  between  the  permeabilities  calculated  from  the  Kozeny- 
Carman  equation  and  measured  permeabilities  as  follows: 

1)  For  void  ratios  greater  than  about  0.4  to  0.5  the 
permeability  decreased  more  than  predicted  with 
decreasing  void  ratio. 

2)  For  void  ratios  less  than  0.4  the  permeability 
decreased  less  than  predicted  with  decreasing 
void  ratio. 

3)  For  rebounded  samples  the  permeability  increased 
at  a  rate  less  than  predicted  with  increasing 
void  ratio. 

Olsen,  after  considering  if  the  above  discrepancies  could  be 
explained  by  non-Darcy  flow,  electrokinetic  coupling,  high  viscosity, 
or  tortuous  flow  path,  concluded  that  only  unequal  pore  sizes  could 
account  for  the  discrepancies.   Olsen  hypothesized  a  cluster  or  packet 
and  domain  model  for  the  soil  structure  which  would  explain  the  unequal 
pore  sizes.   Recent  investigations  of  fine  grained  soil  fabric  have 
confirmed  Olsen' s  cluster  hypothesis,  as  will  be  discussed  below. 


Several  investigators  including  Hansbo  (1960),  Swartzendruber 
(1962),  Miller  and  Low  (1963),  Olsen  (1965)  and  (1966),  Mitchell  and 
Younger  (1967),  and  Russell  and  Swartzendruber  (1971)  have  investigated 
the  validity  of  Darcy's  law  for  fine-grained  soils  at  low  hydraulic 
gradients.   Barring  experimental  errors,  two  explanations  have  been 
proposed  to  account  for  the  non-linear  relation  between  flow  velocity 
and  hydraulic  gradient: 

1)  particle  migration  causing  clogging  and  unclogging  of  the 
soil  pores  (Hansbo,  1960); 

2)  quasi-crystalline  behavior  of  water  near  the  clay  surface 
(Miller  and  Low,  1963). 

Some  insight  can  be  gained  by  examining  hypothetical  soil  structures 
which  might  explain  each  of  the  above  causes  for  non-Darcy  flow. 

For  particle  migration  to  be  a  factor,  a  small  fraction  of  the 
smaller  size  soil  particles  (presumably  non-clay)  must  exist  in  a 
loose  state  within  the  soil  mass.   These  loosely  held  particles  should 
be  roughly  the  same  size  or  slightly  smaller  than  the  larger  pores 
which  control  the  permeability  of  the  porous  medium.   Mitchell  and 
Younger  (1967)  found  that  for  a  compacted  silty  clay,  particle  migration 
caused  changes  in  soil  fabric  during  flow  which  led  to  non-Darcy 
effects;  these  effects  were  most  significant  for  samples  compacted  at 
low  water  contents  (significantly  lower  than  optimum)  and  low  densities. 
Such  soils  probably  have  a  relatively  high  permeability,  as  discussed 
previously,  and  an  unstable  structure  (susceptible  to  collapse  or 
piping). 

For  quasi-crystalline  water  to  be  responsible  for  non-Darcy  flow 
at  low  hydraulic  gradients,  the  pores  controlling  the  permeability 
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must  be  of  size  comparable  to  the  thickness  of  the  adsorbed  water 
layer.  This  would  mean  a  tightly  packed  soil  fabric  (with  small  pore 
size)  containing  a  clay  of  high  activity  (Miller  and  Low,  1963).   Such 
a  soil  would  be  expected  to  have  a  low  permeability. 

The  factors  affecting  the  permeability  of  soils  and  other  porous 
media  have  been  discussed  by  a  number  of  authors,  including  Taylor 
(1948),  Lambe  (1954),  Scheidegger  (1957)  and  Leonards  (1962).   In 
general  terms,  three  factors  control  the  permeability  of  a  soil: 

1)  the  geometry  of  the  porous  network  (i.e.  the  fabric); 

2)  the  properties  of  the  permeating  fluid  (viz.,  viscosity  and 
density) ; 

3)  the  surface  interaction  between  the  permeating  fluid  and 
the  porous  media. 

From  the  various  investigations  just  described,  it  is  clear  that 
the  geometry  of  the  porous  network  is  by  far  the  most  important  factor 
controlling  the  permeability  of  most  fine-grained  soils.   The  second 
factor  above  is  of  little  consequence  in  geotechnical  engineering  since 
water  is  nearly  always  the  pore  fluid.   The  influence  of  the  third 
factor  on  permeability  should  certainly  be  significant  for  soils  of 
high  plasticity  (Schmid,  1957)  or  dispersive  clays;  however,  this  topic 
was  not  addressed  in  this  investigation.   Blends  of  silt  and  kaolin  of 
low  plasticity  were  used  in  an  effort  to  minimize  the  influence  of 
interactions  between  fluid  and  soil  surface  on  the  permeability. 

Because  of  the  importance  of  soil  fabric  on  permeability,  a  brief 
review  of  some  previous  investigations  dealing  with  the  structure  of 
compacted  clays  and  silts  is  appropriate. 
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1-1.2  Fahric 

Prior  to  commencing  this  discussion  of  soil  fabric  it  is  necessary 

to  define  the  terms  "fabric"  and  "structure"  as  used  in  this  study. 

According  to  Yong  and  Sheeran  (1973):  , 

"The  structure  of  a  soil  has  been  defined  as  that 
property  of  soil  which  provides  its  integrity.   An 
important  component  of  structure  is  its  fabric,  i.e.  the 
physical  arrangement  of  soil  particles  including  particle 
spacing  and  pore  size  distribution.   When  soil  fabric  is 
considered  in  conjunction  with  bonding  forces  and 
particle  interaction  mechanisms  developed,  the  structure 
of  soil  is  obtained  ..." 

The  term  "macrofabric",  according  to  Mitchell  (1976),  will  refer  to  "... 

those  features  that  can  be  seen  with  the  unaided  eyes,  or  a  hand 

lens  ..."  including  "stratification,  fissuring,  voids,  and  large 

inhomogeneities  ..."  The  importance  of  "macrofabric"  in  geotechnical 

engineering  is  discussed  by  Rowe  (1972),  and  will  not  be  considered  in 

this  study. 

Lambe  (1958)  was  one  of  the  first  investigators  to  hypothesize  a 

model  for  the  structure  of  compacted  clay.   He  based  his  model  on  the 

principles  of  colloidal  and  crystal  chemistry  and  the  behavior  of 

compacted  clays.   Lambe  postulated  that  individual  clay  particles  are 

the  predominant  units  which  influence  the  compaction  characteristics 

and  behavior  of  a  soil  mass.   At  low  water  contents  a  compacted  soil  is 

in  an  open  and  floccuated  state  because  the  suppressed  double  layer  has 

reduced  particle  repulsion.   As  the  molding  water  is  increased  toward 

optimum,  the  double  layer  expands  which  increases  particle  repulsion, 

but  a  more  orderly  particle  arrangement  is  formed  due  to  increased 

"lubrication".   At  water  contents  greater  than  optimum,  the  diffuse 

double  layer  continues  to  expand,  increasing  particle  repulsion, 
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increasing  the  distance  between  particles,  and  resulting  in  a  nearly 

parallel  particle  arrangement. 

More  recent  investigations  of  soil  structure  employing  scanning 

electron  microscopy,  X-ray  diffraction,  and  pore  size  distribution 

measurements  have  resulted  in  modifications  of  Lambe's  compacted  soil 

model.   As  Yong  and  Sheeran  (1973)  stated: 

"Examination  of  electron  micrographs  and  soil  per- 
formance shows  that  individual  particles  rarely  act  as 
single  particle  units,  except  for  certain  types  of  clays 
and  under  certain  conditions,...   The  different  sizes 
and  arrangements  of  particle  groups  observed  in  fabric 
viewing  suggest  that  response  behavior  might  be 
controlled  by  the  kinds  and  arrangements  of  these 
particle  groups." 

Independent  studies  by  Barden  and  Sides  (1970)  and  Hodek  (1972) 
have  proposed  a  revision  of  Lambe's  compacted  clay  model. 

Barden  and  Sides  (1970)  combined  microscopic  fabric  analysis  with 
laboratory  measurements  to  study:  1)  the  development  and  drainage  of 
pore  pressures  and  2)  the  compression  characteristics  of  compacted 
clays. 

Hodek  (1972)  examined  the  strength  and  particle  orientation 
characteristics  of  aggregations  of  kaolin  formed  by  conventional  mixing 
techniques  prior  to  compaction.   He  then  analyzed  the  influence  of 
these  kaolin  aggregates  at  various  water  contents  and  of  various  sizes 
on  the  compaction  characteristics  and  swelling  pressure  of  compacted 
specimens. 

Both  investigations  reached  the  conclusion  that  the  behavior  and 
characteristics  of  a  compacted  clay  can  be  explained  in  terms  of  a 
deformable  aggregate  soil  model.   Prior  to  compaction,  the  soil 
particles  are  grouped  in  agglomerations  or  "peds"  whose  size  and 
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strength  characteristics  are  influenced  by  the  molding  water  content. 
During  compaction  at  low  water  contents  (below  optimum)  these  peds  or 
aggregates  have  a  high  strength  and  are  better  able  to  resist  the 
compaction  pressures  without  much  distortion.   Thus,  there  exist  two 
networks  of  pore  space  in  the  clay  mass:  a  network  of  large  inter- 
aggregate  pores  and  a  network  of  small  intra-aggregate  pores.  As  the 
molding  water  content  increases  (dry  of  optimum),  the  aggregates 
decrease  in  strength  and  suffer  greater  deformation  during  compaction. 
This  results  in  a  decrease  in  inter-aggregate  pore  space  and  an  increase 
in  dry  density.   Olson  (1963),  Langf elder  et  aL  (1968),  and  Barden  and 
Sides  (1970)  have  all  found  that  clays  compacted  near  their  optimum 
water  contents  have  zero  air  permeability,  which  would  indicate  that 
the  inter-aggregate  pore  network  is  no  longer  continuous.  With 
increasing  water  content  above  optimum,  the  aggregates  are  easily 
distorted  and  fuse  together  making  them  indistinguishable.   Individual 
particle  reorientation  and  dispersion  may  occur  at  this  stage. 

Pore  size  distribution  measurements  of  compacted  clays  by  Diamond 
(1970)  and  (1971),  Sridharan  et  al.  (1971),  Ahmed  et  aL  (1974)  and 
Bhasin  (1975)  have  also  provided  strong  evidence  for  a  deformable 
aggregate  soil  model. 

Diamond  (1971)  employed  pore  size  distribution  methods,  scanning 
electron  microscopy,  and  X-ray  diffraction  to  study  the  fabric  of  impact 
compacted  kaolinite  and  illite.   Pore  size  distribution  measurements 
indicated  that  clays  compacted  dry  of  optimum  have  a  significant  volume 
of  pores  between  1  and  10  Um  in  size  while  samples  compacted  at  or  above 
the  optimum  water  content  contained  very  little  pore  space  greater  than 
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0..1  ym  in  size.   This  conclusion  was  also  reached  by  Ahmed  et  al.  (1974) 
and  Bhasin  (1975)  for  other  types  of  clayey  soils  and  using  improved 
soil  drying  techniques  for  pore  size  measurements.   Scanning  electron 
micrographs,  presented  by  Diamond  (1971),  revealed  individual 
aggregations  of  clay  particles  several  microns  in  size  and  corresponding 
voids  of  the  same  size  for  samples  compacted  dry  of  optimum.   Samples 
compacted  at  or  above  the  optimum  water  content  did  not  show  these 
individual  aggregations  nor  the  large  pore  space.   Using  X-ray 
diffration  measurements  of  oriented  compacted  clay  samples,  Diamond 
found  no  evidence  of  increased  particle  orientation  with  increasing 
molding  water  content,  as  hypothesized  by  Lambe  (1958). 

Sridharan  et  al.  (1971)  performed  pore  size  distribution  measurements 
of  statically  compacted  kaolin,  illite,  and  Boston  blue  clay  and 
artifically  sedimented  kaolin.   Static  compaction  of  a  given  soil  at  a 
given  water  content  to  decreasing  void  ratio  values  was  found  to  be  at 
the  expense  of  the  larger  pores  with  little  or  no  influence  on  pores 
about  0.1  urn  in  size.   They  reported  that  each  type  of  clay  had  its  own 
characteristically  different  pore  size  distribution  when  compacted  at 
similar  water  contents.   There  appeared  to  be  no  obvious  correlation 
between  the  pore  size  and  grain  size  distributions.   The  shape  of  the 
pore  size  distribution  curve  was  also  found  to  be  sensitive  to  the 
method  of  soil  preparation.   Artificially  sedimented  kaolinite  had  a 
pore  size  curve  distinctly  different  from  that  for  statically  compacted 
kaolinite. 

Ahmed  et  al.  (1974)  found  no  significant  difference  in  the  shapes 
of  pore  size  distribution  curves  between  clays  prepared  by  impact  and 
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by  kneading  compaction  to  a  common  water  content  and  dry  density. 

Bhasin  (1975)  found  that  increasing  the  compactive  effort  on  the 

dry  side  of  optimum  for  a  given  soil  decreases  the  total  porosity  and 

diminishes  the  fraction  of  large  pores.   Increasing  the  compaction 

effort  on  the  wet  side  of  optimum  had  little  effect  on  either  the  total 

porosity  or  the  distribution  of  pore  sizes.   He  also  concluded  that: 

"The  gross  differences  in  pore  size  distributions 
occurring  at  the  same  percentage  compactions,  for 
different  soils  and  compaction  efforts,  emphasize  the 
lack  of  control  over  the  compacted  product  exercised  by 
most  end  result  compaction  specifications." 

To  date,  no  pore  size  distribution  measurements  have  been  reported 
for  field  compacted  fine-grained  soils.   Certainly,  the  sensitivity  of 
soil  fabric  to  slight  changes  in  water  content,  compaction  energy  and 
method  of  compaction  would  warrant  such  an  investigation. 

1-1.3  Permeability  Models 
The  aim  of  this  section  is  to  examine  several  of  the  relatively 
simple  and  conventional  models  used  to  describe  flow  through  porous 
media  in  order  to  gain  some  insight  into  the  relation  between 
permeability  and  pore  size  distribution.   The  final  objective  will  be 
to  arrive  at  several  pore  size  parameters,  which  can  be  calculated  from 
a  pore  size  distribution,  and  which  may  correlate  with  permeability 
measurements.   Scheidegger  (1974)  and  Bear  (1972)  both  discuss  the 
complex  nature  of  a  porous  media  and  the  limitations  of  conventional 
models  for  predicting  permeability.   They  conclude  that  statistical 
modeling  would  be  a  more  logical  approach  to  predict  permeability  and 
present  such  models;  however,  owing  to  its  complexity,  statistical 
modeling  of  porous  media  was  considered  beyond  the  scope  of  this  study. 
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Two  classical  models  have  been  employed  to  describe  the  flow 
through  porous  media,  viz.,  the  equivalent  capillary  model  and  the 
equivalent  hydraulic  radius  model.   Both  can  be  derived  from  the 
Hagen-Poiseuille  equation  as  demonstrated  by  Leonards  (1962). 

The  Hagen-Poiseuille  equation  for  laminar  flow  through  a  cylindrical 
capillary  is  given  as: 

71  Y  S  „4 

q  =  ITT-  R  (i-2) 

where  q  is  the  volume  of  flow  per  unit  time 
R  is  the  radius  of  the  capillary 
S  is  the  hydraulic  gradient 
Y  is  the  fluid  unit  weight 
y  is  the  absolute  viscosity  of  the  fluid. 

Purcell  (1949)  employed  a  modified  capillary  model  to  calculate  the 

permeability  of  sandstones  from  pore  size  distribution  measurements. 

Purcell  modeled  the  porous  media  as  a  set  of  parallel  cylindrical 

capillaries  of  varying  pore  diameter  extending  from  one  end  of  the 

porous  medium  to  the  other.   Each  capillary  is  assumed  to  have  a 

constant  radius  throughout  its  length.   The  frequency  of  occurrence  of 

each  capillary  size  is  determined  from  the  pore  size  distribution  curve. 

Schiedegger  (pp.  129-130,  1974)  discusses  a  simplified  method  of 

calculating  the  permeability  from  this  model  and  also  discusses  some  of 

its  limitations.   Beginning  with  the  Hagen-Poiseuille  equation,  the  flow 

rate,  q^  through  a  capillary  of  radius  r.  is  given  as: 


8u 
Equation  (1-3)  can  be  rewritten  as: 


7T  y  S  4 

-8u—  ri  (1-3) 
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Y  S  2    2  ,,  ,n 

q.  =  -V"  r.  tt  r.  (1-4) 

For  more  than  one  capillary  of  the  radius  r . ,  let  a(r. )  represent  the 
number  of  pores  of  size  r..   Thus  the  flow  rate  through  all  capillaries 
of  r .  becomes : 


«i  =  8*  ri  a(ri>  «  4  (1"5> 

2 
The  term  "ct(r.)  tt  r."  is  the  "area-frequency"  of  the  capillary  radius 

r..   It  is  evident  that  for  any  cross-section  taken  perpendicular  to 
the  direction  of  flow,  the  area-frequency  of  r.  will  remain  the  same. 
Thus,  multiplying  the  area-frequency  by  a  unit  length  gives  the 
volumetric  frequency  f(r. )  of  r.  (Harr,  1962  p.  A),  which  is  the 
quantity  measured  by  pore  size  distribution  methods. 

Considering  a  unit  volume  of  pore  space  in  the  porous  medium  with 
a  distribution  of  parallel  capillaries  traversing  it,  the  volumetric- 
frequency  distribution  of  the  capillaries  can  be  given  as: 

I  f(r.)  •  1  (1-6) 

i 

Combining  equations  (1-5)  and  (1-6),  the  total  flow  rate  through  all  the 
capillaries  through  a  unit  volume  of  pore  space  is: 


*  -  ^  ? f  (r±>  A  (1-7) 


Dividing  equation  (1-7)  by  a  unit  area  of  pore  space  gives  the  average 
seepage  velocity  through  the  pores  v  : 
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-  -  *&  ?  f  (ri>  ri  Cl~8) 


a    8y  ± 


The  discharge  velocity  v  through  a  unit  section  of  porous  media  is: 


v  •>   v  n  (1-9) 

s 


v  =  ^|n  I   ff(r±)  v[  (1-10) 


where  n  is  the  porosity. 

Applying  Darcy's  Law,  "v  =  k  "  S",  to  equation  (1-10)  gives: 


1       , 
k=  gjn  I   f(r±)  x\  (1-lD 


where  k  is  the  permeability. 

1       2 
The  value  "£  f(r. )  r. "  is  defined  as  the  second  moment  about  the  origin 

1  2 

of  the  pore  size  distribution  and  is  usually  represented  as  "E(r  )". 

If  the  pore  diameter  is  considered  rather  than  the  pore  radius,  equation 

(1-11)  becomes: 


k  =  ifc  E(d2)  n  C1"12) 


Replacing  the  permeability  k  by  the  "specific  or  physical  permeability, 
K"  and  using  k  =  K  y/v»  equation  (1-12)  becomes: 

K  =  C  E(d2)  n  Cl-13) 

s 


where  C  is  a  shape  factor  which  is  equal  to  1/32  for  cylindrical 
pores. 
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There  are  several  shortcomings  of  this  equation,  as  pointed  out  by 
Scheidegger;  (1974): 

1)  All  pores  are  assumed  to  go  in  one  direction. 

2)  Tortuosity  is  not  considered. 

3)  The  shape  factor  is  considered  constant  for  all  pore  sizes. 

2 

4)  The  value  of  E(d  )  is  extremely  sensitive  to  the  frequency 

of  the  larger  pores. 

5)  The  model  is  an  oversimplification  of  the  porous  media. 

Considering  equation  (1-13)  as  an  empirical  equation  and  solving 

for  C  from  actual  pore  size  and  permeability  measurements,  Purcell 
s 

2 
found  that  changes  in  "E(d  )  n"  (or  its  equivalent)  accurately  reflected 

changes  in  permeability.  However,  the  value  of  C  (or  its  equivalent) 
varied  for  the  different  materials  tested,  and  for  a  given  material  at 
varying  porosities,  C  increased  with  increasing  permeability. 

Childs  and  Collis-George  (1950)  and  Marshall  (1958)  employed 
another  form  of  the  capillary  model  using  the  pore  size  distribution  to 
calculate  permeability.   They  considered  two  unit  cross-sections  of  a 
porous  medium  placed  together  such  that  the  pores  from  one  surface  are 
randomly  connected  to  those  of  the  other  surface.   Each  of  the  cross- 
sections  contain  cylindrical  pores  of  varying  radius.   The  pore  size 
distributions  for  each  of  the  two  cross-sections  are  assumed  to  be 
identical  and  represented  as: 

I   f(r.)  =  n  (1-14) 

i 

where  f(r.)  is  again  the  volumetric  frequency  of  occurrence  of  pore  r. 
and  n  is  the  porosity.  Note  that  a  unit  volume  of  porous  medium  and 
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not  a  unit  volume  of  pore  space  is  being  considered;  therefore,  the  sum 
of  the  volumetric  frequency  is  equal  to  the  porosity. 

The  probability  that  a  pore  r.  from  the  i   cross-section  is 
joined  with  a  pore  r.  from  the  j   cross-section,  according  to 
probability  theory,  can  be  represented  as: 


r±  -*■  r  =  f(r±)  f(rj)  (1-15) 


or  the  product  of  the  two  probabilities  as  each  is  independent  of  the 
other.  The  flow  rate  through  the  capillary  connecting  r.  and  r.  from 
equation  (1-4)  above  can  be  represented  as: 

V  +r.  "  8?  ^  f (ri>  f(rj)  (1~16) 

i    J 

where  r  is  assumed  to  be  the  smaller  of  the  two  pore  sizes  r .  and  r. . 

It  follows  that  the  total  quantity  of  flow  through  a  unit  volume  of 

porous  medium  is: 


YSnv^2 


q-gZ  I  r  f(r  )  f(r .)  (1-17) 

i  j 


The  discharge  velocity  v  is: 


n  n 
v 


£l^2  f(r±)  f(r.)  (1-18) 

ij  J 

and  using  Darcy's   law  to   solve  for  the  permeability,   k,   yields: 

v  n  n     0 
k-gjllr     f(r.)   f(rj)  (1-19) 


i  J 
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Rewriting  equation   (1-19)    in  terms  of  pore  diameter  and  specific 
permeability: 


1^~2 


or 


K-  3J  H  d     f(d  )   f(d  )  (1-20) 

i  J 


K  =   C     d2  (1-21) 

s 


where 


d2  -    II   d2    f 


I   I   «T    f(d.)f(d.) 


n  is  the  porosity,  and 
C  is  the  shape  factor. 

The  value  of  d  will  be  referred  to  as  the  Marshall  model  pore  size 

parameter. 

The  third  model  examined  was  the  Kozeny  hydraulic  radius  model. 

Scheidegger  (1957)  and  Leonards  (1962)  discuss  its  derivation,  and 

Scheidegger  (Chapter  6,  1957)  presents  a  critique  of  the  model.   From 

the  Hagen-Poiseuille  equation,  the  seepage  velocity  v  through  the  pore 

space  is  assumed  to  be: 


51  =  v_  =  C_  *_§.  r^  (i_23) 


s  \i 

where  a  is  the  cross-sectional  area  of  flow  and 
R^  is  the  hydraulic  radius. 

The  discharge  velocity  v  is: 
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v  -  Cs  XI  4  n  (1-24) 


s   |4 
Applying  Darcy's  law,  the  permeability  is  equal  to 


k  =  C  *  rJ  n  (1-25) 

s  y  h 


The  specific  permeability  is  then  equal  to: 


K  =  Cg  R^  n  (1-26) 


The  hydraulic  radius  R^  can  be  expressed  as  (Leonards,  1962): 

"h  -  <r  (1-27) 


where  S  is  the  surface  area  per  unit  volume  or  the  specific  surface, 
s 

Section  1-2.3  will  describe  a  method  of  calculating  the  specific 
surface  S  from  the  pore  size  distribution  (Rootare  and  Prenzlow, 
1967)  as  follows: 


f(d.) 

S  =  4  n  7  — ■—-  (1-28) 

s       v   d . 
l    i 


where  £  f(d±)  =  1. 

The  hydraulic  radius  R^  may  then  be  expressed  as  a  pore  size  parameter 

as  follows: 


\  "  TIE)  cl-29> 

♦J-sr 

x         i 
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and  the  specific  permeability  as: 

2 


K  =  C  n   .  \(A   ,  (1-30) 

s     1  f(d.) 

4^ 


Summarizing,  three  different  models  have  been  described  which 
involve  the  calculation  of  permeability  from  pore  size  distribution. 
The  pore  size  parameters  from  each  of  the  models  will  be  correlated 
with  actual  permeability  measurements  for  the  soils  tested  to  gain 
some  insight  into  the  relationship  between  the  permeability  and  pore 
size  distribution  of  fine-grained  soils.   Table  1  presents  a  summary 
of  the  models. 

1-2  Pore  Size  Distribution 

The  measurement  of  pore  size  distribution  is  one  technique  of 
inferring  the  relative  arrangement  of  particles  and  pores  (i.e.,  the 
fabric)  of  a  soil.   Frequently,  knowledge  of  the  gross  soil  properties 
is  insufficient  information  to  understand  or  predict  the  behavior  of  a 
soil  when  it  is  subjected  to  changed  conditions.   Information  about  the 
pore  size  distribution  and  changes  in  distribution  with  changed 
conditions  can  be  useful  in  interpreting  soil  behavior.   Pore  size 
information  may  also  be  helpful  in  obtaining  information  about  the 
depositional  history  of  a  soil,  the  movement  of  fluids  and  solids 
through  a  soil,  and  the  chemical  reactivity  of  a  soil.   Compression, 
sample  disturbance  or  remolding,  structural  collapse  and  swelling  can 
all  be  detected  by  determining  changes  in  pore  size  distribution. 

This  writer  believes  that  the  measurement  of  pore  size  distri- 
bution has  been  refined  to  such  a  stage  that  it  may  become  a  routine 
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TABLE  1 


SUMMARY  OF  PERMEABILITY  MODELS  USED 


1)   VARIABLE  DIAMETER  CAPILLARY  MODEL 
Equation 


K  =  Cg  E(d  )  n 


Pore  Size  Distribution 
Parameter 


2)   MARSHALL  MODEL 
Equation 


E(d2)  -  I   f(d.)  d2. 
i 

where  £  f(d.)  =  1 


K  =  C  d 
s 


Pore  Size  Distribution 
Parameter 


-2   £  ?  -2 


3)   HYDRAULIC  RADIUS  MODEL 
Equation 

Pore  Size  Distribution 
Parameter 


d  =  H   dZ  f(d.)  f(d  ) 
i  J  J 

where  d  =  d.  if  d.  <  d. 
i     i    J 

or  d  =  d,  if  d±  >   d 

and  I   f(d±)  -  J  f(dj  ■  n 


K  -  Cs  h  n 


^1     1  f(d.) 

where  £  f(d.)  =  1 
i    X 


K  =  specific  permeability,  L 
n  =  porosity 


C  =  shape  factor 


d.    =  pore  diameter,  L 

f(d.)  =  volumetric  frequency  of  d. 
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test  for  research  studies  as  well  as  for  special  practical  projects 
where  knowledge  of  the  soil  fabric  is  required. 

1-2.1  Applications 
Pore  size  distribution  studies  have  found  application  in  several 
fields  including  chemical  engineering,  petroleum  engineering,  materials 
engineering,  soil  science  and  geotechnical  engineering.   Ahmed  (1971) 
presents  a  thorough  review  of  the  application  of  pore  size  distribution 
measurements  primarily  in  the  geotechnical  field.   Several  of  the  pore 
size  distribution  studies  not  discussed  by  Ahmed  (1971)  will  be 
mentioned  here. 

Silveira  (1965)  combined  pore  size  distribution  calculations  with 
probabilistic  modeling  to  study  the  problem  of  fine  particles  passing 
through  protective  sand  filters.   The  pore  size  distribution  is 
calculated  from  the  grain  size  distribution  by  finding  the  probability 
of  the  various  particle  sizes  occurring  in  assemblages  of  three 
spherical  particles  tangent  to  each  other.   From  the  assemblages, 
equivalent  pore  diameters  are  calculated  and  a  pore  size  distribution 
curve  is  generated.   With  the  pore  size  distribution  calculated  for 
the  filter  material,  the  absorbing  Markov  chain  process  is  used  to 
calculate  the  mean  distances  which  a  given  soil  particle  will  travel 
through  the  filter.   This  ingenious  approach  can  be  used  to  determine 
the  suitability  of  a  filter  material  and  the  thickness  of  filter 
required  to  prevent  clogging. 

Sridharan  (1968)  successfully  applied  pore  size  distribution 
studies  to  calculate  the  negative  pore  pressures  in  partly  saturated 
clays.  By  measuring  the  surface  area  of  a  soil  and  its  pore  size 
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distribution,  Sridharan  was  able  to  estimate  the  negative  pore  water 
pressures  resulting  from  adsorption  and  capillary  retention  of  water 
for  compacted  and  artificially  sedimented  soils.   The  objective  of  his 
research  was  to  gain  a  better  understanding  of  the  shear  strength  of 
partly  saturated  soils  by  examining  the  soil  structure. 

Badger  and  Lohnes  (1973)  used  mercury  intrusion  pore  size  distri- 
bution measurements  to  study  the  fabric  of  natural  and  compacted 
loess.   Samples  of  loess  compacted  in  the  laboratory  to  void  ratios 
equal  to  that  of  natural  loess  were  found  to  have  fewer  large  pores  than 
the  undisturbed  loess.   They  also  identified  several  types  of  loess 
soil  structure  by  comparing  the  grain  size  and  pore  size  distributions. 

Lohnes,  Tuncer  and  Demirel  (1976)  studied  the  affect  of  pre- 
cipitation on  the  structure  of  tropically  weathered  basaltic  soils. 
Pore  size  distribution  measurements  with  mercury  intrusion  revealed 
that  increased  precipitation  resulted  in  a  finer  pore  fabric  for  the 
residual  soils  tested.   Void  ratio  measurements  did  not  correlate  with 
these  changes  in  pore  size  distribution. 

Reed  (1977)  demonstrated  the  relation  between  frost  heave  and  pore 
size  distribution  of  compacted  silty  soils.   Conventionally,  the  frost 
susceptibility  of  a  soil  is  predicted  from  the  grain  size  distribution 
and  soil  texture,  and  no  consideration  is  given  to  soil  fabric.   Reed's 
work  demonstrated  that  samples  compacted  dry  of  optimum  display 
significantly  greater  frost  heave  than  those  compacted  at  optimum  or 
wet  of  optimum.   He  also  determined  an  empirical  equation  to  predict 
frost  heave  from  pore  size  distribution  parameters  of  the  soils  tested. 
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Bhasin  (1975)  presented  a  critical  review  of  the  various  techniques 
available  to  measure  pore  size  distribution  of  silty-clayey  soils.  He 
concluded  that  only  the  mercury  intrusion  method  is  capable  of  measuring 
the  entire  range  of  pore  sizes,  up  to  five  orders  of  magnitude,  of  fine 
grained  soils. 

1-2.2  Theory,  Assumptions  and  Limitations 
Pore  size  distribution  determinations  for  this  study  were  performed 
with  mercury  intrusion.   The  mercury  intrusion  equipment  employed  has  a 
pressuring  capacity  of  15,000  psi  and  is  capable  of  intruding  pores  from 
600  pm  to  0.016  um  in  size.   New  commerical  equipment  currently  available 
has  extended  the  lower  pore  size  limit  to  0.004  ]M   by  increasing  the 
pressure  generating  capacity  to  60,000  psi. 

The  mercury  intrusion  technique  is  based  on  the  principle  that  the 
surface  tension  of  a  non-wetting  liquid  (one  which  has  a  contact  angle 
greater  than  90  with  a  given  solid)  will  oppose  the  entry  of  the 
liquid  into  a  small  pore  of  a  solid.   Washburn  (1921)  determined  that 
this  opposition  could  be  overcome  by  external  pressure,  and  that  the 
external  pressure  required  was  inversely  proportional  to  the  pore 
diameter.  Assuming  a  cylindrical  pore,  Washburn  calculated  the 
relation  to  be: 

4  T  cos  9 
P  =  -j (2-1) 


where  P  is  the  absolute  pressure  required 

T  is  the  surface  tensii 
g 

8  is  the  contact  angle 
d  is  the  pore  diameter 


T  is  the  surface  tension  of  the  intruding  liquid 

6  is  the  contact  angle  between  the  solid  and  the  liquid 
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Mercury  was  used  as  the  intruding  fluid  because,  according  to 
Wintlow  and  Diamond  (1970),  "...  its  practical  advantages  include  low 
vapor  pressure,  relative  inertness  in  terms  of  chemical  reactivity,  and 
the  fact  that  it  is  normally  non-wetting  for  most  kinds  of  surfaces." 
The  surface  tension  of  mercury  used  for  this  study  was  that  determined 
by  Kemball  (1946)  and  recommended  by  Diamond  (1970),  of  484  dynes/cm 
at  25°C. 

Diamond  (1970)  demonstrated  that  the  contact  angle  between  mercury 
and  clay  minerals  is  dependent  on  the  mineral  present.   He  measured, 
by  the  sessile  drop  method,  the  contact  angle  for  mercury  on  kaolinite 
or  illite  which  was  within  one  degree  of  147  ;  and  mercury  on  mont- 
morillonite  which  was  within  one  degree  of  139  .   A  contact  angle  of 
147  was  used  for  this  study  since  kaolinite  was  the  principal  clay 
mineral  present  in  the  soils  tested. 

Briefly,  pore  size  distribution  is  determined  in  the  following 
manner.   The  sample  is  initially  evacuated  and  surrounded  by  mercury, 
the  pressure  is  raised  in  small  increments  and  the  volume  of  mercury 
entering  the  sample  after  each  increment  is  recorded.   Each  pressure 
increment  forces  mercury  into  the  accessible  soil  pores  of  a  diameter 
larger  than  or  equal  to  that  calculated  by  the  Washburn  equation  for 
the  given  pressure.   In  this  manner,  the  volume  of  pore  space  between 
pressure  increments,  and  thus  diameter  decrements,  is  recorded, 
generating  a  pore  size  distribution. 

The  assumptions  and  sources  of  error  of  the  mercury  intrusion 
technique  are  discussed  in  detail  by  Ritter  and  Drake  (1945),  Rootare 
(1968) ,  Winslow  (1969) ,  and  Diamond  (1970) ,  but  it  is  worth  reviewing 
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the  more  important  aspects.   The  most  important  consideration  is  that 
the  pore  diameter  measured  at  a  given  pressure  is  the  diameter  of  the 
opening  of  the  pore.   "If  there  exist  in  the  material  pores  which  are 
considerably  larger  than  the  largest  entrances  to  them,  these  pores 
will  be  measured  as  the  size  of  the  largest  opening"  (Ritter  and 
Drake,  1945).   For  this  reason,  the  pore  diameters  measured  are 
referred  to  as  limiting  pore  diameters;  the  "true"  pore  size  distri- 
bution curve  could  contain  more  pore  volume  in  the  larger  diameters  than 
the  pore  size  distribution  measured. 

Other  assumptions  or  factors  to  be  considered  for  mercury  intrusion 
are: 

1)  Only  pores  accessible  from  the  outside  of  the  sample  are 
penetrated. 

2)  The  minimum  pore  size  measured  is  dependent  on  the  pressure 
capacity  of  the  porosimeter. 

3)  The  soil  is  considered  to  be  incompressible  during 
pressurization  and  does  not  break  down. 

A)  The  contact  angle  and  surface  tension  values  are  assumed 
constant  during  pressurization. 

5)  The  mercury  will  not  react  with  the  soil  during  intrusion. 

6)  Cylindrical  pore  shapes  are  assumed  to  calculate  pore 
diameter. 

In  Geotechnical  Engineering,  the  anisotropy  of  the  fabric  is 
frequently  of  interest.   Unfortunately,  the  mercury  intrusion  technique 
is  of  little  value  in  this  respect  because  the  intrusion  process  occurs 
in  all  directions  which  have  large  enough  pore  openings  for  the 
mercury  to  enter. 
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As  Winslow  (1969),  Diamond  (1970),  and  Bhasin  (1975)  all  pointed 
out,  the  proper  determination  and  calculation  of  pore  size  distributions 
is  dependent  on  several  precautions  and  corrections  which  must  be 
considered: 

1)  During  the  intrusion  procedure,  enough  time  must  be  allowed 
for  the  mercury  to  reach  an  equilibrium  intrusion  for  each 
pressure  increment,  prior  to  measuring  the  incremental 
volume  of  intrusion. 

2)  After  the  initial  filling  of  the  sample  container  with 
mercury,  there  remains  a  small  volume  of  air  trapped  outside 
the  sample.   For  pressures  less  than  one  atmosphere,  a 
Boyle's  law  correction  must  be  made  to  account  for  the 
decreasing  volume  of  air  as  the  pressure  increases. 

3)  Another  volume  correction  must  be  included  for  the 
compression  of  mercury  and  the  expansion  of  the  sample 
container  (the  glass  penetrometer).   This  volume 
correction  is  only  significant  for  pressures  higher  than 
100  psi,  and  should  be  determined  experimentally  by  a 
mercury  intrusion  run  without  a  sample. 

1-2.3  Presentation  of  Data 
Two  methods  have  been  used  to  graphically  present  pore  size 
distributions:  the  differential  distribution  curve  and  the 
cumulative  distribution  curve.   Both  are  useful  for  interpreting  the 
shape  and  changes  in  shape  of  pore  size  distributions.   The 
differential  distribution  shows: 
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1)  the  type  of  distribution, 

2)  the  most  frequent  or  modal  pore  sizes, 

3)  the  symmetry  of  the  distribution, 

4)  any  gaps  or  irregularities  in  the  data,  and 

5)  permits  comparison  of  various  curves. 

The  cumulative  distribution  displays  the  fraction  of  pore  space  greater 

or  smaller  than  a  given  size,  and  the  fraction  of  pore  space  between 

any  two  unmeasured  pore  sizes  can  be  quickly  evaluated.   Ideally,  to 

fully  illustrate  the  pore  size  distribution  of  a  material,  both  the 

differential  and  the  cumulative  pore  size  curves  should  be  presented. 

The  differential  distribution  curve  consists  of  some  differential 

pore  volume  function  plotted  against  its  corresponding  pore  diameter. 

Drake  and  Ritter  (1945)  were  the  first  to  present  such  curves.   The 

method  they  proposed  has  been  found  to  be  suitable  when  the  range  of 

pore  sizes  extend  over  one  or  two  orders  of  magnitude.   For  fine-grained 

soils,  however,  the  range  of  pore  sizes  may  extend  across  three  to  five 

orders  of  magnitude.   As  Winslow  and  Diamond  (1971)  stated, 

"...the  differential  presentation  is  severely  distorted 
if  the  range  of  diameters  covers  several  orders  of  magnitude; 
it  vastly  overemphasizes  the  apparent  importance  of  the  finest 
pore  sizes  at  the  expense  of  the  coarser  sizes." 

An  alternate  method  of  determining  the  differential  distribution  which 

may  be  free  of  the  distortion  problems  is  presented  in  Section  1-2.4. 

The  cumulative  distribution  curve  consists  of  a  plot  of  a 

cumulative  void  volume  parameter  versus  its  corresponding  pore 

diameter.   The  cumulative  void  volume  parameter  commonly  used  is  the 

volume  of  voids  per  dry  weight  of  solid.   For  this  study,  however,  the 

cumulative  volume  of  voids  per  volume  of  sample,  or  cumulative 
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porosity,  was  used  because  of  its  familiarity  to  geotechnical  engineers. 
With  this  exception,  the  cumulative  pore  size  distribution  curves 
presented  in  this  study  are  presented  in  the  same  fashion  as  those 
shown  by  Diamond  (1970),  Sridharan  et  al.  (1971),  Ahmed  et  al. (1974) 
and  Bhasin  (1975). 

If  pore  size  distribution  measurements  are  to  be  correlated  with 
soil  properties  either  to  understand  or  predict  soil  behavior,  it  is 
desirable  to  express  the  pore  size  distribution  in  mathematical  terms. 
Drake  and  Ritter  (1945)  have  used  curve  fitting  techniques  to  quantify 
relatively  simple  pore  size  distributions,  but  for  most  soils  the 
distributions  are  generally  very  complex  and  frequently  multimodal, 
which  makes  curve,  fitting  a  very  arduous  and  impractical  task.   An 
alternative  is  to  determine  descriptive  measures  of  the  pore  size 
distribution  which  reflect  the  shape  of  the  entire  distribution  and 
which  accurately  reflect  changes  in  soil  behavioral  properties. 
Section  1-1.3  examined  three  models  for  predicting  permeability  which 
employ  pore  size  distribution  parameters,  as  shown  on  Table  1. 

The  hydraulic  radius  permeability  model  requires  knowledge  of  the 
surface  area  of  the  porous  medium.   Conveniently,  Rootare  and  Prenzlow 
(1967)  presented  a  method  of  calculating  surface  area  from  pore  size 
distribution  measurements.   The  formula  is  given  as: 


..Jol 


P  dV 

-  £4^  «-2) 
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where  S.A.  is  surface  area 

P  is  the  pressure  applied 

T  is  the  surface  tension  of  the  intruding  liquid 

0  is  the  contact  angle  between  the  intruding  liquid  and  solid 

dV  is  the  differential  volume  of  pores  intruded 

V  is  the  total  volume  of  pores, 
v  r 

Substituting  the  Washburn  equation  (2-1)  into  equation  (2-2)  for  P 
yields: 


s.a.  -  a  £  «£  (2-3) 


where  dV  is  the  differential  volume  of  pores  intruded  between  r  and  dr. 

The  specific  surface  area  S  ,  by  definition,  is  equal  to  the 
surface  area  S.A.  divided  by  the  total  volume  V,  therefore: 


S  -  £  fnV  ^  C2-4) 

s   V  I  0   r 


Factoring  the  volume  of  voids  out  of  the  intergral  results  in: 


2  V 


XdV 
0  r 


S  =  2n 

s 
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where  n  is  the  porosity. 

Expressing  equation  (2-5)  in  a  discrete  form  and  substituting  f(r.), 

the  volumetric  frequency,  for  dV  yields: 


1  f(r.) 

ss  =  2n  I  —r-  C2_6) 

s         i     ri 
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Rewriting  equation  (2-6)  in  terms  of  pore  diameter  d  yields: 


1  f(d.) 
s      i   di 


Theoretically,  this  method  of  determining  the  specific  surface 
area  is  valid  only  if  there  are  no  "ink-bottle  pores"  and  if  the 
pressure  applied  is  sufficient  to  intrude  all  of  the  accessible  pore 
volume.   Neither  of  these  assumptions  were  fully  satisfied  in  this 
study.   As  mentioned  previously,  the  mercury  intrusion  technique 
measures  "limiting  pore  diameters";  consequently,  ink  bottle  pores  are 
measured  as  pores  of  the  neck  diameter  in  size.   Also,  the  intrusion 
pressure  capacity  of  the  equipment  used  was  only  able  to  penetrate 
about  95%  of  the  pore  volume.   These  two  limitations  may  result  in 
some  disparity  if  the  specific  surface  area,  as  calculated  by  equation 
(2-7),  is  compared  to  surface  area  measurements  determined  by  other 
means.   This  does  not,  however,  preclude  the  use  of  the  equation  as  a 
descriptive  parameter  of  the  pore  size  distribution. 

1-2.4  Differential  Distribution  and 
Calculation  of  Parameters 

Soils  are,  basically,  an  agglomeration  of  various  size  particles 
with  a  pore  network  present  where  particles  are  not.   The  pore  system, 
is  extremely  complex  in  geometry;  therefore,  to  speak  of  a  pore  network 
in  terms  of  pore  sizes  and  diameters  is  a  rather  abstract  concept. 

Scheiddeger  (1957),  in  an  effort  to  clarify  the  nebulous  nature  of 
a  pore  diameter  and  its  distribution,  has  considered  a  probabilistic 
approach  to  define  a  pore  system.  He  defines  a  pore  diameter  at  any 
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one  point  within  the  pore  space  as  "...the  diameter  of  the  largest 
sphere  which  contains  this  point  and  remains  wholly  within  the  pore 
space."  In  this  way,  a  pore  "diameter"  can  be  assigned  to  each  point 
in  the  pore  space,  and  a  pore  size  distribution  can  be  defined  by 
determining  what  fraction  f  of  the  pore  space  has  a  pore  diameter 
between  6  and  6  +  d6  (6  denotes  pore  diameter).   Scheiddeger,  in  effect, 
considers  the  pore  diameter  6  to  be  a  continuous  random  variable  with 
some  probability  distribution  of  occurrence  given  by: 


f 


f (6)  d6  =  1  (2-8) 

o 


Although,  Scheiddeger 's  definition  of  a  pore  diameter  is  more 
precise  than  that  measured  by  mercury  intrusion,  the  concept  of  a  pore 
diameter  being  a  continuous  random  variable  is  very  useful  for  analysis 
purposes.   Statistical  methods  and  probability  theory  may  be  used  to 
analyze  pore  size  distributions. 

The  usual  procedure  for  determining  the  density  function  f(x)  of  a 
continuous  random  variable  "x"  is  to  estimate  it  from  limited 
statistical  data.   The  procedure  involves  grouping  the  data  into 
classes  (or  intervals)  of  x  and  forming  a  frequency  distribution,  by 
determining  the  proportion  of  measurements  in  each  of  the  classes.   The 
probability  distribution  of  x  is,  thus,  approximated  by  a  discrete 
frequency  distribution  across  intervals  of  x.   Details  of  the  procedure 
are  presented  in  Walpole  and  Myers  (1972,  pp.  41  to  45)  and  Harr 
(1977). 

The  same  procedure  may  be  applied  to  pore  size  distribution 
measurements  by  mercury  intrusion.   Each  pressure  increment  corresponds 
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to  a  decrease  in  pore  diameter  from  6.  ,  to  6.  which  represents  a  class 

or  interval  of  the  pore  diameter.   Associated  with  each  of  these 

intervals  is  the  fraction  of  porosity  f(d, )  measured  between  6  and  6.  - , 

where  d.  is  the  midpoint  between  6.  and  6j  ,. 
i  l      i-1 

To  approximate  the  differential  distribution  a  frequency  histogram 
is  plotted.   The  histogram  is  conventionally  drawn  with  equal  class 
widths  in  order  to  prevent  distortions  in  the  data.   Since  pore  size 
distributions  extend  across  several  orders  of  magnitude,  plotting  a 
histogram  with  classes  of  constant  width  is  not  practical.   The  frequency 
of  the  smaller  pores  is  overemphasized,  as  previously  discussed  in 
Section  1-2.3. 

To  overcome  these  distortions  in  a  pore  size  frequency  histogram, 
the  class  widths  may  be  set  at  logarithmically  equal  intervals  such 
that : 


log  6i_1  -  log  6±  ■  c  (2-9) 


6i-l 

-~-  =  10  (2-10) 

i 


where  6  is  the  pore  diameter  corresponding  to  an  intrusion  pressure  and 
c  is  the  logarithmic  interval  constant. 

From  equation  (2-1),  it  is  evident  that  the  intrusion  pressure  is 

inversely  proportional  to  the  pore  diameter.   For  each  pressure 

increment  to  correspond  to  a  logarithmically  equal  class  width,  the 

pressure  P  must  be  raised  as  follows: 


P. 

^T—  -  10°  (2-11) 

i-1 
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By  increasing  the  pressure  of  intrusion  according  to  equation  (2-11)  and 
maintaining  c  constant  the  pore  size  data  can  be  presented  in  a  histogram 
as  shown  by  Figure  2 . 

The  differential  pore  size  distribution  is  approximated  from  the 
frequency  histogram  by  connecting  the  values  of  d . ,  f(d. )  as  shown  in 
Figure  2,  where  d.,  the  midpoint  of  each  class,  is  calculated  as  follows: 


di  = 


1 
6i  '  6i-l)2 


and  f(d. )  is  fraction  of  porosity  intruded  between  6  1  and  6..   Such 
a  curve  is  referred  to  as  a  frequency  polygon  (Harr,  1977  pp.  6  and  7). 
All  the  differential  distributions  presented  in  this  study  were 
determined  directly  from  the  frequency  polygon  of  the  pore  size 
measurements  as  shown  in  Figure  2. 

The  pore  size  parameters  used  in  the  permeability  predicting  models, 
shown  on  Table  1,  were  also  calculated  from  d.  and  f(d.).   No  graphical 
intergration  was  necessary.   The  parameters  were  calculated  directly 
from  the  mercury  intrusion  data  by  computer. 

1-3  Soil  Drying  Methods 

The  measurement  of  pore  size  distribution  by  mercury  intrusion 
requires  dry  samples.   The  dehydration  of  soil  samples  without  changing 
the  fabric  presents  a  serious  problem.   Several  methods  have  been 
developed  to  remove  the  water  from  soils.   Various  means  have  also  been 
used  to  check  for  fabric  changes  caused  by  dehydration. 

Before  discussing  the  drying  techniques,  it  would  be  fruitful  to 
mention  how  fabric  changes  may  be  detected.   The  first  requirement  is 
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to  have  an  idea  of  the  kind  of  fabric  change  which  may  take  place  as  a 
result  of  a  particular  drying  technique.   The  dehydration  process  may 
cause  shrinkage,  expansion,  cracking,  fabric  reorientation  or  even 
mineralogical  changes.   Knowing  what  to  expect  will  dictate  what  to 
monitor  before  and  after  water  removal. 

Volume  changes  usually  accompany  fabric  changes,  and  since  volume 
change  is  relatively  easy  to  measure,  it  is  commonly  monitored. 
Diamond  (1970),  Sridharan,  Altschaeffl  and  Diamond  (1971),  Ahmed  (1971), 
Bhasin  (1975),  Zimmie  and  Almaleh  (1976)  and  Reed  (1977)  all  measured  a 
volumetric  characteristic  such  as  total  volume,  porosity  or  void  ratio 
before  and  after  dehydration.   This  should  be  a  requirement  for  all 
mercury  porosimetry  studies.   In  this  way,  any  large  shrinkage  or 
expansion  of  dehydrated  samples  can  be  detected  and  the  efficacy  of 
the  drying  procedure  can  be  judged. 

Another  technique  for  monitoring  fabric  changes  during  dehydration 
may  be  used  to  contrast  separate  dehydration  methods.   This  method 
involves  comparing  the  pore  size  distribution  curves  of  soils  prepared 
in  identical  fashion  but  dried  by  different  techniques.   Diamond 
(1970),  Sridharan  et  al.  (1971)  and  Ahmed,  Lovell  and  Diamond  (1974) 
present  such  comparisons.   These  types  of  comparisons  are  useful 
because  they  show  both  the  total  volumetric  difference  and  in  what  pore 
size  range  these  differences  are  taking  place. 

If  mineralogical  change  is  suspect  as  a  result  of  a  dehydration 
process,  more  sophisticated  soil  fabric  and  composition  analysis  may 
be  required.   Differential  thermal  analysis,  X-ray  diffraction  and 
scanning  electron  microscopy  techniques  are  useful  in  determining  what 
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mineralogical  changes  are  being  caused  by  the  dehydration  technique. 
The  applications  and  suitability  of  these  methods  is  discussed  in 
greater  detail  by  Mitchell  (1976)  and  Yong  and  Warkentin  (1975).   These 
methods  have  seldom  been  used  in  pore  size  distribution  studies; 
although,  Bhasin  (1975)  made  extensive  use  of  them  to  check  the 
reliability  of  critical-region  drying. 

The  three  most  frequently  used  methods  of  drying  soils  are  air  or 
oven  drying,  freeze  drying,  and  critical  region  drying.   The  proper 
technique  is  dependent  on  the  soil  to  be  tested  and  the  time  and  costs 
one  is  willing  to  expend. 

1-3.1  Air  and  Oven  Drying 

Simple  air  or  oven  drying  has  been  used  by  Diamond  (1970)  and 
Sridharan  et  al.  (1971).   The  success  of  this  drying  technique  is  largely 
dependent  on  the  degree  of  saturation,  strength,  and  compressibility 
characteristics  of  the  soil.  Air  or  oven  drying  results  in  the 
evaporation  of  water  from  the  soil  through  the  air-water  menisci.   As 
the  water  exits,  meniscus  radii  become  smaller,  increasing  the  surface 
tension,  which  results  in  increasing  compressive  forces  on  the  soil 
fabric.   If  these  compressive  forces  exceed  the  soil  structural 
"strength",  shrinkage  will  result.   Some  clays  may  crack  and  increase 
in  volume  because  of  the  shrinkage  forces.   Yong  and  Warkentin  (1975, 
pp.  197-206)  discuss  the  shrinkage  process  in  greater  detail. 

Sridharan  et  al.  (1971)  and  Ahmed  et  al.  (1974)  concluded  that  clays 
compacted  at  low  initial  degrees  of  saturation  suffered  only  minor 
shrinkage  with  oven  drying.   This  was  not  the  case  for  compacted  clays 
near  or  above  their  optimum  water  content,  however.   Their  explanation 
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for  this  behavior  was  that  the  drier  soil  had  a  greater  strength 
because  of  the  large  negative  pore  pressures.   In  general,  air  or  oven 
drying  is  not  an  effective  drying  method  for  saturated  clayey  soils 
because  of  the  large  shrinkage  that  occurs  with  drying.   Evidence 
gathered  by  Zimmie  et  al,  (1976),  Diamond  (1970)  and  Bhasin  (1975) 
confirm  this. 

1-3.2  Freeze  Drying 

Freeze  drying  is  another  method  of  removing  the  water  from  soils. 
The  method  consists  of  quickly  freezing  a  small  soil  specimen  to 
cryogenic  temperatures,  followed  by  vacuum  drying  the  frozen  sample  to 
remove  the  water  by  sublimation.   This  process  eliminates  the  air-water 
menisci-shrinkage  effect  replacing  it  with  an  air-solid  system.   Ahmed 
(1971)  and  Zimmie  et  al.  (1976)  describe  the  theory  and  process  in 
detail.   However,  the  more  salient  aspects  of  the  process,  and  how  the 
soil  fabric  can  be  disturbed  during  dehydration  are  worth  reviewing. 

Quick  freezing  to  very  low  temperatures  is  done  primarily  for  two 
reasons.   The  rapid  freezing  process  minimizes  ice  crystal  nucleation 
and  water  migration,  thus  preventing  frost  heaving  and  swelling.   There 
is  also  some  evidence  that  freezing  water  to  cryogenic  temperatures 
reduces  the  degree  of  ice  crystallization  (Ahmed,  1971).   By  reducing 
the  degree  of  ice  crystallization,  the  10%  volume  expansion  which 
normally  results  from  freezing  liquid  water  to  ice,  just  below  0  C,  can 
be  reduced.   If  the  freezing  technique  is  not  effective,  significant 
swelling  of  the  soil  sample  can  take  place. 

The  vacuum  drying  process  is  also  critical  to  the  success  of  the 
freeze  drying  technique.   When  the  soil  sample  is  withdrawn  from  the 
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freezing  media,  it  is  quickly  placed  in  a  vacuum  chamber  and  allowed  to 
heat  up.   The  pressure  and  temperature  in  the  chamber  must  be  maintained 
below  the  triple  point  of  water  to  prevent  liquid  water  from  forming 
in  the  sample.   If  liquid  water  should  form,  the  same  mechanism  which 
causes  shrinkage  of  air  dried  specimens  will  result  in  shrinkage  of  the 
freeze  dried  specimens. 

In  effect,  deficiencies  in  the  freeze  drying  method  may  result  in 
either  sample  swelling  or  shrinkage,  depending  on  whether  the  freezing 
process  or  the  drying  process  is  ineffective.   It  is  conceivable  that 
both  swelling  and  shrinkage  may  occur  in  a  freeze  dry  run.   If  such 
were  the  case  and  the  net  volume  change  was  small,  any  significant 
fabric  changes  taking  place  would  go  undetected  by  volume  monitoring 
the  drying  procedure. 

Freeze  drying  of  soils,  according  to  Diamond  (1970),  Tovey  and 
Yan  (1973),  Ahmed  et  al.  (1974)  and  Zimmie  et  al.  (1976),  if  properly 
conducted,  is  far  superior  to  air  or  oven  drying  for  most  clayey  soils. 
Though  some  volume  changes  may  take  place  as  a  result  of  freeze  drying, 
they  will  generally  be  less  than  5%. 

1-3.3  Critical  Region  Drying 
Critical  region  drying  is  the  third  and  most  complex  procedure  of 
dehydrating  soils.   The  method  consists  of  elevating  the  temperature 
and  pressure  of  a  saturated  soil  specimen  in  a  controlled  manner,  to 
the  critical  region  of  water.   In  the  critical  region  both  the  gas 
phase  and  the  liquid  phase  have  the  same  specific  volume  and  thus  exist 
in  a  single  phase.  Once  the  critical  region  is  reached,  the  pressure 
and  temperature  of  the  sample  are  slowly  reduced  in  a  way  that  the 
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water  in  the  sample  remains  in  the  gas  phase.  When  atmospheric 
conditions  are  reached,  the  sample  has  been  "dried"  without  imposing 
the  air-water  menisci-shrinkage  forces  on  the  soil  fabric.   Bhasin 
(1975)  presents  a  more  detailed  explanation  of  the  process. 

Several  problems  may  arise  when  using  the  critical  region  dehydra- 
tion technique.   The  temperature  and  pressure  regulation  schedule  must 
be  gradual  enough  so  that  water  in  the  sample  is  not  subject  to  large 
gradients.   Flow  into  or  out  of  the  sample  as  a  result  of  too  rapid  a 
pressure-temperature  change  can  cause  either  compression  or  swelling, 
depending  on  whether  flow  is  out  of  or  into  the  sample. 

To  reach  the  critical  region  of  water,  the  pressure  must  be  raised 
to  3800  psi  and  the  temperature  to  380°C.   The  alteration  of  the  soil 
composition  at  such  extreme  conditions  must  be  considered.   Bhasin 
(1975),  after  careful  mineralogical  analysis  of  various  natural  and 
artificial  soils,  found  that  no  mineralogic  change  took  place  as  a 
result  of  critical  region  drying.   Bhasin* s  conclusion  that  critical 
region  drying  does  not  alter  the  mineralogic  composition  of  soils  was 
valid  for  the  soils  he  tested.   However,  to  extrapolate  this  conclusion 
for  all  soils  would  not  be  correct. 

This  study  was  originally  designed  to  use  critical  region  drying 
for  dehydrating  the  soil  specimens.   However,  during  preliminary 
critical  region  runs,  void  ratio  increases  of  35  to  40%  were  measured 
after  drying.   Mineralogical  analysis  of  the  soil  before  and  after 
critical  region  drying  revealed  that  the  soil  contained  large  fractions 
of  kaolin  and  dolomite.   During  the  drying  process  the  kaolin  and 
dolomite  reacted  to  form  a  plagioclase  feldspar;  Appendix  C  contains 
further  details  of  this  reaction. 
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Critical  region  drying  is  a  complicated  and  arduous  procedure  in 
its  present  state  of  development,  but  if  oven  drying  and  freeze  drying 
are  not  suitable  for  a  particular  soil,  it  may  be  the  only  alternative 
available.   Based  on  the  experiences  of  this  study,  critical  region 
drying  should  always  be  monitored  with  both  volume  measurements  and 
mineralogical  analysis  to  check  the  efficacy  of  the  drying. 

1-3.4  Fluid  Replacement 
Tovey  and  Yan  (1973)  discuss  fluid  replacement  as  a  possible 
drying  technique.   The  method  involves  impregnating  a  sample  with  an 
organic  liquid  to  replace  the  pore  fluid.   After  replacement,  the 
sample  is  air  or  oven  dried  or  critical  region  dried.   For  air  drying, 
the  replacement  fluid  is  usually  a  liquid  of  low  surface  tension  such 
as  methanol  or  acetone.   This  will  reduce  the  shrinkage  forces  on  the 
soil  during  drying.   For  critical  region  drying,  two  fluid  replacements 
are  necessary.   The  sample  must  initially  be  impregnated  with  acetone, 
then  placed  in  a  high  pressure  cell  and  impregnated  again  with  carbon 
dioxide.   Liquid  carbon  dioxide  is  a  more  suitable  fluid  for  critical 
region  drying  than  water  because  of  its  lower  critical  state  (31.1  C 
and  1050  psi).   The  lower  temperatures  and  pressures  reduce  the 
possibility  of  mineral  alterations. 

Fluid  replacement  is  not  without  its  problems,  however.   Tovey 
(1970)  and  Smart  (1966)  (as  reported  by  Tovey  and  Yan,  1973)  both 
found  that  swelling  frequently  accompanies  fluid  replacement.   Therefore, 
volume  monitoring  would  be  necessary  prior  to  replacement,  after 
replacement,  and  after  drying.   Tovey  and  Yan  (1973)  conclude  that 
fluid  replacement  and  air  or  oven  drying  does  reduce  the  amount  of 
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shrinkage  if  compared  to  air  or  oven  drying  alone,  but  the  method  does 
not  eliminate  shrinkage  completely.   Fluid  replacement  and  critical 
region  drying  may  prove  to  be  fruitful  if  a  suitable  replacement 
technique  is  used.  i 

Of  the  various  dehydration  methods  available,  only  freeze  drying 
and  critical  region  drying  are  suitable  for  most  clayey  soils,  if 
fabric  disturbance  is  to  be  minimized.   Soil  drying  should  always  be 
monitored  with  at  least  volumetric  measurements.   Finally,  regardless 
of  the  drying  method  used,  the  need  for  precise  and  careful  laboratory 
technique  cannot  be  overemphasized.   Sample  trimming,  volumetric 
measurements  and  insufficient  drying  times  can  be  significant  sources 
of  errors. 

1-3.5  Method  Used 

The  drying  technique  used  for  this  study  was  the  freeze  drying 
method  recommended  by  Zimmie  and  Almaleh  (1976).   Liquid  nitrogen  was 
used  as  the  freezing  agent,  and  vacuum  drying  was  performed  at  room 
temperature  without  a  refrigerant.   Zimmie  et  aL  (1976)  measured 
negligible  volume  expansion  of  kaolinitic  soils  after  freezing  in 
liquid  nitrogen  and  felt  that  no  intermediate  cooling  fluids  would  be 
necessary  if  the  samples  were  kept  small  (less  than  8  mm  on  a  side). 
This  recommendation  was  followed  in  this  study. 

Vacuum  drying  at  room  temperature  was  used  because  as  Zimmie  et  aL 
stated:  (1)  "...sublimation  is  an  evaporative  process  accompanied  by 
cooling,  so,  although  the  specimen  container  is  exposed  to  room 
temperature,  the  sample  temperature  will  remain  below  freezing;"  and 
(2)  the  rate  of  sublimation  is  a  function  of  sample  temperature,  thus 
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the  lower  the  temperature,  the  lower  the  rate  of  sublimation.   By 
vacuum  drying  at  room  temperature  without  using  a  refrigerant,  the 
drying  technique  was  simplified  and  accelerated.   Zimmie  et  al. 
determined  that  5  to  6  hours  of  vacuum  drying  would  be  sufficient  for 
kaolinitic  soils  and  9  to  10  hours  would  be  sufficient  for  montmorilloni- 
tic  soils.   To  be  conservative,  8  to  10  hours  of  vacuum  drying  was  felt 
to  be  sufficient  for  this  study.  After  freeze  drying,  the  samples  were 
stored  in  a  dessicator  containing  anhydrous  magnesium  perchlorate  to 
prevent  the  samples  from  adsorbing  moisture. 

The  success  of  the  f reeze-drying  technique  used  in  this  study  was 
strongly  influenced  by  the  size  of  the  samples.   Too  large  a  sample 
resulted  in  cracking  and  large  volume  changes  after  freeze  drying.   The 
method  was  not  100  percent  successful.   A  small  but  significant  percent- 
age of  the  samples  which  were  freeze  dried  had  to  be  discarded  because 
they  cracked  or  crumbled  apart,  as  was  also  reported  by  Zimmie  et  al. 
For  the  successful  runs,  however,  the  maximum  volume  changes  experienced 
as  a  result  of  freeze  drying  was  about  3%,  which  agrees  with  what 
Reed  (1977)  determined. 


47 


2  -  SOILS  STUDIED 

The  soils  were  selected  based  on  certain  desirable  characteristics 
and  within  the  limitations  imposed  by  mercury  intrusion  testing.   The 
soils  are  intended  to  be  representative  of  natural  silty  clays  which 
are  frequently  employed  in  earthwork  construction.   It  was  hoped  that 
the  method  of  compaction  and  the  types  of  soils  used  would  give  a  wide 
range  of  permeabilities.   However,  since  the  affect  of  chemical  changes 
in  the  pore  water  on  soil  permeability  was  beyond  the  scope  of  this 

study,  it  was  necessary  to  use  clays  of  low  activity.  This  limited  the 

—8 
minimum  permeability  measured  to  about  10   cm/sec.   The  pore  size 

distribution  measurements,  because  of  the  small  samples  necessary, 

limited  the  maximum  particle  size  to  a  fine  sand  and  required  as 

homogeneous  a  soil  as  practicable.   For  these  reasons,  the  soils  tested 

were  blends  of  natural  silt  and  commercial  kaolin  clay,  artificially 

prepared  in  controlled  proportions. 

The  silt  was  taken  from  a  natural  loess  deposit,  near  the  Wabash 

River,  in  southwestern  Indiana.   The  soil  was  excavated  from  the  second 

bench  level  on  the  east  side  of  U.  S.  Highway  41,  about  1  mile  south  of 

Patoka,  Indiana.   After  air  drying,  the  loessial  blocks  were  passed 

through  a  No.  40  U.  S.  standard  sieve  to  remove  any  large  particles  and 

to  break  down  the  initial  soil  structure.   X-ray  diffraction  analysis 

revealed  that  the  silt  was  composed  primarily  of  quartz  and  dolomite 

with  small  amounts  of  calcite  and  feldspar  also  present.  Although  the 
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loessial  soil  contained  only  about  1  to  3%  clay  size  particles,  a 
separate  X-ray  diffraction  analysis  was  also  performed  on  the  clay 
fraction.  The  preparation  technique  developed  by  Kinter  and  Diamond 
(1956)  was  used  for  this  analysis.   Illite,  kaolinite  and  montmorillonite 
were  found  to  be  present,  with  a  possible  trace  of  vermiculite.  As 
shown  on  Figure  3,  the  silt  has  a  very  uniform  gradation  with 
approximately  85%  of  the  particles  between  0.07  and  0.01  mm  in  size. 
The  silt  is  non-plastic,  as  shown  in  Table  2. 

Edgar  plastic  kaolin  was  blended  with  the  silt  to  prepare  the  soil 
combinations  to  be  tested.   The  kaolin  is  commercially  processed  to 
remove  material  coarser  than  40  ym  in  size,  dried  in  a  tunnel  dryer  to 
300  to  450  F,  and  then  pulverized.   The  manufacturer  claims  that  99.5% 
of  the  clay  mineral  composition  is  kaolin  with  the  remainder  being 
quartz  or  mica.   The  classification  properties  of  the  kaolin,  as 
determined  by  Bhasin  (1975),  are:  specific  gravity  of  solids  =  2.65, 
liquid  limit  =  59%,  plasticity  index  =  22%,  and  the  clay  fraction 
(%  less  than  2  vim)  =  76%. 

The  three  artificial  blends  of  silt  and  kaolin  were  tested: 

90%  silt  -  10%  kaolin 

70%  silt  -  30%  kaolin 

50%  silt  -  50%  kaolin 

Grain  size  distribution  curves  and  soil  classification  properties  for 

each  of  the  mixes  are  shown  on  Figure  3  and  Table  2,  respectively. 
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3  -  APPARATUS  AND  EXPERIMENTAL  PROCEDURE 

3-1  Soil  Mixing  and  Curing 

All  soil  mixing  for  this  study  was  done  with  a  Patterson-Kelley 
liquid  solid  mechanical  blender.   The  blender,  shown  in  Figure  4, 
consists  of  a  slowly  rotating  V  shell  housing  a  high  speed  rotating  bar. 
The  rotating  bar  contains  several  dispersion  blades  and  two  sprayers. 
Connected  to  the  rotating  bar  is  a  graduated  cylinder  which  serves  as 
the  reservoir  for  the  deionized  water.   A  stopcock  at  the  base  of  the 
cylinder  regulates  the  quantity  and  rate  of  flow  into  the  shell. 

The  mixing  procedure  was  begun  by  placing  the  desired  proportions 
of  air  dried  silt  and  kaolin  into  the  blender.   It  was  found  that  2  kg 
of  soil  was  a  sufficient  amount  of  soil  to  prepare  one  sample.   The  soil 
was  dry  mixed  for  approximately  10  minutes.   Deionized  water  was  then 
added  at  a  rate  of  50  to  75  ml  per  minute  while  the  soil  was  being  mixed. 
Once  the  desired  quantity  of  water  was  added,  the  end  caps  of  the 
blender  were  removed,  and  the  soil  was  scraped  off  the  inside  walls  of 
the  shell.   After  the  scraping,  mixing  was  resumed  for  about  7  minutes, 
followed  again  by  scraping  and  another  7  minutes  mixing  period. 
Following  the  final  mixing,  the  soil  was  removed  from  the  shell  and 
sealed  inside  two  polyethylene  bags.   The  bagged  samples  were  placed  in 
a  high  humidity  barrel  and  cured  for  two  days. 

Experimentation  revealed  that  in  order  to  achieve  the  desired 
water  contents  it  was  necessary  to  add  0.5%  less  water  to  the  air-dry 


52 


FIGURE  4     PATTERSON- KELLEY  BLENDER 
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soils  than  calculated.  This  procedure  gave  total  water  contents  within 
0.25%  of  the  desired  water  contents.   Less  than  the  calculated  amount  of 
water  was  added  because  the  silt  and  kaolin  were  at  some  small  initial 
water  content  in  the  air  dried  state.  , 

3-2  Compaction 

3-2.1  Apparatus 
Kneading  compaction  and  special  sample  preparation  procedures  were 
used  in  order  to  obtain  a  high  degree  of  homogeneity  in  the  prepared 
samples.   An  electrically  driven,  semi-automatic  kneading  compactor, 
manufactured  by  the  August  Company,  was  used  to  prepare  the  samples. 
The  kneading  compactor  is  shown  on  Figure  5.   The  compactor  foot  pressure 
is  controlled  by  a  pneumatic-hydraulic  system  with  a  standard  air 
regulator  and  pressure  gage.   The  triangular  shaped  foot  rotates  60 
degrees  between  each  tamp  providing  uniform  coverage  of  the  sample 
surface.   Two  minor  problems  were  encountered  with  the  kneading  compactor 
during  the  course  of  this  study.  First,  the  foot  pressure  imparted  on 
the  soil  is  not  directly  determinable.   It  is  controlled  by  an  air 
regulator,  and  a  pressure  gage  measures  the  hydraulic-pneumatic  pressure 
in  the  compaction  system.   Gaudette  (1960)  and  Reed  (1977)  have 
determined  calibration  curves  for  converting  gage  pressure  to  the  actual 
foot  pressure  being  applied  to  the  soil;  however,  these  curves  are  not 
unique.   Seasonal  changes  of  temperature  and  humidity  can  cause  a 
significant  change  in  foot  pressure  for  a  given  gage  pressure, 
particularly  at  the  lower  pressures.   In  an  effort  to  overcome  this 
problem,  the  groups  of  samples  of  a  given  soil  type  to  be  compacted  at 
the  same  energy  level,  were  usually  prepared  within  the  same  two  week 
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FIGURE  5    KNEADING  COMPACTOR 
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period.   Thus  a  group  of  samples  was  compacted  at  approximately  the  same 
.  energy  level  and  gage  pressure,  though  the  peak,  foot  pressure  during 
compaction  was  unknown.   It  is  suggested  that  before  future  research 
using  the  kneading  compactor  is  undertaken,  a  simple  and  accurate 
method  of  calibrating  gage  pressure  to  foot  pressure  be  devised,  which 
can  be  repeated  on  a  periodic  basis.  One  other  problem  was  encountered 
when  using  the  kneading  compactor.   If  the  foot  was  not  at  the  bottom 
of  its  stroke  when  a  series  of  tamps  was  begun,  the  first  few  tamps 
gave  exceedingly  high  foot  pressures.   Bailey  (1976),  Gaudette  (1960) 
and  Aughenbaugh  et  al.  (1963),  present  further  details  about  the 
operation  of  the  kneading  compactor. 

3-2.2  Procedure 
All  samples  were  compacted  using  the  same  procedure  with  the 
exception  of  the  varying  energy  levels.   Following  the  two  day  curing 
period,  the  samples  were  compacted  in  a  cylindrical  lucite  mold  4.0 
inches  in  diameter  and  4.89  inches  in  length  which  was  coated  with 
silicon  oil.   A  circular  base  plate  extended  0.31  inches  into  the 
bottom  of  the  mold,  thus  making  the  compacted  sample  4.58  inches  long. 
These  dimensions  are  identical  to  those  specified  by  AASHTO  T  99-70 
(i.e.,  1/30  ft  in  volume).  A  circular  strip  of  filter  paper  4  inches 
in  diameter  was  placed  on  the  base  plate  prior  to  placing  the  soil  in 
the  mold.   The  filter  paper  served  two  purposes:  1)  to  prevent  the 
soil  from  adhering  to  the  base  plate,  and  2)  to  prevent  the  soil  from 
plugging  the  porous  stones  during  permeation.   The  soil  was  placed  in 
the  mold  and  compacted  in  five  layers  of  approximately  equal  thickness. 
Thirty  tamps  of  the  compaction  foot  were  applied  per  layer.  Prior  to 
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the  mechanical  compaction  of  each  layer,  the  layer  below  was  scarified, 
fresh  soil  was  added  and  manually  tamped  three  times  with  a  piston  tamper, 
The  scarification  of  each  layer  prevented  cracking  at  the  contact  be- 
tween the  layers.     The  initial  manual  tamping  was  necessary  to  densify 
the   fresh  soil  just  enough  to  prevent  the  mechanical  foot   from  punching 
through  the  layer.     Once  the  compaction  procedure  was  completed,  the 
collar  was   removed  and  the  soil  sample  was  trimmed  flush  with  the  top 
of  the  lucite  mold.     The  weight  of  the  sample  and  mold  was   then  taken, 
and  the  remaining  soil  trimmings  were  used  for  water  content   determin- 
ations .     After  placing  another  circular  strip  of  filter  paper  k  inches 
in  diameter  on  the  top  of  the  sample,  the  mold  and  sample  were  ready  to 
be  placed  in  the  permeameter  cell. 

Several  difficulties  were  encountered  during  the  preliminary  com- 
paction runs.     Soils  with  water  contents   3  to  k  percent   greater  than 
the  optimum  water  content  could  not  be  compacted  with  the  kneading 
compactor.     The  soil  would  stick  to  the  compaction   foot  during  the 
upward  stroke  thus  preventing  effective  compaction.     Soils  with  a  water 
content  slightly  greater  than  the  optimum  water  content  presented 
another  problem.      During  the  downward  stroke  of  the  compaction   foot, 
the  soil  adjacent  to  the  foot  would  heave  upward.     This  heaving  created 
sufficient   forces  to  unseat  the  mold  from  the  base  plate.     This  prob- 
lem was  overcome  by  building  a  collar  which  could  be   fastened  to  the 
base  plate  assembly.     Figure  6  shows  the  mold,   collar  and  tamper. 
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FIGURE  6     COMPACTION  MOLD, COLLAR  AND  TAMPER 
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3-3  Saturation  and  Permeation 

3-3.1  Apparatus 
As  was  discussed  in  Section  1-1.1,  the  permeability  of  a  soil  is  a 
function  of  the  degree  of  saturation.   To  maximize  the  permeability  and 
minimize  the  effect  of  the  degree  of  saturation,  samples  were  brought  as 
close  to  full  saturation  as  practicable  prior  to  permeation. 

The  permeameter  cells  used  to  house  the  soil  samples  during 
saturation  and  permeation  are  shown  on  Figure  7.   The  cells  consist  of  a 
lucite  tube  and  two  lucite  end  caps  containing  brass  porous  stones.   The 
transparent  lucite  parts  allow  one  to  view  the  soil  during  permeation. 
This  was  very  advantageous  because  changes  in  sample  volume,  cracking, 
and  washing  out  of  soil  particles  could  be  detected.   The  top  end  cap  is 
inset  with  a  porous  stone  and  rests  flat  on  the  sample  top.   The  base 
cap  has  a  1/8  inch  slot  for  the  mold  to  rest  in,  and  the  porous  stone  is 
inset  in  the  bottom  of  the  lucite  tube.   The  drainage  lines  in  both  caps 
are  1/4  inch  in  diameter.   Four  grooves  radiate  outward  from  the  center 
of  the  caps.   The  grooves  reduce  the  constriction  of  flow  between  the 
end  cap  and  the  porous  stone  and  prevent  air  bubbles  from  blocking  the 
drainage  line.   Quicklock  fittings  connect  the  drainage  lines  to  the 
permeameter  or  saturation  devices.   The  end  caps  and  mold  are  fastened 
together  by  four  bolts  and  wing  nuts.   "0"  rings,  inset  in  the  caps, 
are  compressed  against  the  edges  of  the  cylindrical  mold,  providing  a 
pressure  tight  seal. 

The  soil  was  compacted  directly  in  the  lucite  mold  to  reduce 
seepage  along  the  walls.  The  permeameter  cell  was  designed  so  the  mold 
fit  directly  into  the  two  caps  without  compressing  or  shifting  the 
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FIGURE  7     PERMEAMETER  CELLS 
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sample  within  the  lucite  tube.  Once  the  sample  was  placed  in  the 
perraeameter  cell,  it  was  confined  to  prevent  an  increase  in  volume. 

The  saturation  device,  built  by  Reed  (1977),  consists  of  two 
reservoirs,  two  air  pressure  regulators  and  one  air  pressure  gage. 
The  reservoirs  could  be  pressurized  independently  or  together,  and  a 
vacuum  system  could  be  connected  to  the  outflow  reservoir. 

The  permeameters  used  are  closed  system,  falling  head  permeameters 
with  provisions  for  back  pressuring.   Because  of  the  wide  range  of 
permeabilities  to  be  measured,  two  permeameters  were  built.   The  two 
permeameters  are  nearly  identical  except  for  the  inside  diameters  of 
the  standpipes.   The  permeameter  used  for  measuring  the  more  permeable 
soils  has  1/2  inch  I.D.  standpipes,  while  the  other  has  1/4  inch  I.D. 
standpipes.   The  standpipe  reservoirs  could  be  pressurized  simultaneously 
from  one  regulator  or  independently  from  two  regulators.   A  vacuum 
could  be  applied  to  one  or  both  reservoirs.  Meter  sticks  were  used  to 
measure  the  water  elevation  heads.   Connected  to  both  permeameters  was 
one  common  reservoir  containing  deionized,  distilled  water.  The 
permeameters  were  designed  to  be  used  for  both  saturation  and  permeation, 
but  in  order  to  accelerate  the  testing  program,  the  saturation  device, 
described  above,  was  used  for  saturation.   Figure  8  shows  the  saturation 
device  and  permeameters  used. 

3-3.2  Procedure 

Following  compaction,  the  lucite  tube  holding  the  sample  was 

connected  to  the  two  end  caps.   The  bottom  cap  was  filled  with  distilled 

water  and  held  a  saturated  porous  stone.   The  top  cap  and  porous  stone 

were  attached  dry.  The  assembled  permeameter  cell  was  then  connected 
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FIGURE  8     SATURATION  DEVICE    (RIGHT)  AND 
PERMEAMETERS  (LEFT) 
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to  the  saturation  apparatus.  A  vacuum  of  about  5  psi  absolute  pressure 
was  placed  on  the  top  of  the  sample  while  the  bottom  end  was  under 
approximately  18  inches  of  water  head.   This  condition  was  maintained 
for  about  8  hours  permitting  deionized,  deaired  water  to  flow  upward 
through  the  sample  and  allowing  air  bubbles  to  escape.  Vacuum 
saturation  of  samples  was  found  to  be  very  effective  for  samples 
compacted  on  the  dry  side  of  optimum.   Simply  passing  water  through  dry 
side  samples  under  a  low  head  difference  seldom  gave  degrees  of 
saturation  greater  than  88%;  although  vacuum  saturation  consistently 
resulted  in  saturation  levels  greater  than  95%. 

After  vacuum  saturation,  back  pressuring  was  begun.   The  pressures 
on  the  two  ends  of  the  sample  were  raised  in  5  psi  increments  with  the 
pressure  at  the  base  of  the  sample  3  to  5  psi  greater  than  the  pressure 
at  the  top.   This  was  done  with  the  hope  that,  as  the  back  pressuring 
increased,  a  fraction  of  the  air  remaining  in  the  sample  would  go  into 
solution  and  flow  out  the  top  of  the  sample.   Fifteen  to  30  minutes 
were  allowed  to  pass  between  the  application  of  pressure  increments, 
when  the  bottom  end  pressure  reached  50  psi,  both  ends  of  the  sample 
were  connected  to  a  common  air  regulator  maintaining  a  50  psi  back 
pressure  overnight. 

Preliminary  testing  to  establish  a  reliable  saturation  technique 
revealed  that  it  was  important  not  to  exceed  a  5  psi  pressure 
differential  across  the  sample.   In  preliminary  tests,  a  10  to  20  psi 
pressure  differential  was  used,  and  the  sample  could  be  seen 
consolidating  in  the  permeameter  cell.  A  gap  of  about  1/4  inch  was 
visible  across  the  base  of  the  sample.  When  the  pressure  differential 
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was  eliminated,  the  gap  would  fill  in,  but  large  horizontal  cracks 
extending  completely  across  the  sample  would  appear.   Such  samples  were 
obviously  unsatisfactory  for  permeability  testing  but  would  not  have 
been  detected  if  the  samples  were  not  visible.  Appendix  B  provides 
further  details  concerning  seepage  forces  and  consolidation. 

With  the  samples  at  50  psi  back  pressure  and  the  saturation 
process  complete,  the  permeability  tests  were  then  conducted.   The 
quicklock  connections  on  the  permeameter  cell  permitted  the  easy 
transfer  of  the  cell  from  the  saturation  apparatus  to  the  permeameter 
while  maintaining  the  high  pressure  in  the  cell.  A  one  meter  head 
difference  in  the  water  columns  was  established  between  the  two 
standpipes,  and  a  50  psi  back  pressure  was  placed  on  both  of  the 
columns.   The  cell  was  then  connected  to  the  permeameter  and  permeation 
was  begun.   The  direction  of  flow  as  from  top  to  bottom  across  the 
sample.   Figure  59  in  Appendix  B  shows  a  sketch  of  the  permeability 
test.   Periodic  readings  of  the  elevation  heads  in  the  columns,  the 
time  of  the  readings,  and  the  air  temperature  were  taken.   The 
frequency  of  the  readings  depended  on  the  rate  of  flow.   It  was 
assumed  that  the  air  temperature  was  equal  to  the  water  temperature 
to  normalize  the  permeability  readings  to  20  C.   The  permeability  tests 
varied  in  duration  from  one  to  six  days,  depending  on  the  permeability 
of  the  soil.   At  least  ten  to  twelve  readings  were  taken  for  each 
sample.   Some  samples  were  repermeated  to  check  the  reproducibility  of 
the  test.   This  was  done  by  disconnecting  the  permeameter  cell  from 
the  permeameter  apparatus,  reestablishing  the  one  meter  head  difference 
between  the  standpipes,  and  reconnecting  the  cell  to  the  permeameter. 
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Following  permeation,  the  samples  were  dissected  to  check  for 
uniformity,  and  trimmed  for  freeze  drying  and  pores  size  distribution 
measurements.   Because  of  this,  it  was  important  to  prevent  sample 
disturbance  and  volume  changes  when  removing  the  permeated  soil  from 
the  lucite  mold.  At  the  completion  of  permeation,  the  sample  was 
gradually  depressurized  from  50  psi  in  5  psi  decrements  per  15  (or 
more)  minutes.   When  atmospheric  pressure  was  reached,  the  sample  was 
kept  in  the  mold  (at  constant  volume)  for  at  least  2  hours.   If  the 
back  pressure  on  the  sample  was  released  too  rapidly,  the  sample  would 
swell  noticeably  when  the  end  caps  of  the  permeameter  were  removed. 

3-4  Freeze  Drying 

Freeze  drying  was  used  to  dehydrate  samples  prior  to  mercury 
intrusion.  The  procedure  used  by  Zimmie  et  al  (1976)  and  Reed  (1977) 
was  used  for  this  study.   The  effectiveness  of  the  drying  procedure  was 
checked  by  the  volumetric  method  discussed  in  Section  1-3. 

3-4.1   Apparatus 
The  apparatus  used  in  this  study  was  assembled  by  Reed  (1977)  and 
this  description  is  taken  directly  from  his  work  (p.  44). 

Figures  9  and  10  show  the  basic  components  of  the  freeze  drying 
equipment  which  include: 

"1)  Wire  cage  to  hold  the  samples  during  freezing  and 

sublimation. 
2)   Large  Dewar  flask  containing  liquid  nitrogen  into 
which  the  sample  and  cage  are  dipped. 
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3)  Dessicator  with  inside  wire  support  from  which 
sample  cage  is  suspended  during  sublimation. 

4)  Vacuum  pump,  capable  of  evacuating  to  less  than 
0.01  mm  of  mercury. 

5)  Condenser,  which  is  hooked  between  the  dessicator  and 
vacuum  pump  in  order  to  prevent  vapor  from  the  soil 
entering  the  pump.  Liquid  nitrogen  contained  in  a 
Dewar  flask  was  used  to  cool  the  condenser." 

3_4.2   Procedure 

After  depressurizing  the  permeated  soil,  the  samples  were  ready 
for  dissection  and  freeze  drying.   The  permeameter  cell  was 
disassembled,  and  the  mold  and  soil  were  weighed  to  determine  the 
final  degree  of  saturation.  The  soil  was  then  extruded  by  an  hydraulic 
extruder.   Quickly  and  carefully,  small  cubes  were  trimmed  from  the 
•  compacted  sample  for  freeze  drying.   The  small  specimens  were  trimmed 
with  a  razor  blade  and  were  approximately  8  mm  long  on  each  side.   Four 
or  five  samples  were  trimmed.   Two  or  three  were  taken  from  about  1  inch 
below  the  top  of  the  sample,  and  the  remaining  ones  were  taken  from  1 
inch  above  the  sample  bottom.   The  trimmed  samples  were  placed  in  the 
cage  and  lowered  into  liquid  nitrogen.   They  remained  in  liquid 
nitrogen  for  approximately  ten  minutes.  While  the  samples  were  immersed 
in  the  liquid  nitrogen,  bubbles  would  form  and  adhere  to  the  sides  of 
the  soil  cubes.   Although  it  was  not  done  in  this  study,  it  would  have 
been  beneficial  to  vibrate  the  specimens  while  they  were  immersed  to 
prevent  insulation  from  the  cold  source. 

While  the  small  specimens  were  being  frozen,  the  remaining  soil 
was  utilized  to  determine  the  water  content  and  porosity  within  the 


68 

compacted  sample.   Samples  from  the  top,  middle  and  bottom  of  the 
remaining  soil  were  trimmed  into  cuhes  3  to  4  cm  on  a  side.   Wax 
displacement,  as  described  in  U.  S.  Dept.  of  the  Army,  EM  1110-2-1906 
(1970),  was  used  to  measure  the  porosity.   The  results  of  these  measure- 
ments were  averaged  and  compared  to  both  the  gross  sample  porosity  after 
compaction,  and  the  porosity  measured  after  freeze  drying.   The  post 
freeze  dry  porosity  was  measured  by  mercury  displacement  in  the 
porosimetry  run.   Thus,  three  different  porosity  measurements  were  made 
to  check  the  effectiveness  of  freeze  drying;  one  of  the  gross  compacted 
sample,  another  of  the  dissected  sample  after  permeation,  and  one  of  the 
freeze  dried  sample. 

After  the  small  samples  had  been  in  the  liquid  nitrogen  the  proper 
length  of  time,  they  were  quickly  transferred,  in  the  wire  cage,  to  the 
dessicator  for  vacuum  drying.   The  samples  were  vacuum  dried  for  eight 
to  ten  hours  then  stored  in  a  dessicator  containing  anhydrous  magnesium 
perchlorate  to  remove  any  remaining  moisture.   The  samples  remained  in 
the  desiccator,  which  was  kept  under  a  slight  vacuum,  until  they  were 
used  for  pore  size  determinations. 

Reed  (1977)  explains  the  importance  of  the  more  salient  procedural 
details.   The  most  important  detail  was  to  trim  the  freeze  dry  samples 
no  larger  than  specified  above  because  the  larger  samples  often  cracked 
and  experienced  large  volume  changes. 

3-5  Pore  Size  Distribution 

The  apparatus  and  procedure  used  in  this  study  to  measure  pore  size 
distribution  by  mercury  intrusion  is  basically  the  same  as  that  employed 
by  Ahmed  (1971),  Bhasin  (1975)  and  Reed  (1977).  For  this  reason  the 
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descriptions  and  details  which  follow  are  kept  brief.   The  references 
listed  above  should  be  consulted  if  further  details  are  desired. 

3-5.1  Apparatus 

Three  primary  components  are  required  to  measure  pore  size  distri- 
bution by  mercury  intrusion:  a  penetrometer,  a  filling  device  and  a 
porosimeter.   The  penetrometer  is  a  glass  capillary  stem  with  a  bulb 
on  one  end  which  is  used  to  house  the  sample.   The  capillary  stem  is 
calibrated  in  0.002  ml  increments  and  has  a  capacity  of  0.20  ml.   A 
steel  cap  and  teflon  fastener  fit  over  the  bulb  which,  when  greased  and 
sealed,  maintains  a  pressure  tight  seal.   The  filling  device  is  a  two 
chambered  glass  tube  which  holds  the  penetrometer  during  filling  and 
low  pressure  intrusion.   The  filling  device  is  connected  to  a  pressure 
control  board  which,  in  turn,  is  connected  to  a  vacuum  pump,  McLeod 
gauge,  mercury  manometer  and  bleeder  valve.   The  vacuum  pump  is  for 
pumping  down  the  pressure  in  the  filling  device  to  the  required  vacuum 
prior  to  mercury  filling.   The  McLeod  gauge  is  for  measuring  pressures 
less  than  1  mm  of  mercury,  while  the  manometer  is  used  to  measure 
pressures  greater  than  1  mm  of  mercury  but  less  than  atmospheric 
pressure. 

Pressures  greater  than  atmospheric  are  applied  by  an  Aminco- 
Winslow  Porosimeter.   The  porosimeter  has  a  pressure  capacity  of  15,000 
psi  and  an  electrical-mechanical  sensing  device  which  measures 
volumetric  increments  of  0.0001  ml  in  the  penetrometer  stem.   Intrusion 
pressures  are  measured  by  two  Bourdon  pressure  gauges  on  the 
porosimeter:  one  which  has  a  capacity  of  1,000  psi  and  the  other  a 
capacity  15,000  psi. 
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3-5.2  Procedure 
Following  freeze  drying,  the  samples  to  be  lntruded  ^  ^^ 
with  a  hammer  and  chisel  along  several  faces  then  bio™  clean  with  a 
dry  air  stream   The  tricing  method  was  intended  to  cause  a  tensile 
breafc  in  the  soil  to  prevent  surfacfi  pores  from  ^  ^  ^ 

The  sample  was  weighed,  then  sealed  in  the  penetrometer  and  weighed 
again.   The  penetrometer  and  sample  were  placed  in  the  filling  chamber 
and  evacuated  to  0.010  to  0.015  mm  of  mercury.  Evacuation  generally 
required  about  30  minutes  of  pumping  to  reach  the  required  pressure. 
At  this  point,  the  penetrometer  was  filled  with  mercury. 

Filling  the  penetrometer  with  mercury  was  accomplished  by  rotating 

the  filling  device  so  that  the  levpl  -in  ,-k~ 

du  tne  level  m  the  mercury  reservoir  would  just 

cover  the  opening  in  the  penetrometer  stem.   The  vacuum  pump  was  turned 
off,  and  air  was  bled  into  the  assembly  to  a  pressure  of  20  mm  of 
mercury.   Raising  the  pressure  in  the  filling  device  forced  mercury  into 
the  penetrometer,  enveloping  the  soil  specimen.   In  some  cases  an  air 
bubble  would  remain  trapped  between  the  top  of  the  sample  and  the  glass 
surface  of  the  penetrometer.   When  this  occurs,  Bhasin  (1975) 
recommended  that  the  pressure  be  increased  to  40  ™  or  80  mm  of  mercury 
to  surround  the  soil  completely  with  mercury.   A  small  air  bubble  was 
always  present  after  filling  due  to  the  volume  of  air  in  the  bulb  prior 
to  filling.   If  the  bubble  appeared  to  be  trapped  above  the  sample  and 
if  it  was  rather  large,  the  filling  pressure  was  raised,  but  if  the 
bubble  was  small,  the  filling  pressure  was  left  at  20  mm  of  mercury. 

After  filling  the  penetrometer,  the  filling  device  was  rotated 
back  to  its  original  position  breaking  off  the  mercury  column  at  the 
Penetrometer  stem.   The  stem  reading  and  the  filling  pressure  were 
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recorded.   The  pressure  was  then  raised  in  increments,  noting  the 
position  of  the  stem  reading  after  each  increment.   This  process  is 
referred  to  as  the  low  pressure  intrusion  and  was  continued  to 
atmospheric  pressure. 

When  atmospheric  pressure  was  reached,  the  penetrometer  was 
removed  from  the  filling  device,  weighed  and  placed  in  the 
porosiraeter.   Because  the  penetrometer  fit  into  the  porosimeter  in  a 
"stem  down"  position,  the  specimen,  located  at  the  top  of  the  penetro- 
meter bulb,  was  under  tension  induced  by  the  suspended  column  of 
mercury  below  it.   To  bring  the  sample  at  atmospheric  pressure,  the 
porosimeter  pressure  was  raised  to  4  psi.   The  probe  reading  was 
recorded  at  4  psi,  and  the  pressure  was  increased  incrementally  to 
about  14,000  psi.   Probe  readings  were  taken  after  each  increment  of 
pressure  when  the  incremental  intrusion  was  complete. 
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4  -  RESULTS  AND  DISCUSSION  OF  RESULTS 

4-1  Compaction 

The  three  blends  of  clayey  silt  used  were  prepared  by  three  levels 
of  kneading  compaction  effort.   The  medium  effort  was  selected  such 
that  the  optimum  dry  density  of  the  90%  silt-10%  kaolin  mix  was  nearly 
equal  to  the  optimum  dry  density  of  the  same  soil  compacted  by  the 
AASHTO  T99-70  procedure  (i.e.  standard  Proctor).   The  high  compaction 
effort  was  intended  to  produce  an  optimum  dry  density  about  8  to  10% 
higher  than  the  medium  effort,  thus  approximating  the  modified  Proctor 
compaction  effort.   A  low  compaction  effort  was  also  used. 

After  generating  preliminary  compaction  curves  for  the  soils  and 
efforts  to  be  used,  the  samples  which  would  be  used  for  permeability 
and  pore  size  measurements  were  prepared  at  selected  water  contents. 
Three  molding  water  contents  were  generally  used  to  represent  a 
compaction  effort:  one  sample  was  prepared  near  the  optimum  water 
content,  one  was  prepared  about  1.5  to  3%  dry  of  optimum  and  one  was 
prepared  about  1.5  to  3%  wet  of  optimum  water  content. 

Table  3  summarizes  the  code  designation  used  to  identify  the  soils 
tested  in  this  study.   The  code  designation  gives  the  soil  type, 
compaction  effort,  and  relative  water  content  with  respect  to  its 
optimum.   Table  4  presents  the  compaction  parameters  (dry  density, 
water  content  and  degree  of  saturation)  measured  after  compaction  and 
again  after  permeation  for  the  twenty-two  samples  tested.   Permeability 
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TABLE  3  SAMPLE  CODE  DESIGNATION 


Soil  Type 

S9  -  90%  Silt  and  10%  Kaolin 

S7  -  70%  Silt  and  30%  Kaolin 

S5  -  50%  Silt  and  50%  Kaolin 

Compaction  Effort    (Hydraulic  System  Pressure) 

L     -  Low  Compaction  Effort ;  4.0  psi 

M    -  Medium  Compaction  Effort;  8.5  psi 

H    -  High  Compaction  Effort;  40.  psi 

Moisture  Condition 

0    -  Near  Optimum  Water  Content 

0-D  -  Slightly  Dry  of  Optimum  Water  Content 

D    -  Dry  of  Optimum  Water  Content 

W   -  Wet  of  Optimum  Water  Content 


Example  S9H0:  90%  Silt  and  10%  Kaolin  compacted  at 
40  psi  gauge  pressure  near  the  optimum 
water  content 
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TABLE   4     COMPACTION  PARAMETERS 

OF  SOILS  TESTED 

SAMPLE 

Dry 
Density 

Water 
Content 

Saturation 
After 

Water 
Content 

Saturation 
After 

CODE 

After 
Compaction 

Compactior 

i        After 
Permeation 

Permeation 

pcf 

% 

% 

% 

% 

S9M0 
S9M0-0 
S9MD 
S9MW 

109.7 
109.2 
106.7 
108 

16.2 
14.4 
11.7 
17.1 

81 
70 
54 
81 

18.6 
20.4 
21.4 
18.3 

93 
100 
99 
87 

S9HO 
S9HD 
S9HW 

120.4 
117.4 
115.4 

11.9 

7.5 

13.8 

79 
46 
80 

13.7 
17.3 
15.9 

91 

105 

92 

S7LO 
S7LD 
S7LW 

107.8 
103.3 
106.4 

19.0 
16.9 
20.3 

91 
72 
94 

19.9 
23.6 
20.8 

95 
101 
96 

S7M0 
S7MD 
S7MW 

I  13.9 
108.1 

II  1.4 

14.6 
13.5 
16.9 

82 
65 
89 

17.4 

21.4 
18.0 

93 
103 
94 

S7H0 
S7HD 
S7HW 

120.1 
118.0 
116.3 

12.7 
10.5 
14.8 

85 
66 
88 

14.9 
16.3 
16.3 

99 
103 
98 

S5LO 
S5LD 
S5LW 

105.6 
99.3 
103.9 

19.8 
18.3 
21.2 

90 

71.3 

92.6 

21.3 
26.1 
22.1 

97 
102 
97 

S5M0 
S5MD 
S5MW 

108.4 
99.5 
106.8 

18.1 
16.6 

19.9 

89 
65 
93 

19.9 
26.4 
20.3 

98 
103 
96 
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and  pore  size  distribution  measurements  were  made  for  each  of  the 
samples.   Figures  11,  12  and  13  present  the  compaction  curves,  the 
location  of  each  of  the  sample  points,  and  the  permeability  measured  for 
each  of  the  samples. 

4-2  Saturation 

As  mentioned  previously,  it  was  necessary  to  saturate  samples  prior 
to  permeation;  three  successive  steps  were  followed  for  the  saturation 
process.   Water  was  first  passed  through  the  sample  under  a  vacuum 
pressure  differential  of  about  10  psi.   Water  was  then  passed  through 
the  soil,  while  applying  back  pressure,  under  a  pressure  differential 
of  3  to  5  psi.   Finally,  the  permeability  tests  were  conducted  at  a  50 
psi  back  pressure.   The  vacuum  saturation  technique  was  very  effective 
in  increasing  the  degree  of  saturation  of  the  dry  side  samples.   Table 
4  shows  that  all  samples  compacted  dry  of  optimum  had  a  final  saturation 
of  about  100%.   However,  the  vacuum  saturation  may  have  caused  a  slight 
consolidation  in  some  of  the  samples  because  of  the  high  pressure 
differential  across  the  sample.   As  shown  in  Appendix  B,  seepage  forces 
cause  increases  in  effective  stress  which  can  equal  the  water  pressure 
differential  across  a  sample.   In  retrospect,  the  saturation  procedure 
could  have  been  improved  by  limiting  the  pressure  differential  during 
vacuum  saturation  to  5  psi  and  increasing  the  back  pressure  to  75  to 
100  psi. 

Table  4  shows  the  saturation  of  each  of  the  samples  after  per- 
meation as  measured  from  the  gross  sample.   It  is  important  to  note 
that  these  values  were  measured  after  the  50  psi  back  pressure  was 
removed;  thus,  the  saturation  during  permeation  was  not  known.   From 
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FIGURE    II     COMPACTION  CURVES  FOR    90%  SILT  AND 
10%  KAOLIN 
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FIGURE   12     COMPACTION  CURVES   FOR  70%  SILT  AND 
30%    KAOLIN 
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FIGURE  13    COMPACTION  CURVES  FOR  50%  SILT  AND 
50%  KAOLIN 
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the  theoretical  relationship  between  degree  of  saturation  and  back 
pressure  presented  by  Lowe  and  Johnson  (1960),  the  degree  of  saturation 
during  permeation  can  be  estimated.   At  50  psi  back  pressure,  an  initial 
degree  of  saturation  of  S .  =  93%  is  necessary  to  bring  a  sample  to  full 
(100%)  saturation,  and  S.  =  91%  is  required  to  bring  a  sample  to  99% 
saturation  (for  back  pressuring  at  constant  volume).   If  this  criterion 
is  used,  all  but  one  of  the  samples  were  99%  saturated  and  all  but  three 
of  the  samples  were  100%  saturated  during  permeation. 

4-3  Freeze  Drying 
It  was  critical  to  monitor  changes  in  volume  during  the  different 
phases  of  the  investigation  to  detect  sample  disturbance.   Table  5 
presents  porosity  measurements  made  after  compaction,  after  permeation, 
and  after  freeze  drying.   As  shown,  some  of  the  samples  experienced  a 
1  to  3%  change  in  volume  between  the  time  of  compaction  and  the  measure- 
ment of  pore  size  distribution.   Several  sources  of  error  are  likely  to 
have  caused  the  volume  changes.   Contrasting  columns  1  and  2  of  Table 
5,  it  is  apparent  that  some  samples  experienced  a  1  to  2%  increase  in 
volume  after  permeation.   This  probably  occurred  when  the  samples  were 
extruded;  a  slight  swelling  may  have  resulted  when  the  confining  mold 
was  removed.   Contrasting  columns  2  and  3  of  Table  5  indicates  that 
freeze  drying  may  have  also  resulted  in  slight  volume  changes  (both 
shrinkage  and  expansion)  of  1  to  3%  for  some  of  the  samples.   Reasons 
for  the  volume  changes  due  to  freeze  drying  were  discussed  previously 
in  Section  1-3.   Random  scatter  of  the  porosity  measurements  and 
sample  fabric  may  also  explain  some  of  the  volume  changes. 
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TABLE    5       POROSITY  MEASUREMENTS  BEFORE  AND  AFTER 
FREEZE   DRYING 


1 

2 

3 

SAMPLE 
CODE 

Porosity 
After 

Porosity  * 
After 

Porosity 
After 

Porosity 
Change 

Compaction 

Permeation 

Freeze  Drying 

(3-1) 

(Gross) 

(Avg) 

(Avg) 

S9M0 
S9M0-D 
S9MD 
S9MW 

.35 
.36 
.37 
.36 

.36 
.36 
.37 

.36 
.38 
.38 
.36 

0.01 
0.02 
0.01 
0.00 

S9K0 
S9HD 
S9HW 

.29 
.31 
.32 

.31 

.32 

.33 

.31 
.33 
.34 

0.02 
0.02 
0.02 

S7L0 
S7LD 
S7LW 

.36 
.39 
.37 

.36 
.39 
.37 

.36 
.39 
.37 

0.00 
0.00 
0.00 

S7M0 
S7MD 
S7MW 

.33 
.36 
.34 

.34 
.35 

.35 
.38 
.35 

0.02 
0.02 
0.01 

S7H0 
S7HD 
S7HW 

.29 

.30 
.31 

.29 
.31 

.29 
.33 
.30 

0.00 

0.03 

-0.01 

S5L0 
S5LD 
S5LW 

.37 
.41 
.38 

.38 
.38 

.37 
.40 
.37 

0.00 
-0.01 
-0.01 

S5MO 

S5MD 
S5MW 

.35 
.41 
.36 

.37 
.42 
.36 

.36 
.41 
.36 

0.01 

0.00 

0.00 

*"  Measured  by  wax  displacement.  Values  missing  were  too  soft  to  trim 
or  were  not  measured. 
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The  final  column  of  Table  5,  however,  indicates  that  the  maximum 
volume  change  from  compaction  to  pore  size  measurement  was  3%.   This 
was  felt  to  be  a  tolerable  amount . 

4-4  Permeability 

The  permeability  was  measured  by  a  closed  system,  falling  head 
apparatus  under  a  back  pressure  of  50  psi.   The  hydraulic  gradient 

during  the  tests  varied  from  10  to  1.   The  range  of  permeability  tested 

-4      -8 
extended  from  10   to  10   cm/sec.   Appendix  B  discusses  the 

permeability  test  and  calculation  procedures  used. 

Since  the  temperature  was  not  controlled  during  the  permeability 
tests,  it  was  monitored  with  each  water  level  reading.   The  permea- 
bility readings  were  normalized  to  a  20  C  temperature  as  suggested  by 
Lambe  and  Whitman  (1969)  as  follows: 


k-°c  ■  T^  k*  <4-» 


where  k  o   is  the  permeability  at  20  C 

k     is  the  permeability  at  temperature  T 
y„no   is  the  absolute  viscosity  of  water  at  20  C 
]i  is  the  absolute  viscosity  of  water  at  T  C 

T     is  the  average  temperature  between  two  water  level 
readings  in  C. 

Shown  on  Figures  14  through  17  are  the  permeability  measurements 

plotted  versus  mean  hydraulic  gradient  for  each  of  the  samples.   It  is 

apparent  that  some  of  the  samples  did  not  perform  in  accordance  with 

Darcy's  law,  which  requires  that  the  permeability  remain  constant  with 

changes  in  hydraulic  gradient.   Some  of  the  samples,  such  as  S7LW,  S7L0, 

S7MW  etc.,  initially  showed  a  decrease  in  permeability  with  decreasing 
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FIGURE  14     PERMEABILITY  VS   HYDRAULIC  GRADIENT  FOR 
S9  SAMPLES 
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FIGURE   15    PERMEABILITY  VS  HYDRAULIC  GRADIENT  FOR 
S7L  SAMPLES 
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FIGURE    16      PERMEABILITY  VS  HYDRAULIC  GRADIENT  FOR 
S7M   AND  S7H  SAMPLES 


85 


100 


10 
to 

o 

x 

o 

a* 

< 
E 

u 
c 

o 

o 

O 

C\J   |.0 

o 


o 


S5LD   74 


/ 


S5MD,  3.3 


S5L0,0.I6 


S5LW 
0.14 


S5M0.0.I2 

LEGEND 
•  Initial  Permeation 

S5M0      Sample  Code 

0.12      Mean  Permeability  in 
cm/sec  x  10"6 


S5MW,  0.063 


0.015 


4  6 

Mean  Hydraulic  Gradient 


10 


FIGURE    17     PERMEABILITY  VS  HYDRAULIC  GRADIENT  FOR 
S5  SAMPLES 
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hydraulic  gradient  then  tended  to  level  out  to  a  constant  permeability 
value.   This  type  of  behavior  was  initially  thought  to  be  caused  by  non- 
Darcy  type  flow  discussed  previously  (Miller  and  Low,  1963  or  Hansbo, 
1960);  however,  a  closer  examination  of  the  data  revealed  that  this  was 
not  so.   Sample  S7LW  will  be  used  as  an  example.   This  sample  was 
permeated  twice,  and  as  shown  in  Figure  15,  the  permeability-hydraulic 
gradient  relation  was  different  for  each  of  the  runs.   Figure  18  shows 
the  relation  between  permeability  and  net  cumulative  volume  of  inflow 
plotted  versus  time  for  sample  S7LW.   The  net  cumulative  volume  of 
inflow  was  determined  as  described  in  Appendix  B.   The  Figure  clearly 
indicates  that  initially,  steady-state  flow  did  not  exist  in  the 
sample,  thus  rendering  the  permeability  calculations  inaccurate.   A 
probable  cause  for  the  initial  net  inflow  of  water  into  the  sample  was 
that  the  sample  may  have  lost  part  of  its  50  psi  back  pressure  when  it 
was  transferred  from  the  saturation  device  to  the  permeameter.   When 
permeation  began,  several  hours  were  required  to  equilibrate  the  flow 
in  the  sample.   As  shown  in  Figure  18,  once  steady-state  flow  conditions 
were  approached,  the  permeability  measurements  tended  to  stabilize, 
agreeing  with  Darcy's  law. 

The  permeability  values  assigned  to  each  of  the  samples  were 
calculated  from  an  average  of  six  to  eight  measurements  from  the 
stabilized  portion  of  the  permeability-hydraulic  gradient  curves.   The 
scatter  in  the  permeability  values  was  partially  due  to  the  accuracy 
of  the  water  level  readings  and  temperature  fluctuations  between 
readings.   Only  one  sample,  S5LW,  shown  on  Figure  17,  exhibited  some 
peculiar  permeability  fluctuations  which  could  not  be  explained. 
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FIGURE   18      PERMEABILITY  AND  NET  VOLUME  INFLOW  VS 
TIME  FOR  SAMPLE  S7LW 
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The  relation  between  the  compaction  curves  and  the  permeability 
values  are  shown  on  Figures  11,  12  and  13.   The  results  generally  agree 
with  trends  found  by  Lambe  (1954),  Bjerrum  arid  Huder  (1957),  and 
Mitchell  et  al.  (1965),  which  were  previously  discussed.   The  dry  side 
samples  generally  had  permeabilities  one  to  three  orders  of  magnitude 
greater  than  the  optimum  and  wet  side  samples.   There  was  little 
difference  between  the  permeabilities  of  the  optimum  and  wet  side 
samples  for  a  given  compaction  effort  and  soil  type. 

4-5  Pore  Size  Distribution 


Pore  size  distribution  curves  are  shown  on  Figures  19  through  32 
and  in  Appendix  A  on  Figures  37  through  58.   Each  of  the  figures 
consist  of  a  differential  and  cumulative  graph  of  the  porosity  (volume 
of  voids  intruded  divided  by  the  total  volume)  versus  the  logarithm  of 
the  pore  diameter.   Also  shown  on  the  Figures  are  the  compacted  water 
contents  and  the  permeabilities  of  each  of  the  samples.   The  differential 
plots  were  determined  as  described  in  Section  1-2.4.   The  large  data 
points  along  the  ordinate  axis  of  the  cumulative  plots  represent  the 
porosity  of  the  pore  size  specimen  as  measured  by  mercury  displacement. 
Several  general  observations  are  evident  from  the  pore  size  distribution 
figures . 

The  differential  plots  provide  a  useful  complement  to  the 
cumulative  plots  for  representing  pore  size  distributions.   The  most 
volumetrically  frequent  or  "modal"  pore  sizes  are  clearly  depicted  as 
peaks  on  the  differential  plots.   When  contrasting  separate  pore  size 
distribution  curves,  the  location  and  size  (both  height  and  width)  of 
these  peaks  serve  as  useful  comparative  descriptive  parameters. 
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Hereafter,  the  most  frequent  pore  diameters  will  be  referred  to  by  their 
proper  statistical  name,  modes. 

The  pore  size  distribution  measurements  of  the  samples  tested  in 
this  study  displayed  two  modes  (bimodal)  when  considered  on  a  log 
diameter  scale.   The  relative  position  and  size  of  these  modes  gave 
clear  indications  of  the  nature  of  the  soil  fabric.   One  of  the  modes 
was  between  10  and  1.0  ym  and  will  be  referred  to  as  the  large  pore  mode. 
The  other  mode  occurred  consistently  at  0.1  \xn   and  will  be  referred  to 
as  the  small  pore  mode. 

All  of  the  pore  size  distribution  curves  showed  little  or  no  volume 
of  pores  larger  than  10  \sm. 

4-5.1  Replication 

The  reliability  and  reproducibility  of  pore  size  distribution 
measurements  was  one  of  the  topics  addressed  in  this  study.   As 
mentioned  previously,  several  precautions  were  taken  to  insure 
homogeneous  compacted  samples.   A  minimum  of  two  pore  size  distributions 
were  measured  from  each  sample.   The  plots  of  these  pore  size  measure- 
ments are  presented  in  Appendix  A,  Figures  37  through  58.   Overall 
the  replication  of  the  pore  size  distribution  measurements  was  felt  to 
be  quite  good.   Only  one  sample,  S7M0  (shown  on  Figure  47),  produced 
significant  scatter  between  pore  size  measurements. 

The  most  obvious  reason  for  scatter  in  replicate  pore  size 
distribution  curves  is  natural  variability  of  the  soil  fabric.   Some  of 
the  specimens  may  have  suffered  more  compression  or  disturbance  than 
others.   Such  differences  in  fabric  may  be  detected  by  comparing  the 
location  and  size  of  the  large  pore  modes  and  comparing  the  porosities 
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of  each  of  the  specimens.  Examples  of  these  discrepancies  are  shown 
by  S9H0  (Figure  41),  S9HW  (Figure  43),  S7H0  (Figure  50)  S7HW  (Figure 
52),  and  S5LD  (Figure  54). 

Measurement  errors  may  also  account  for  scatter  in  replicate  pore 
size  distribution  curves.   Significant  differences  in  the  cumulative 
pore  size  curves  may  result  because  of  a  high  initial  intrusion  in  some 
of  the  specimens.   This  intrusion  at  the  larger  diameters  is  probably 
caused  by  cracks  in  the  specimen  due  to  poor  drying  techniques  or  by  a 
trapped  air  pocket  on  the  surface  of  the  specimen  when  it  was  initially 
surrounded  with  mercury.   Such  specimens  will  give  erroneous  porosity 
measurement  which  are  too  high.   The  cumulative  plots  of  these  samples 
will  not  coincide,  but  the  differential  plots  should  be  nearly 
identical.   Examples  are  S9HD  (Figure  42)  and  S5LW  (Figure  55). 

Measurement  errors  may  occur  when  the  pore  volume  of  a  specimen 
exceeds  the  intrusion  capacity  of  the  porosimcter  (the  probe  limit). 
Such  specimens  yield  a  negative  intrusion  at  the  higher  pressures 
(smaller  pore  diameters)  as  shown  by  samples  S7LD  (Figure  45)  and 
S5LD  (Figure  54).   The  apparent  negative  intrusions  result  from  the 
corrections  for  the  compression  of  mercury. 

If  differential  pore  size  curves  are  plotted,  it  is  imperative 
that  the  pressure  increments  be  taken  at  logarithmically  equal 
intervals;  otherwise,  serious  distortions  of  the  differential  curves 
will  result,  as  disscussed  in  Section  1-2.4.   This  fact  was  not  fully 
appreciated  in  this  study  until  after  the  pore  size  distribution 
measurements  were  made.   As  mentioned  previously,  a  proper  pressure 
increment  schedule  for  mercury  intrusion  should  be: 
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where  P  is  the  intrusion  pressure  and  "c"  is  the  logarithmic  interval 
constant.   The  value  of  c  generally  varied  between  0.16  and  0.23 
with  an  average  of  0.2  for  pores  smaller  than  50  ym  in  size.   In 
retrospect,   c  should  have  been  kept  constant  during  all  the  pore  size 
measurements.   Samples  S7L0  (Figure  44)  and  S7LW  (Figure  46)  exhibit 
examples  of  the  distortions  which  can  occur.   Both  figures  show  slight 
discrepancies  of  the  small  pore  mode  because  of  poor  pressure  spacings 
in  the  small  pore  diameter  range. 

Several  of  the  samples,  including  S9M0  (Figure  37),  S9M0-D 
(Figure  38),  S9MD  (Figure  39)  and  S5M0  (Figure  56),  show  excellent 
replication  of  pore  size  distribution  measurements. 

The  specimen  with  the  porosity  closest  to  the  gross  compacted 
sample  porosity  was  chosen  to  represent  a  given  sample,  if  no  major 
measurement  errors  were  evident  in  its  pore  size  distribution  curves. 
The  first  specimen  which  appears  in  the  legend  of  the  pore  size  figures 
in  Appendix  A  was  the  one  selected  to  represent  each  of  the  samples. 

4-5.2  Contrasting  Pore  Size  Distributions 
This  section  discusses: 

1)  changes  in  pore  size  distribution  produced  by  varying 
compactive  water  content  and  compactive  effort  for  a 
given  soil,  and, 

2)  pore  size  distributions  of  different  soils  prepared  at 
similar  compactive  effort  and  moisture  conditions 
relative  to  their  optimum  water  contents. 
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The  90%  silt  -  10%  kaolin  is  the  first  soil  examined.   Its  compac- 
tion curves  are  presented  on  Figure  11,  and  its  pore  size  distribution 
curves  are  shown  on  Figures  19,  20  and  21.   From  the  pore  size  figures 
it  is  evident  that  the  pore  size  distributions  are  bimodal  with 
differences  in  the  large  pore  mode  accounting  for  the  changes  in  soil 
fabric,  while  the  small  pore  modes  remain  nearly  constant.   Figures 
19  and  20  show  the  large  pore  mode  decreasing  with  increasing  water 
content  on  the  dry  side  of  optimum.   Only  minor  differences  in  the 
large  pore  mode  are  present  between  optimum  and  wet  side  samples  for 
a  given  compactive  effort.   The  permeability  values  of  each  of  the 
samples  increase  with  increasing  large  pore  mode.   Figure  21  shows 
that  as  the  compactive  effort  is  increased,  the  large  pore  mode 
decreases  for  both  dry  side  and  optimum  samples.   The  permeability, 
again,  decreases  with  decreasing  large  pore  mode.   This  same  Figure 
also  indicates  that  the  small  pore  mode  is  not  affected  by  changes  in 
compactive  effort  or  moisture  condition;  its  frequency  and  location 
remain  fixed. 

Figure  12  shows  the  compaction  curves  of  the  70%  silt-30%  kaolin 
soil  and  Figures  22  through  26  show  its  pore  size  distribution  curves. 
For  each  of  the  three  compaction  efforts  the  dry  side  samples  have  a 
significantly  greater  large  pore  mode  than  the  optimum  and  wet  side 
samples,  as  shown  by  Figures  22,  23  and  24.   The  wet  side  and  optimum 
samples  have  nearly  identical  differential  distributions  for  the  high 
and  medium  compactive  efforts;  however,  for  the  low  compactive  effort, 
the  wet  side  sample  has  a  slightly  higher  large  pore  mode  than  its 
optimum  counterpart.   The  permeabilities  of  the  optimum  and  wet  side 
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FIGURE   19     DIFFERENTIAL  AND  CUMULATIVE    PORE  SIZE  DISTRIBUTION 
CURVES    FOR  90%  SILT -10%  KAOLIN  COMPACTED  AT 
MEDIUM  EFFORT 
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FIGURE  20     DIFFERENTIAL  AND   CUMULATIVE  PORE  SIZE  DISTRIBUTION 
CURVES    FOR  90%  SILT- 10%  KAOLIN  COMPACTED  AT 
HIGH  EFFORT 
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FIGURE  21 


PORE  SIZE  DISTRIBUTION  CURVES  FOR  90%  SILT-  10% 
KAOLIN  COMPACTED  AT  OPTIMUM   AND  DRY  OF  OPTIMUM 
WATER  CONTENTS  FOR  TWO  EFFORTS 
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FIGURE  22     DIFFERENTIAL  AND  CUMULATIVE  PORE  SIZE  DISTRIBUTION 
LOW^FORT  7°%  S,LT"30°/o  KA0LIN  COMPACTED  AT 
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FIGURE  24     DIFFERENTIAL    AND  CUMULATIVE  PORE  SIZE  DISTRIBUTION 
CURVES   FOR  70%  SILT-  30%  KAOLIN   COMPACTED  AT 
HIGH  EFFORT 
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FIGURE  25    PORE   SIZE  DISTRIBUTION  CURVES  FOR   70%  SILT-  30% 
KAOLIN  COMPACTED  AT  OPTIMUM  WATER  CONTENT  FOR 
THREE  EFFORTS 
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FIGURE  26    PORE  SIZE  DISTRIBUTION  CURVES  FOR  70%  SILT-  30% 

KAOLIN  COMPACTED  DRY  OF  OPTIMUM  FOR  THREE  EFFORTS 
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samples  are  nearly  equal  at  each  of  the  three  compactive  efforts.   The 
effect  of  increasing  compactive  effort  on  optimum  and  dry  side  samples 
is  evident  on  Figures  25  and  26,  respectively.   The  magnitude  of  the 
large  pore  mode  remains  unchanged  for  the  optimum  sample;  however,  the 
frequency  (or  height)  of  the  mode  decreases  with  increasing  compactive 
effort.   The  same  is  true  of  the  dry  side  samples;  increasing  the 
compactive  effort  decreases  the  frequency  of  the  large  pore  mode  but 
does  not  change  its  value.   A  small  discrepancy  is  shown  in  Figure  25 
between  the  permeability  and  pore  size  distribution  of  samples  S7L0  and 
S7M0.   From  the  differential  distribution  curves,  sample  S7L0  would  be 
expected  to  have  a  slightly  higher  permeability  than  S7M0,  but  the 
opposite  is  true.   The  small  pore  mode  of  the  70%  silt  -  30%  kaolin 
soil,  unlike  the  90%  silt  -  10%  kaolin  soil,  does  seem  to  undergo  a 
very  slight  decrease  with  increased  compactive  effort,  as  shown  in 
Figures  25  and  26. 

The  50%  silt  -  50%  kaolin  soil  exhibits  a  distinctly  different 
fabric  pattern  than  the  other  two  soil  types  tested.   Figure  13  shows 
the  compaction  curves,  and  Figures  27,  28  and  29  display  the  pore  size 
curves.   The  pore  size  distributions  indicate  that  samples  compacted 
at  optimum  or  wet  of  optimum  water  content  have  little  or  no  large  pore 
mode,  while  dry  side  samples  exhibit  quite  a  distinct  large  pore  mode. 
This  explains  why  the  50%  -  50%  soil  was  the  most  sensitive  to  changes 
in  permeability  between  dry  side  and  wet  side  samples.   Figure  29  is 
inconclusive  about  the  affect  of  increased  compaction  effort  on  the 
small  and  large  pore  mode.   It  would  appear  that  the  increased  effort 
did  little  to  alter  the  dry  side  and  optimum  sample  pore  size  distri- 
butions, respectively. 
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FiGURE  27     DIFFERENTIAL  AND  CUMULATIVE  PORE  SIZE  DISTRIBUTION 
CURVES   FOR  50%  SILT -50%  KAOLIN  COMPACTED   AT 
LOW  EFFORT 
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FIGURE  28     DIFFERENTIAL  AND  CUMULATIVE    PORE  SIZE  DISTRIBUTION 
CURVES  FOR  50%  SILT- 50%  KAOLIN  COMPACTED  AT 
MEDIUM  EFFORT 
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A  contradiction  between  the  permeabilities  and  pore  size  distri- 
butions of  the  two  dry  side  samples,  S5LD  and  S5MD,  is  shown  on 
Figure  29.   The  two  samples  appear  to  have  nearly  identical  pore  size 
distributions;  however,  the  permeability  of  S5LD  is  an  order  of 
magnitude  greater  than  that  of  S5MD.   The  cause  of  the  discrepancy  is 
unknown.   However,  sample  S5LD  was  the  only  sample  which  was  too  fragile 
to  trim  for  wax  displacement  measurements;  it  literally  crumbled  apart 
when  it  was  extruded  from  the  permeameter  mold.   The  low  compactive 
effort  and  low  initial  degree  of  saturation  resulted  in  a  sample  with 
an  extremely  small  soaked  strength.   This  made  the  task  of  trimming 
pore  size  specimens  rather  difficult.   Because  of  these  factors,  it  is 
believed  that  the  true  pore  size  distribution  of  sample  S5LD  was  under- 
estimated. 

Figure  30  is  a  particle  size  frequency  plot  of  the  three  soil 
types.   The  frequency  was  determined  by  the  same  methods  used  to 
calculate  the  differential  pore  size  curves.   Logarithmically  equal 
particle  size  intervals  were  taken  from  the  cumulative  grain  size 
distributions  (Figure  3).   The  specific  gravity  of  solids  was  assumed 
to  be  constant,  which  introduced  a  slight  error  in  the  relative 
frequencies  of  the  clay  and  silt  fractions.   However,  the  error  was 
not  felt  to  be  significant. 

Figures  31  and  32  contrast  the  pore  size  distributions  of  the 
three  soil  types  all  compacted  at  medium  effort  and  optimum  and  dry  of 
optimum  water  contents,  respectively.   Comparing  the  particle  size 
and  pore  size  volumetric  frequency  plots  it  is  evident  that  both  are 
bimodal;  although,  the  complete  distribution  of  the  clay  size 
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FIGURE  31       PORE  SIZE    DISTRIBUTION  CURVES  FOR  THREE  SOIL  TYPES 
COMPACTED  AT  MEDIUM  EFFORT  AND  OPTIMUM  WATER 
CONTENT 
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fraction  was  not  obtained.   The  large  pore  mode  of  the  soils  increases 
with  increasing  silt  content,  while  the  small  pore  mode  remains  at 
0.1  ym,  but  increases  in  frequency  with  increasing  clay  content.   The 
larger  pores  are  probably  inter-silt  and  inter-aggregate  pores,  and  the 
smaller  pores  are  intra-aggregate  pores. 

Comparing  the  permeability  values  and  pore  size  distributions  of 
the  three  soils  shown  on  Figures  31  and  32,  the  permeability  decreases 
with  a  decrease  in  the  large  pore  mode. 

For  the  soils  tested  the  position  and  frequency  of  the  large  pore 
mode  appears  to  be  the  best  indicator  of  changes  in  soil  fabric  caused 
by  varying  water  content  and  compactive  effort  for  a  given  soil  type. 
The  distribution  of  the  smaller  pores  remained  relatively  fixed  for  each 
of  the  soil  types. 

Of  the  three  soil  types,  only  the  50%  silt  -  50%  kaolin  behaved  in 
accordance  with  the  deformable  aggregate  model  of  clay  fabric  proposed 
by  Barden  and  Sides  (1970)  and  Hodek  (1972).   For  water  contents  at 
or  greater  than  optimum  very  few  large  pores  remained  and  no  significant 
large  pore  mode  was  evident. 

The  other  two  soil  types  were  strongly  influenced  by  the  large 
fraction  of  silt  present.   A  significant  volume  of  large  pores  always 
remained,  regardless  of  the  compactive  effort  or  water  contents  used. 

Based  on  the  rather  limited  number  of  soils  used  in  this  stud);  it 
is  hypothesized  that  the  shape  of  the  differential  pore  size  distri- 
bution in  the  small  pore  region  (<  0.02  Urn)  may  serve  as  an  identifi- 
cation or  "signature"  of  a  given  soil  type.   The  value  of  the  small 
pore  mode  is  no  doubt  related  to  the  size  of  clay  particles  present, 
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while  the  frequency  is  related  to  the  amount  of  clay  in  the  soil. 
Factors  such  as  plasticity,  swelling  potential  and  clay  mineral  content 
are  probably  related  to  the  characteristics  of  the  small  pore  mode. 
Further  research  would  be  required  to  test  this  hypothesis. 

4-6  Relation  Between  Permeability  and  Pore  Size  Distribution 

From  the  qualitative  observations  discussed  in  the  previous  section 
it  is  evident  that  the  permeability  is  related  to  the  distribution  of 
pores  about  the  large  pore  mode.   The  objective  of  this  section  is  to 
quantify  this  relationship. 

Previous  investigators  including  Taylor  (1948)  and  Lambe  and 
Whitman  (1969)  have  shown  that  the  log  of  permeability  is  linearly 
related  to  the  void  ratio  of  some  soils.   From  Figure  33  it  is  apparent 
that  the  relationship  is  not  valid  for  compacted  fine-grained  soils. 
This  agrees  with  the  fact  that  the  void  ratio  had  little  or  no 
relationship  to  the  large  pore  mode  and  its  frequency.   This  Figure 
clearly  shows  that  gross  volumetric  parameters  cannot  account  for 
changes  in  permeability  of  compacted  soils;  consequently,  any 
permeability  predicting  equations  based  on  these  volumetric  parameters 
is  not  suitable  for  compacted  soils. 

Section  1-1.3  examined  three  models  which  relate  permeability  to 
pore  size  distribution  parameters.   These  are  summarized  in  Table  1. 
The  specific  permeability  K  can  be  expressed  in  terms  of  each  of  these 
three  parameters  as: 


K  =  C   •  PSP  C4-2) 
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2    A9 
where  PSP  is  one  of  the  three  pore  size  parameters,  E(d  )n,  d  or 

R^  n,  and 

C   is  a  shape  factor  which  differs  for  each  of  the  models. 
s 

Equation  (4-2)  can  be  written  in  terms  of  the  permeability  k,  since 

v 
k  ■  K  -k  with  only  a  change  in  the  value  of  the  shape  factor,  as 

follows: 


k  =  C  -PSP  (4-3) 

s 


The  PSP  values  were  calculated  for  each  of  the  samples  from  the 
pore  size  distribution  data  according  to  the  formulae  shown  on  Table  1 
and  the  method  discussed  in  Section  1-2.4.   These  values  are  shown  on 
Table  6  along  with  the  permeability  of  each  of  the  soils.   An  important 
assumption  had  to  be  made  to  calculate  the  Capillary  model  and  Marshall 
model  parameters.   As  shown  on  Table  1,  the  pore  size  parameters  for 
these  two  models  are  very  sensitive  to  the  larger  pore  diameters.   In 
order  to  obtain  a  relationship  between  permeability  and  the  parameters 
it  was  necessary  to  consider  only  those  pores  less  than  10  urn  in  size. 
The  selection  of  10  ym  was  arbitrary.   However,  as  shown  by  the 
differential  pore  size  distributions,  none  of  the  soils  experienced 
significant  or  continuous  intrusion  for  pores  greater  than  10  ym. 

The  calculation  of  the  Hydraulic  Radius  pore  size  parameter  also 
required  an  assumption.   As  shown  on  Table  1,  the  Hydraulic  Radius 
parameter  is  most  influenced  by  the  smaller  pores.   The  mercury 
intrusion  equipment  used,  left  about  1  to  5%  of  the  pore  volume 
unintruded.   This  unintruded  pore  volume  was  assigned  a  diameter  of 
0.01  Um  for  calculation  purposes. 
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TABLE    6        PERMEABILITY  AND  EQUIVALENT  PORE  SIZE 
PARAMETERS 


Sample 
Code 


Permeability 

k  at  20°  C 

(xlO"6cm/sec) 


Hydraulic 
Radius  Modes 

_     2 

Rh    n 
(xlO'V*2) 


Capillary 
Model  E(d2)-n 

for  d<  10/xm 
(^m2) 


Marshall 
Moael 
d<IO^.m 
d~2   (/tm2) 


S9M0 

S9M0-D 

S9MD 

S9MW 

S9H0 

S9HD 

S9HW 

S7L0 

S7LD 

S7LW 

S7M0 

S7MD 

S7MW 

S7H0 

S7HD 

S7HW 

S5L0 

S5LD 

S5LW 

S5M0 

S5MD 

S5MW 


4.3 
23 
39 
1.7 

0.81 

II 

0.50 

0.25 
16 

0.25 

0.40 

7.3 

0.12 

0.035 

0.89 

0.049 

0.16 
74 

0.14 

0.12 

3.3 

0.063 


4.3 

7.8 

8.1 

4.9 

3.8 

4.7 

2.4 
1.78 

3.4 

2.7 
1.03 

2.1 

0.83 

0.81 
1.23 
1.05 
1.16 
1.89 
1.42 
1.04 
1.74 

0.98 


2.8 
5.5 
7.5 
2.4 
1.39 
3.8 
2.0 
0.86 
2.1 

1.13 
0.85 

1.58 
0.83 
0.39 

1.00 
0.61 
0.48 

1.30 
0.61 
0.62 

1.23 
0.56 


0.51 
1.02 
1.39 
0.39 
0.15 
0.73 
0.29 
0.090 
0.27 
0.14 
0.084 
0.19 
0.078 
0.020 
0.1 1 
0.032 
0.041 
0.1  I 
0.063 
0.043 
0.12 
0.048 
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The  values  of  the  pore  size  parameters  of  the  Capillary,  Marshall 
and  Hydraulic  Radius  models  are  plotted  versus  the  permeability  on 
Figures  34,  35  and  36,  respectively.   All  axes  are  on  a  logarithmic 
scale.   A  linear  trend  in  the  data  is  apparent  from  these  figures.   The 
sample  with  the  largest  permeability,  which  is  farthest  removed  from  the 
linear  trend  on  all  these  Figures,  is  sample  S5LD.   As  discussed 
previously,  the  pore  size  distribution  measurement  of  this  sample  is 
believed  to  be  in  error. 

Simple  linear  regression  by  the  method  of  least  squares  was  used 
to  estimate  the  relationship  between  the  permeability  and  the  pore  size 
parameters.   Details  of  this  method  are  presented  in  Wapole  and  Myers 
(1972).   The  log  of  the  permeability  was  taken  as  the  dependent 
variable  and  the  log  of  the  pore  size  parameter  was  taken  as  the 
independent  variable.   The  linear  relation  is  expressed  as: 

log  k  =  b  log  PSP  +  a  (4-4) 

where  b  and  a  are  the  linear  regression  parameters.   Equation  (4-4)  may 
be  simplified  by  raising  it  to  the  power  of  10,  resulting  in: 


k  =  10a  PSPb  (4-5) 


Note  that  equation  (4-3),  determined  from  the  permeability-pore  size 
models,  is  a  form  of  equation  (4-5)  with  b  =  1  and  C  =  10  .   Thus, 


equation  (4-5)  may  be  expressed  as; 


k  =  C  PSPb  (4-6) 

s 


where  C  =  10  . 

s 
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The  C   and  b  parameters  and  the  coefficient  of  determination 

2 
"r  ''  were  calculated  for  each  and  all  of  the  soil  types  tested.   Table 

2 
7  presents  the  linear  regression  equations  and  r  values  for  each  of 

the  three  pore  size  parameters.  i 

The  accuracy  of  the  regression  equations  can  be  judged  by  the  b 

2 

and  r  parameters.   Comparing  equations  (4-3)  and  (4-6),  it  is  evident 

that  the  b  parameter  reflects  the  sensitivity  of  the  permeability  to 
the  pore  size  parameter  PSP.   The  closer  b  is  to  1,  the  more  sensitive 

the  parameter  is  to  the  permeability.   The  coefficient  of  determination, 

2 
r  ,  is  a  measure  of  the  degree  of  association  between  the  log  of  the 

permeability  and  the  log  of  the  pore  size  parameter  for  a  given  equation, 

2 
The  closer  r  is  to  1.0 'the  higher  the  degree  of  association  between  the 

2 
two  terms.   Consequently,  the  closer  b  is  to  1,  and  the  higher  the  r 

value,  the  closer  the  prediction  equation  fits  a  given  theoretical 

model. 

Table  7  indicates  that  for  the  Capillary  and  Marshall  models,  the 

2 
b  parameter  increases  and  the  r  value  decreases  with  increasing  clay 

2       ~2 
content.   The  two  pore  size  parameters,  E(d  )n  and  d  ,  are  most  affected 

by  the  larger  pores.   It  is  believed  that  this  decrease  in  sensitivity 

of  the  pore  size  parameters  with  increasing  clay  content  is  a  result  of 

the  d  <  10  ym  cut-off  diameter  used  to  calculate  the  parameters.   The 

soils  with  the  higher  clay  content  did  not  have  a  significant  intrusion 

for  pore  diameters  greater  than  about  6  ym;  therefore,  the  pore  size 

parameters  were  not  as  sensitive  to  changes  in  permeability.   The  90% 

silt-10%  kaolin  samples  display  an  excellent  linear  trend  between  the 

Capillary  and  Marshall  parameters  and  the  permeability  as  shown  in 

Figures  34  and  35  and  by  the  regression  equations  on  Table  7.   None  of 
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TABLE  7  LINEAR  REGRESSION  EQUATIONS  FOR  PERMEABILITY  AND  PORE  SIZE 
PARAMETERS  (PSP) 


Capillary  Model 

90%  silt  -  10%  kaolin 
70%  silt  -  30%  kaolin 
50%  silt  -  50%  kaolin 
all  samples 
all  except  S5LD 


PSP  =  E(d2)n  in  pro2 


k  in  cm/ sec  x  10 


-6 


k  =  0.194  (PSP) 

k  =   0.567  (PSP) 

k  =   3.487  (PSP)' 

k  =   0.532  (PSP) 

k  =   0.434  (PSP) 


2.73 
3.945 


2.42 
2.42 


0.91 
0.87 
0.83 
0.69 
0.83 


Marshall  Model 

90%  silt  -  10%  kaolin 

70%  silt  -  30%  kaolin 

50%  silt  -  50%  kaolin 

all  samples 

all  except  S5LD 

Hydraulic  Radius  Model 

90%  silt  -  10%  kaolin 
70%  silt  -  30%  kaolin 
50%  silt  -  50%  kaolin 
all  samples 
all  except  S5LD 


~2       2 
PSP  =  d   in  wm 

k  in  cm/scc  x  10 


-6 


k  =   17.41    (PSP)2,06 

0.88 

k  =   112.4    (PSP)2-28 
k  =   4.6  x  105    (PSP)A 

.95 

0.80 
0.73 

k  =   23.60    (PSP)1,62 

0.64 

k  =   20.57    (PSP)1'67 

0.82 

2                     2 
PSP  -   RH  n      in   pin     x 

10"A 

k  in   cm/sec  x   10 

k  =   0.015    (PSP)3'60 

0.82 

k  =   0.141    (PSP)2'95 

0.55 

k  =   0.046    (PSP)8'91 

0.79 

k  =   0.166    (PSP)2'43 
k  =  0.129    (PSP)2,A7 

0.57 
0.71 
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the  pore  size  parameters  were  as  sensitive  to  permeability  as  predicted 
by  the  theoretical  equations,  where  b  =  1. 

The  regression  equations  calculated  from  all  of  the  samples  except 

S5LD  for  the  Capillary  and  Marshall  models  are  shown  on  Figures  3A  and 

2 
35,  respectively.   The  b  parameter  and  r  value  indicate  that  the 

equation  is  fairly  accurate;  however,  several  systematic  errors  are 

shown  by  the  Figure.   At  the  higher  permeabilities  the  90%  silt  -  10% 

kaolin  soils  generally  fall  along  or  above  the  regression  line,  and 

the  soils  with  higher  clay  content  generally  fall  below  the  regression 

line. 

The  Hydraulic  Radius  model  is  not  as  useful  as  the  previous  two 

models  discussed.   Each  of  the  soil  types  result  in  a  regression 

2 
equation  with  either  a  very  high  b  parameter  or  a  low  r  value. 

Figure  36  shows  the  regression  line  calculated  from  all  the  samples 

2 
except  S5LD.   The  b  parameter  and  r  value  of  this  equation  (from 

Table  7)  are  completely  deceiving.   None  of  the  soil  types  fit  the 

regression  line  individually,  but  when  considered  en  masse  the  relation 

appears  adequate.   The  Hydraulic  Radius  model  is  not  believed  to 

produce  an  accurate  pore  size  parameter.   As  mentioned  previously,  the 

parameter  is  mostly  influenced  by  the  small  pore  mode,  which  is  not 

directly  related  to  the  large  pore  mode. 

As  discussed  previously,  the  permeability  is  strongly  influenced 

by  the  distribution  of  pores  about  the  large  pore  mode,  which  occurred 

between  10  and  1  yra  for  the  three  soils  tested.   It  is  not  a  coincidence 

that  the  Marshall  and  Capillary  pore  size  parameters  correlated  well 

with  the  permeability  only  when  pore  diameters  less  than  10  ym  were 

considered.   This  indicates  that  the  permeability  is  controlled  by  the 
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larger  pores  that  occur  in  sufficient  volume  or  frequency.   The 
frequency  is  probably  not  given  enough  weight  in  the  models  considered. 
The  Marshall  model  is  a  step  in  the  right  direction,  however.   As  shown 
in  Section  1-1.3,  the  Marshall  model  is  based  on  random  junctions  of 
pores  between  two  cross-sections  of  a  porous  medium;  consequently, 
the  frequency  term  appears  as  a  squared  term  in  the  pore  size  parameter. 
By  expanding  this  model  and  considering  "m"  cross-sections  in  sequence, 
the  frequency  term  would  be  given  greater  weight  and  the  permeability 
model  may  be  improved.   Further  research  would  be  required  to  confirm 
this. 

In  summary,  a  relationship  between  the  pore  size  distribution  and 
the  permeability  is  clearly  shown  on  Figures  34,  35  and  36.   The 
Marshall  or  Capillary  model  pore  size  parameters  may  provide  an 
excellent  empirical  predictive  equation  for  a  given  soil  type,  if  a 
proper  cut-off  diameter  is  selected  for  the  calculations  of  the 
parameters.   More  sophisticated  theoretical  permeability  models, 
accounting  for  the  complexities  of  the  pore  network,  are  necessary  to 
improve  the  prediction  of  permeability  from  pore  size  distribution. 
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5  -  CONCLUSIONS  AND  RECOMMENDATIONS 


5-1  Conclusions 


1.  A  method  is  described  for  approximating  and  presenting  the 
differential  pore  size  distribution  on  a  log  diameter  scale  which 

is  relatively  free  from  distortions.   The  differential  plot  provides 
a  useful  complement  to  the  cumulative  pore  size  distribution  curve 
by  showing: 

a.  the  type  of  distribution, 

b.  the  modal  or  most  frequent  pore  diameters, 

c.  any  gaps  or  irregularities  in  the  data, 

d.  simply,  how  different  distributions  compare. 

2.  Freeze  drying  was  an  effective  method  of  drying  the  soils  used. 
The  maximum  volume  change  caused  by  drying  was  3%. 

3.  Critical  region  drying  produced  a  chemical  alteration  of  the  soils 
tested.   Dolomite,  calcite  and  kaolinite  in  the  soil  reacted  to 
form  a  feldspar,  resulting  in  void  ratio  increases  of  35  to  40%. 

4.  The  pore  size  distributions  of  the  soils  tested  were  found  to  be 
bimodal  on  a  log  diameter  scale.   A  large  pore  mode  occurred 
between  10  and  1.0  urn,  and  a  small  pore  mode  occurred  at  0.1  ym. 

5.  Changes  in  soil  fabric  caused  by  varying  water  content  and  com- 
paction effort  were  generally  at  the  expense  of  the  large  pore 
mode. 
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6.  The  shape  of  the  differential  pore  size  distribution  about  the 
small  pore  mode  was  not  affected  by  the  compaction  variables  for 
a  given  soil  type. 

7.  The  relation  between  the  compaction  variables  and  the  permeability 
found  in  this  study  agree  with  the  trends  found  by  Lambe  (1954), 
Bjerrum  and  Huder  (1957),  and  Mitchell  et  al.  (1965). 

8.  The  permeability  generally  varied  with  the  magnitude  and  frequency 
of  the  large  pore  mode.   Decreasing  the  large  pore  mode  or  its 
frequency  resulted  in  a  decrease  of  the  permeability. 

9.  Three  theoretical  models  which  relate  pore  size  distribution  to 
permeability  were  examined.   Although,  the  models  were  not  completely 
successful  in  predicting  permeability,  the  pore  size  parameters  from 
two  of  the  models,  with  several  modifications,  did  reflect  the 
dependence  of  permeability  on  pore  size  distribution.   These  pore 
size  parameters  were  successfully  used  to  determine  empirical 
permeability  predicting  equations  for  any  one  of  the  soils  tested. 
The  general  form  of  the  predictive  equations  is: 


k  =  C  PSPb 
s 


where  k  is  the  permeability, 

PSP  is  the  pore  size  parameter  calculated  from  the  pore  size 
distribution,  and 

C  and  b  are  regression  constants. 

s  ° 

10.   Conventional  permeability  predictive  equations  based  on  volumetric 
and  grain  size  parameters  do  not  accurately  reflect  the  nature  of 
the  pore  system  of  compacted  fine-grained  soils,  and  therefore,  are 
unsuitable  for  estimating  permeability.   The  empirical  predictive 
equations  determined  in  this  study,  (1)  demonstrate  the  dependence 
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of  permeability  on  pore  size  distribution  and,  (2)  are  reflective 
of  changes  in  soil  structure  caused  by  altering  the  compaction 
variables. 
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5-2  Recommendations  for  Future  Research 

This  study  demonstrates:  1)  the  dependence  of  permeability  on 
pore  size  distribution,  and,  2)  the  utility  of  pore  size  distribution 
measurements  to  infer  soil  fabric.   In  light  of  the  work  reported,  the 
following  topics  should  be  considered  for  further  research: 

1)  Improvement  of  theoretical  models  for  permeability  to 
better  account  for  the  complex  nature  of  the  pore  system 
and  the  pore  size  distribution. 

2)  Particle  migration,  piping,  and  filter  criteria  based  on 
pore  size  distribution  parameters.   Improved  pore  size 
measurement  techniques  would  be  required  for  coarse  soils. 

3)  Non-Darcy  flow  and  the  physio-chemical  factors  which  affect 
permeability. 

4)  The  influence  of  pore  size  distribution  on  other  important 
engineering  properties  such  as  strength  and  compressibility. 

5)  The  pore  size  distribution  of  field  compacted  versus 
laboratory  compacted  samples  to  assess: 

a.  the  replication  of  the  field  soil  structure  in 
the  laboratory, 

b.  the  variability  of  the  field  compacted  soil 
structure,  and, 

c.  current  field  compaction  specifications. 

6)  The  use  of  pore  size  distribution  as  a  routine  classification 
test  for  undisturbed  soils.   Conventional  classification  tasts 
of  remolded  samples  are  of  little  value  in  identifying 
structurally  sensitive  soils. 
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7)  The  relation  between  the  small  pore  size  range  Cor  parameters 
calculated  therefrom),  soil  minerology,  and  physio-chemical 
properties  such  as  plasticity  and  swelling. 
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Replicate  Pore  Size  Distribution  Curves 
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FIGURE  40    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAM.^LE  S9MW 
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FIGURE   41      REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S9HO 
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FIGURE  42    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S9HD 
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FIGURE  43    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S9HW 
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FIGURE  44    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7L0 
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FIGURE  45   REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7LD 
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FIGURE  46     REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7LW 
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FIGURE  47     REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7M0 
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FIGURE  48  REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7MD 
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FIGURE  49   REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7.VJW 
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FIGURE  50    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7H0 
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FIGURE  51     REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE   S7HD 
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FIGURE  52     REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S7HVV 
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FIGURE  53   REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S5LO 
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FIGURE  54   REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S5LD 
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FIGURE  55     REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S5LW 
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FIGURE  56    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
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FIGURE  57    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S5MD 
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FIGURE  58    REPLICATE  PORE  SIZE  DISTRIBUTION  CURVES 
FOR  SAMPLE  S5MW 
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APPENDIX  B 

Permeability  Tests 

Laboratory  permeability  measurements  of  fine-grained  soils  may 
involve  several  sources  of  errors.   Aside  from  the  major  difficulty  of 
replicating  field  conditions  in  the  laboratory,  errors  or  fluctuations 
in  permeability  measurements  may  result  from: 

1)  excess  flow  around  the  sample 

2)  leaks  in  the  system 

3)  inadequate  hydraulic  capacity  of  the  permeation  equipment 

4)  evaporation 

5)  volume  changes  during  permeation  (consolidation  or  swelling) 

6)  changing  degree  of  saturation 

7)  particle  migration 

8)  variation  in  electrolyte  concentration 

9)  temperature  fluctuations 

10)  formation  of  bacteria 

11)  general  measurement  errors. 

As  mentioned  previously,  samples  tested  in  this  study  were  compact- 
ed directly  in  the  permeation  molds  in  order  to  eliminate  excess  flow 
along  the  boundaries.   Mitchell  et  al.  (1965)  compared  the  permeability 
of  two  samples  of  compacted  silty  clay  prepared  by  the  same  technique. 
One  sample  was  compacted  directly  in  the  permeation  mold  and  tested; 
the  other  sample  was  tested  in  a  triaxial  apparatus,  surrounded  by  a 
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rubber  membrane.   Both  samples  were  tested  under  identical  back 
pressures,  and  showed  nearly  equal  permeability  values. 

A  closed  system  falling  head  permeability  test  was  used  because  it 
was  felt  to  have  several  advantages  over  the  constant  head  test.   The 
net  volume  of  flow  through  the  sample  could  be  monitored  with  the 
permeability  measurements,  allowing  for  the  detection  of  leaks, 
consolidation  or  swelling.   It  was  necessary  to  assure  that  steady  state 
flow  conditions  were  achieved  in  the  sample  for  accurate  permeability 
determinations.   The  falling  head  test  applied  low  gradients  across  the 
sample,  thus  minimizing  seepage  forces  and  consolidation  during 
permeation.   Back  pressure  could  be  readily  applied  to  the  samples 
during  permeation  to  achieve  a  high  degree  of  saturation.   The  only 
significant  disadvantage  of  the  system  was  the  long  duration  of  tests 
on  samples  of  low  permeability. 

The  permeability  tests  were  usually  run  over  a  period  of  three  to 
five  days  to  check  if  particle  migration,  variations  in  electrolyte 
concentrations  or  other  factors  caused  permeability  fluctuations  with 
time. 

B-l   Calculations 
The  permeabilities  were  calculated  from  the  falling  head  equation 
presented  by  Taylor  (1948)  as  follows: 


"  -  Ita  *"  ^  CB"» 
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where  L  is  the  length  of  sample 

a  is  the  cross-sectional  area  of  the  standpipes 

A  is  the  cross-sectional  area  of  sample 

At  is  the  time  increment 

hQ  is  the  initial  total  head 

h.  is  the  total  head  at  time  At. 

Figure  59  shows  a  sketch  of  the  closed  system  falling  head  test.   It  is 
apparent  from  the  figure  that: 


h  =  (h   +  u  )  -  (h   +  u  ) 
o     xo    o    v  yo    o 


h  =  h   -  h  CB-2") 

o    xo    yo  v  *' 


where  u  is  the  back  pressure  on  the  system, 
and 


\   "  hxl  -  hyl  (B-3) 


The  permeability   k  is   then   calculated  as; 


i  a   L      0 

k=  ATA*n 


h       -  h 
xo  yo 


ih        -   h   . 

*■   xl  yl 


(B-A) 


The  ratio  of  absolute  viscosity  y  at  temperature  T  (°C)  to  the  absolute 
viscosity  at  20  C  for  water  was  determined  from  the  following  correlation 
which  is  accurate  between  20  and  100°C  (CRC,  Handbook  of  Chemistry  and 
Physics,  1970-71,  p.  F-36) : 

loe   yT  _  1.327  (20-T)  -  0.001053  (T-20)2 
g  Honop  ~  T  +  105 


20  C 


therefore, 
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1.327  (20-T)-  0.001053  (T-20)      .        . 
yT  T  +  105  l   } 

=  10 


»2Q°Q 


The  normalized  permeability  was  then  calculated  as  follows: 


yT 


k20°C  =  y2Qoc 


K^  (B-6) 


The  net  volume  of  inflow  AV  during  permeation  is: 


AV  =  V.   -  V  (B-7) 

in    out 


AV  =  {h   -  h  .  -  (h  .  -  h  ) }  a 

xo    xl     yl  yo 

AV  =  a{h   +  h   -  h  ,  -  h  ,}  (B-8) 

xo    yo    xl  yl 


where  a  is  the  area  of  the  standpipes. 

Because  of  the  lengthy  duration  of  the  permeability  tests, 
evaporation  during  permeation  became  a  concern.   Assuming  that  the  rate 
of  evaporation  in  the  two  standpipes  is  equal,  the  evaporation  in  each 
of  standpipes  after  time  At  may  be  expressed  as  Ah.   Entering  equation 
(B-4)  with  a  correction  for  the  evaporation  of  Ah  after  time  At  yields 
the  following: 


h   -  h 

k  =  fLL.  zn  22 vo 

K   At  A    (h   +  Ah)  -  (h   +  Ah) 


which  is   equal   to: 


k-frAte 


h       -  h 
xo yo_ 


hxl  "  V 
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FIGURE  59    CLOSED  SYSTEM  FALLING  HEAD  PERMEABILITY 
TEST 
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It  is  evident  that  evaporation  has,  no  influence  on  the  calculated 
permeability  values  if  the  rate  of  evaporation  is  equal  in  both 
standpipes.   The  same  is  not  true  for  the  net  volume  of  inflow, 
however.  Applying  the  evaporation  correction  Ah  to  equation  (B-8) 
yields  the  following: 


AV  =  a{h   +  h   -  (h  .  +  Ah)  -  (h  ,  +  Ah)} 
xo    yo     yl  yl     " 

AV  =  a{h   +  h   -  h  .  -  h  ,  -  2  Ah}         (B-9) 
xo    yo    xl    yl 


Therefore,  the  net  volume  of  inflow  does  reflect  evaporation  losses 
during  permeation. 

B-2  Consolidation  During  Permeation 

For  several  of  the  preliminary  permeability  tests  performed  during 
this  investigation  a  10  to  20  psi  water  pressure  differential  was 
placed  across  the  samples  to  accelerate  the  saturation  process. 
These  samples  consolidated  up  to  1/4  inch  because  of  the  high  seepage 
forces  on  the  soil. 

Figure  60  shows  an  effective  stress  analysis  of  a  saturated  sample 
under  a  constant  back  pressure  in  a  "no  flow"  condition,  and  the  same 
sample  with  a  constant  pressure  differential  Au  across  it  to  induce 
flow.   From  the  figure,  the  change  in  effective  stress  Aa'  with  depth 
z  between  the  constant  pressure  differential  and  the  no  flow  cases  is 
equal  to: 


Ao^  =  {zy'  +  |  (ux  -  u2)}  -  {zy'} 


A°z  =  I  (ul  "  U2) 
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Act'  =  f  Au  Cb-10) 

Z     u 


where  u1  and  u„  are  the  boundary  water  pressures  on  the  top  and  bottom 
of  the  samples,  respectively 
y'  is  the  submerged  unit  weight  of  the  soil  and  equal 

t0  (Ym  "  V 

Y  is  the  saturated  unit  weight  of  the  soil 

Y  is  the  unit  weight  of  water  and 
L  is  the  sample  length. 

From  equation  (B-10) ,  at  the  base  of  the  sample  (z  =  L) ,  the  effective 
stress  increase  is  equal  to  the  pressure  differential  Au,  while  at  the 
top  of  the  sample  the  effective  stress  remains  unchanged.   For  large 
pressure  differentials,  the  sample  will  experience  increasing  consolida- 
tion and  decreasing  permeability  with  depth  and  time  until  equilibrium 
conditions  are  reached. 

Constant  head  tests  employing  high  pressure  differentials  are 
frequently  used  to  accelerate  permeability  measurements  of  clays.   Such 
tests  are  likely  to  produce  doubtful  results  as  demonstrated  above. 
High  water  pressure  differentials  are  present  in  some  earth  structures 
such  as  the  core  of  earthen  dams.   In  these  cases  it  may  be  advantageous 
to  account  for  the  consolidation  and  for  the  permeability  decreases 
which  will  take  place  with  time  as  a  result  of  the  seepage  forces. 
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APPENDIX  C 

Critical  Region  Drying  Problems 

Critical  region  drying  was  to  be  used  to  prepare  samples  for 
mercury  intrusion  pore  size  distribution  measurements.   However,  during 
preliminary  critical  region  runs,  the  soils  used  in  this  study  under- 
went a  35  to  40%  increase  in  void  ratio  as  a  result  of  the  dehydration 
procedure.   A  mineralogic  analysis  of  the  soil  revealed  that  the 
critical  region  process  caused  a  chemical  reaction  in  the  soil. 

The  critical  region  drying  procedure  employed,  was  that  used  by 
Bhasin  (1975)  and  described  by  Reed  (1977).   As  mentioned  in  Section 
1-3.3,  the  process  involves  elevating  the  pressure  and  temperature  of 
the  soil  to  3800  psi  and  380°C,  respectively. 

Figure  61  shows  a  comparison  of  X-ray  patterns  for  two  soil  types, 
each  air  and  critical  region  dried.   The  air  dried  soils  consist 
primarily  of  quartz  and  dolomite,  with  fractions  of  kaolinite,  calcite, 
and  feldspar  also  present.   Critical  region  drying  eliminated  the  kaolin 
peaks,  significantly  reduced  the  dolomite  and  calcite  peaks,  and 
increased  the  feldspar  peaks. 

Detailed  analysis  of  the  changes  in  interplanar  or  d  spacing  of 
the  peaks  between  the  air  dried  and  critical  region  dried  specimens 
indicated  that  the  mineral  being  formed  by  critical  region  drying  was  a 
form  of  plagioclase  feldspar  (ASTM  10-359  or  9-465)  .   The  reaction 


1.  ASTM  Inorganic  Powder  Diffraction  File  numbers. 
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FIGURE  61      X-RAY  DIFFRACTION  PATTERNS  OF  CRITICAL 
REGION  AND  AIR  DRIED  SOILS 
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between  the  carbonates  and  kaolintte  which  formed  this  mineral  was,  no 
doubt,  responsible  for  the  structural  changes  undergone  by  the  samples 
during  critical  region  drying. 

In  light  of  the  above  findings  it  would  be  prudent  to  monitor 
critical  region  drying  of  soils  with  mineralogic  determinations, 
particularly  when  carbonates  are  present. 
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APPENDIX  D 


Pore  Size  Distribution  Computer  Programs 
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C 
C 

C     THIS  PROGRAM  CALCULATES  AND  TABULATES  PORE  SIZE  DISTRIBUTION 
C     PARAMETERS,  ALONG  WITH  PLOTTING  PORE  DIAMETER  UERSUS  CUMULATIUE 
C     INTRUSION/GRAM,  CUMULATIUE  INTRUSION/UOLUME  AND  INCREMENTAL 
C     INTRUSION  /  UOLUME.   PROGRAM  AFTER  REED  (1377). 
C 

PROGRAM  PORE  ( INPUT, OUTPUT, TAPE5=INPUT, TAPEG=OUTPUT) 

DIMENSION  S(50),P(50),D(50),CU(50),SP(50),CUSP(50),GRAPH(451), 
CPORO(451),HISTO(451),C(50) 

READ,  NDATA 

DO  400  J=l, NDATA 
C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  BE  ABLE  TO  CALCULATE  RESULTS  FOR 
C     MORE  THAN  ONE  SET  OF  DATA  IN  A  SINGLE  RUN.   NDATA  IS  THE  NUMBER 
C     OF  SETS  OF  DATA  TO  EE  CALCULATED.   NDATA  IS  AN  INTEGER  UALUE,  AND 
C     IS  THE  FIRST  DATA  CARD  FOR  THE  SET  OF  DATA. 
C     DETERMINATION  OF  PORE  SIZE  DISTRIBUTION 
C 

READ, DATE, SN 
C  ' 

C     DATE  AND  SAMPLE  NUMBER  (SN)  ARE  TYPED  ON  THE  FIRST  DATA  CARD  FOR 
C     EACH  CURUE. 
C 

WRITE  (G,14) 
14    FORMAT  (1H1) 

PRINT,  DATE,  SN 

READ, US, DS, DM, UP, USP, USPM, PE, PF, SR 
C 

C     US=UEIGHT  OF  SAMPLE,  DS=SPECIFIC  GRAUITY  OF  SOLIDS,  DM=DENSITY  OF 
C     MERCURY,  UP^UOLUME  OF  PENETROMETER,  USP=UEIGHT  OF  SAMPLE  AND  PEN- 
C     ETROMETER,  USPN=UEIGHT  OF  SAMPLE  ,  PENETROMETER  AND  MERCURY, 
C     PE  =EUACUATING  PRESSURE,  PF=FILLING  PRESSURE, IN  MM  OF  MERCURY, 
C     SR  =STEM  READING  AT  FILLING  PRESSURE.   THE  ABOUE  UALUES  ARE  ALL  ON 
C     THE  SAME  DATA  CARD 
C 

PRINT,  US,  DS,  DM,  UP,  USP,  USPM,  PE,  PF,  SR 

US=US/DS 

UM=USPM  -  USP 

UM=UM/DM 

USA=UP-UM* ( 1 . +PE/PF)-SR 

UUO=USA-US 

UOIDR=UUO/US 

POROS=UUOAJSA 

PU=PE«UM 

UAI=PU/PF 

PRINT,  VOLUME  SOLIDS   ISA  US,/,  VOLUME  UOIDS   IS*,UUO 

PRINT, *UOID  RATIO   IS^.UOIDR./.r^POROSITY    IS*,POROS 

READ,N,M,SI 
C 

C     N=INTEGER  UALUE  FOR  THE  NUMBER  OF  LOU  PRESSURE  READINGS  MADE,  NOT 
C     COUNTING  THE  FILLING  PRESSURE  READING.   M=INTEGCR  UALUE  FOR  THE 
C     NUMBER  OF  HIGH  PRESSURE  READINGS  MADE,  NOT  COUNTING  THE  INITIAL 
C     POROSIMETER  READING.   SI=INTRUSION  READING  AT  INITIATION  PRESSURE 
C     IN  ML.   THE  THREE  UALUES  ARE  ALL  ON  ONE  DATA  CARD. 
C 

C     LOU  PRESSURE  CALCULATIONS 
C 
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WRITE  (6,5) 

5  FORMAT  (//,37X,/L0W  PRESSURE  CALCULATIONS*) 
WRITE  (G,G) 

6  FORMAT  C//,3X.*tiri  OF  HG  PRESSURE  UOL  AIR   STEM  RDG  INTRUS    I 
ANTRU/GM  CUMUL  IN  INTR/USA  CU  IN/US  DIAMETER  /) 

CU(1)  =  0. 

CUSP(l)  =  0. 

DO  100  1=2, N+l 

READ,P(I).S(I) 
C 

C     P(I)  AND  S(I)  ARE  THE  PRESSURE  AND  INTRUSION  UALUES  FOR  EACH 
C     READING,  P(I)  IS  GIUEN  IN  MM  OF  MERCURY  FOR  THE  LOW  PRESSURE 
C     INTRUSION  AND  PSI  FOR  THE  HIGH  PRESSURE.   SCI)  IS  IN  ML.   THE  TWO 
C     UALUES  ARE  TYPED  ON  ONE  DATA  CARD,  WITH  A  CARD  FOR  EACH  READING. 
C 

PS=P(I)*. 01934 

UA=PU/P(I) 

CUA=UAI-UA 

CSR=S(I)-SR 

SN=CSR-CUA 

IF  (SN.LT.O.)  GO  TO  10 

SG  =SN/WS 

GO  TO  11 
,  10    SG=0. 

SN  =0. 
11    OKI)  =  CU(I-l)  +  SG 

SP(I)  =  5N/USA 

cuspci)  -  cur,p(i-n  +spm 

D(I)^((-4.)«434.*C0S(2.rJG5G)«.145)/PS 

PRINT  1,P(I),PS,UA,S(I).SN,SG,CU(I),SP(I),CUSP(I),D(I) 
C 

C     P(I)=INTRU5I0N  PRESSURE  IN  MM  OF  MERCURY,  PS=INTRUSION  PRESSURE  IN 
C     PSI,  UA=U0LUME  OF  AIR  IN  ML.  5(I)=STEM  READING,  SN=CORRECTED 
C     INTRUSION  FOR  EACH  INCREMENT  IN  ML,  SG=CORRETED  INTRUSION  /GRAM 
C     FOR  EACH  INCREMENT,  CU( I )=CUMULATIUE  CORRECTED  INTRUSION^GRAM, 
C     SP(I)=INCREMENTAL  POROSITY,  CUSP(I)=CUMULATIUE  POROSITY, 
C     D(I)=LIMITING  PORE  SIZE  DIAMETER 
C 
1     FORMAT(X,10F10.4) 

UAI=UA 
100   SR=S(I) 
C 

C     HIGH  PRESSURE  CALCULATIONS 
C 

WRITE  (G.7) 

7  FORMAT  (//,3GX,/HIGH  PRESSURE  CALCULATIONS/) 
WRITE  (G.8) 

8  FORMAT  C//.3X, /PRESSURE  STEM  RDG  HG  COR    CHG  HGC   INTRUS    I 
BNTRU/GM  CUMUL  IN  INTRAJSA  CU  IN/US  DIAMETER  /) 

HCI=11.«.0011*UM/15000. 

DO  200  I=N+2,M+N+1 

READ, P( I). SCI) 

P(I)=P(I)+11. 
C 

C     HC  IS  CORRECTION  FOR  COMPRESSION  OF  MERCURY 
C 

HC=P(I)».0011*UiV15000. 

CHC=HC-HCI 
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CSR=S(I)-SI 
SN=CSR-CHC 

srci)  =  sn/usa 

cusp  (i)  =  cuspci-i)  +  spm 

SGH=SN/WS 

CU(I)=CU(I-1)+SGH 
C 

C     D(I)  PORE  DIAMETER  CALCULATIONS  BASED  OM  A  SURFACE  TENSION  UALUE 
C     OF  484  DYNES  /  CM  AMD  A  COMTACT  AMGLE  OF  147  DEGREES. 
C 

D(I)=((-4.)*4e4.*C0S(2.5S5S)*.145)/P(I) 

PRINT  2, P( I ) . S( I ) . HC, CHC. SN» SGH, CU( I ) , SP( I ) , CUSP( I) , DC  I ) 
C 

C     P(I)=INTRUSION  PRESSURE  IN  PSI,  S(I)=STEM  READING.  HC=MERCURY 
C     CORRECTION,  CHC=CHANGE  IN  MERCURY  CORRECTION  IN  ML,  SN=CORRECTED 
C     INTRUSION  IN  ML,  SGH=CORRECTED  INTRUSION  /GRAM,  CU(I )=CUMULATIUE 
C     CORRECTED  INTRUSION/GRAM,   SP( I ^INCREMENTAL  POROSITY,  CUSP(I)= 
C     CUMULATIUE  POROSITY) 
C 

2  FORMAT  (X, 10F10.4) 
HCI=HC 

200  si=sm 

DO  300  I=2,N+M+1 
300   C(I)=ALOG10(D(D) 

C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  CONUERT  THE  PORE  SIZE  DIAMETER 

C     TO  THE  LOG  OF  THE  DIAMETER 

C 

CALL  PLOT1  (0,4,10.5.20) 

CALL  PLOT  2  (GRAPH. 451,3. ,-2. , .4, 0. ) 

CALL  PL0T3  (1H*,C,CU.N+M) 

3  FORMAT  ClHi.35Xi54HPL.0T  OF  LIMITING  PORE  DIAMETER  US  CUMULATIUE  IN 
8TRUSIQN//) 

WRITE(G.3) 

CALL  PLOT4  (5.INTRU) 
WRITE  (G.20) 
URITE(G.4) 
20    FORMAT  (IOX.^0.01*. 17X, *0. 1*.  18X./1.0A 17X.*10.*, 17X,*100.*, 1GX 
C^IOOOO.*) 

4  FORMAT  (1H0.45X.32HLIMITING  PORE  DIAMETER  IN  MICRONS) 
CALL  PLOT1  (0,4,10,5,20) 

CALL  PLOT2  (PORO, 451.3. ,-2. , .4. 0. ) 
CALL  PL0T3  (1H», C, CUSP, N+M) 

12  FORMAT  (1H1,2GX,72HPL0T  OF  LIMITING  PORE  DIAMETER  US  CUMULATIUE  IN 
DTRUSION/UOLUNE  OF  SAMPLE//) 

WRITE  (G,12) 

CALL  PL0T4  (5.INTRU) 

WRITE  (6,20) 

WRITE(G,4) 

CALL  PLOT1  (0.4,10.5,20) 

CALL  PL0T2  (HISTO.451.3. ,-2. , .2,  0. ) 

CALL  PL0T3  (1H»,  CSP.N+M) 

WRITE  (G, 13) 

13  FORMAT  (1H1.35X.55HPL0T  OF  LIMITING  DIAMETER  US  INTRUSION/UOLUME  0 
EF  SAMPLE//) 

CALL  PL0T4  (5, INTRU) 
WRITE  (G.20) 
WRITE(G,4) 


170 


D(1)=((-4.)»4S4.*C0S(2.5G5S)».145)/(PF*. 01934) 

PRINT, /PORE  DIAMETER  AT  FILLING  PRESSURE  IS/,D(1) 
C 

C     BELOW  -  CALCULATION  OF  SPECIFIC  SURFACE  AREA  AND  HYDRAULIC  RADIUS 
C     DM  IS  MIDPOINT  DIAMETER  OF  D(I)  AND  D(I-l). 
C 

SUM1  =  0.0 

DO  500  I  =  2, N+M+l 

DM  =  SQRT(D(I)»D(I-1)) 
500  SUM1=  SUM1  +  SP(I)*4.  /DM 

1=  N+M+l 

SUM1=  SUMl+(POROS-CUSP( I ) )*4./. 01 

PRINT,  /SPECIFIC  SURFACE  AREA  IN  1/UM  IS  Z.SUM1 

SUM2=  POROS/SUM1 

PRINT, /HYDRAULIC  RADIUS  IN  UM  IS  Z.SUM2 
C 

C     BELOW  -  CALCULATION  OF  MEAN  PORE  DIAMETER  AND  SECOND  MOMENT  ABOUT 
C     THE  ORIGIN  FOR  PORES  LESS  THAN  10  MICRONS. 
C 

SUM1  =0.0 

SUM2  =0.0 

DO  GOO  1=2, N+M+l 

DM  =  SQRT(D(I)*D(I-1)) 

IF  (DM.  GT.  10.)  DM=0.0 

SUM1  =  SUM1  +  SP(I)«DM 
GOO   SUM2  =  SUM2  +  SP(I)«DM»»2 

DMEA1=  SUM1/POROS 

DMEA2=  SUM2/P0R0S 

WRITE  (6,21) 
21  FORMAT   (/,1X,/THE  STATISTICAL  MEASURES  BELOW  NEGLECT  PORES/,/, IX, 
A*WITH  A  DIAMETER  GREATER  THAN  10.  MICRONS/,/) 

PRINT,  /ADJUSTED  MEAN  PORE  SIZE  IS/, DMEA1,/, /ADJUSTED  SECOND 
BMOMENT  ABOUT  THE  ORIGIN  IS*,DMEA2 
C 

C     BELOW  -  CALCULATION  OF  MARSHALL  MODEL  PORE  DIAMETER  FOR  PORES 
C     LESS  THAN  10  MICRONS. 
C 

Cl=  0.0 

C2=  0.0 

DO  31  1=2, N+M+l 

C3=  SP(I)»*2 

AB=  I 

IF  (AB.  LT.  2.0)  GO  TO  29 

DO  28  K=1,I-1 
28  Cl=  2.«SP(I)*SP(K)+C1 
2S  Cl=  C1+C3 

DM  =  SQRT(D(I)*D(I-1)) 

IF  (DM.  GT.  10.)  DM=0.0 

C2=  C2+DM**2  »C1 
31  C1=0.0 

SC2  =  S0RT(C2) 
400  PRINT,  /MARSHALL  PORE  DIAMETER  SQUARED  IN  UM#«2  IS  Z.C2,/, 
BZMARSHALL  PORE  DIAMETER  IN  UM  IS  Z.SC2 

STOP 

END 
C     DATA  FOLLOWS  THE  7/8/9  CARD 
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C 
C 

C     THIS  PROGRAM  IS  ESSENTIALLY  THE  SAME  AS  THE  PRECEEDING  PROGRAM, 

C     EXCEPT  IT  PLOTS  A  SET  OF  CUMULATIUE  PORE  SIZE  DISTRIBUTIONS  USING 

C     THE  UERSATEC  OR  CALCOMP  PLOTTER.   FOR  EXPLANATION  OF  THESE 

C     PLOTTERS,  SEE  PURDUE  DOCUMENT  J5  CALCOMP.   PROGRAM  AFTER 

C     REED  (1S77). 

C 

PROGRAM  PORE  (INPUT, OUTPUT, TAPE5=INPUT.TAPEG=0UTPUT, PLOT) 

INTEGER  TITL 

INTEGER  TITLE 

DIMENSION  S(50),P(50),D(100),CU(100),SP(50).CUSP(100) 

CALL  PLOTS 

A  =  0.0 

READ,  NSET 
C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  PLOT  MORE  THAN  ONE  SET 
C     OF  CURUES.   NSET  IS  THE  NUMBER  OF  SETS  OF  CURUES  TO  BE  PLOTTED  AND 
C     IS  AN  INTEGER  UALUE  THAT  IS  TYPED  ON  THE  FIRST  DATA  CARD. 
C 

DO  500  K=1,NSET 

B  =  0.0 

CALL  FACTOR  (1.0) 

CALL  SYMB0L(A+1.85..5,.7/G.,  33HLIMITING  PORE  DIAMETER  IN  MICRONS, 
GO. 0.33) 

CALL  SYMBOL(A+l., 0.75,0. 7/G.. 
F50H.01      0.1       1.0      10.0     100.0    1000,0.0,50) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  X  AXIS 

CALL  SYMBOL  (A+2.0, 1.0, 0. 1,3, 0.0,-1) 

CALL  SYMBOL  (A+3.0, 1. 0, 0. 1.3, 0.0,-1) 

CALL  SYMBOL  (A+4.0, 1.0, 0. 1,3, 0.0,-1 ) 

CALL  SYMBOL  (A+5. 0,1.0.0.1,3,0.0,-1) 
C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  X  AXIS 

CALL  SYMBOL  ( A+0 .5, 1 .5, .7/G. ,  1SHCUMULATIUE  POROSITY, 30.0, IS) 

CALL  SYMBOL  (A+0. 6, 1.0, .7/G. 0.3H. 00, 0.0,3) 

CALL  SYMBOL  (A+0.G.2. 0, .7/G.0.3H. 15, 0.0,3) 

CALL  SYMBOL  (A+0. G, 3.0, .7/G. 0, 3H. 30, 0.0,3) 

CALL  SYMBOL  (A+O.G, 3.3, .7/G.0, 3H.45. 0. 0.3) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  Y  AXIS 

CALL  SYMBOL  (A+l. 0,2. 0,0. 1,3, 0.0,-1) 

CALL  SYMBOL  (A+l. 0, 3.0, 0. 1, 3, 0. 0,-1 ) 

CALL  SYMBOL  (A+l. 0, 4.0, 0. 1,3, 0. 0,-1 ) 
C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  Y  AXIS 

CALL  SYMBOL  (A+3.85, 3.8G, .085,GHLEGEND, 0.0.G) 
C     THE  ABOUE  WRITES  LEGENB  ON  THE  PLOT 

CALL  PLOT  (A+l. 0,4. 0,3) 

CALL  PLOT  (A+G.0,4.0,2) 

CALL  PLOT  (A+G. 0,1.0,2) 

CALL  PLOT  (A+l. 0,1.0,2) 

CALL  PLOT  (A+l. 0,4. 0,2) 
C     THE  ABOUE  5  CARDS  BRAN  THE  OUTLINE  OF  THE  GRAPH 
C     THESE  CARDS  FURTHER  LABEL  THE  GRAPH  AND  PROUIDE  THE  GRID. 

CALL  PLOT  (A+l. 0, 11. /3., 3) 

CALL  PLOT  (A+3.5,ll./3.,2) 

CALL  PLOT  (A+3.5,10./3.,3) 

CALL  PLOT  (A+1.0,10./3.,2) 

CALL  PLOT  (A+1.0,9./3.,3) 

CALL  PLOT  (A+1.0,3./3.,3) 
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CALL  PLOT  Cft+3.5»9./3.t2) 

CALL  PLOT  (A+3.5,8./3.,3) 

CALL  PLOT  (A+1.0.8./3..2) 

CALL  PLOT  (A+l. 0.7./3. .3) 

CALL  PLOT  (A+G.0.7./3..2) 

CALL  PLOT  (A+G.0tG./3.»3) 

CALL  PLOT  (A+1.0.G./3..2) 

CALL  PLOT  (A+1.0.5./3..3) 

CALL  PLOT  (A+6.0.5./3..2) 

CALL  PLOT  (A+G.0.4./3..3) 

CALL  PLOT  (A+1.0,4./3.,2) 

CALL  PLOT  (A+l. 0,1. 0,3) 

CALL  PLOT  (A+2. 0,1.0,3) 

CALL  PLOT  (A+2. 0,4. 0,2) 

CALL  PLOT  (A+3.0,4.0,3) 

CALL  PLOT  (A+3. 0,1.0,2) 

CALL  PLOT  (A+4. 0,1.0,3) 

CALL  PLOT  (A+4.0.7./3..2) 

CALL  PLOT  (A+5.0.7./3..3) 

CALL  PLOT  (A+5. 0,1.0,2) 

CALL  PLOT  (A+l. 0,1.0,3) 

DO  G  I0U=1,5 

CIOU=IOU 

DO  5  10=2,9 

CIOA=IO 

CIO=ALOG10(CIOA) 

CALL  SYMBOL  (A+CIOU+CIO, 1.0, .05,3,0.0.-1) 
5  CONTINUE 
G  CONTINUE 

CALL  PLOT  (A+l. 0,4. 0,3) 

READ,  NDATA 
C 

C     NDATA  IS  THE  NUMBER  OF  CURUES  THAT  WILL  BE  PLOTTED  FOR  EACH  SET  OF 
C     CURUES.   NDATA  IS  AN  INTEGER  UALUE,  AND  IS  TYPED  ON  THE  FIRST  DATA 
C     CARD  FOR  THE  SET  OF  CURUES. 
C 

C     DETERMINATION  OF  PORE  SIZE  DISTRIBUTION 
C 

DO  400  J=l, NDATA 

READ, DATE, SN 
C 

C     DATE  AND  SAMPLE  NUMBER  (SN)  ARE  TYPED  ON  THE  FIRST  DATA  CARD  FOR 
C     EACH  CURUE. 
C 

READ, US, DS, DM, UP, WSP, USPN, PE, PF,  SR 
C 

C     US=NEIGHT  OF  SAMPLE,  DS=SPECIFIC  GRAUITY  OF  SOLIDS,  DM=DENSITY  OF 
C     MERCURY,  UP^UOLUME  OF  PENETROMETER,  USP=WEIGHT  OF  SAMPLE  AND  PEN- 
C     ETROMETER,  USPM=UEIGHT  OF  SAMPLE  ,  PENETROMETER  AND  MERCURY, 
C     PE  =EUACUATING  PRESSURE.  PF=FILLING  PRESSURE, IN  MM  OF  MERCURY. 
C     SR  =STEM  READING  AT  FILLING  PRESSURE.   THE  ABOUE  UALUES  ARE  ALL  ON 
C     THE  SAME  DATA  CARD 
C 


US=WS/DS 

UM=USPM  -  WSP 

UM=WM/DM 

USA=UP-UM* ( 1 . +PE/PF )-SR 

UUO=USA-US 
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UOIDR=UUO/US 
P(jROS=UUO/'USA 

pu=pe*um 

UAI=PU/PF 

READ,N,M,SI 
C 

C     N=INTEGER  UALUE  FOR  THE  NUMBER  OF  LOW  PRESSURE  READINGS  MADE,  NOT 
C     COUNTING  THE  FILLING  PRESSURE  READING.   M=IHTEGER  UALUE  FOR  THE 
C     NUMBER  OF  HIGH  PRESSURE  READINGS  MADE.  NOT  COUNTING  THE  INITIAL 
C     POROSIMETER  READING.   SI=INTRUSION  READING  AT  INITIATION  PRESSURE 
C     IN  ML.   THE  THREE  UALUES  ARE  ALL  ON  ONE  DATA  CARD. 
C 

C     LOW  PRESSURE  CALCULATIONS 
C 

CU(l)  =  0. 

cusp a )  =  o. 

DO  100  1=2, N+l 

read. pen, sen 
c 

C     PCI)  AND  SCI)  ARE  THE  PRESSURE  AND  INTRUSION  UALUES  FOR  EACH 
C     READING,  PCI)  IS  GIUEN  IN  MM  OF  MERCURY  FOR  THE  LOW  PRESSURE 
C     INTRUSION  AND  PSI  FOR  THE  HIGH  PRESSURE.   S(I)  IS  IN  ML.   THE  TWO 
C     UALUES  ARE  TYPED  ON  ONE  DATA  CARD,  WITH  A  CARD  FOR  EACH  READING. 
C 

PS=P(I)«. 01934 

UA=PU/P(I) 

CUA=UAI-UA 

CSR=SCI)-SR 

SN=CSR-CUA 

IF  (SN.LT.O.)  GO  TO  10 

5G  =SN/WS 

GO  TO  11 
10    SG=0. 

SN  =0. 

'   ii   cum  =  cuci-n  +  sg 

SP(I)  =  SN/USA 

CUSPCI)  =  CUSP(I-l)  +SP(I) 
C 

C     DC  I)  PORE  DIAMETER  CALCULATIONS  BASED  ON  A  SURFACE  TENSION  UALUE 
C     OF  484  DYNES  /  CM  AND  A  CONTACT  ANGLE  OF  14?  DEGREES. 


C 


DCI)=(C-4.)*484.«C0SC2.565G)*.145)/PS 
UAI=UA 


100   SR=SCI) 

C 

C     HIGH  PRESSURE  CALCULATIONS 

C 

HCI=U.».0011*UM/15000. 

DO  200  I=N+2»M+N+1 

READ, PCI), SCI) 

PCI)=PCI)+11. 

HC=P(I)*.0011*UM/15000. 

CHC=HC-HCI 

CSR=SCI)-SI 

SN=CSR-CHC 

SPCI)  =  SN/USA 

CUSPCI)  =  CUSPCI-1)  +  SPCI) 

SGH=SN/WS 
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CU(I)=CU(I-1)+SGH 

DCI)=((-4.)*484.*C05(2.5G5B)*.145)/P<I) 

HCI=HC 
SOO   BI=S(I) 

READC5.201)  NL,A1,A2.W,H 
201  FORMAT  (II. IX, A4, A2, 1X.F4. 1. 1X.F3.B) 
C 

C     NL  IS  AN  INTEGER  UALUE  WHICH  DESIGNATES  WHAT  SYMBOL  WILL  EE 
C     PLOTTED  ON  THE  CURUE.   Al  AND  A2  ARE  FOR  SAMPLE  CODE.  W  IS  FOR 
C     WATER  CONTENT  AND  H  IS  FOR  PERMEABILITY  IN  CM/SEC. 
C 

DO  300  I=2,N+M+1 

D(I-l)  =ALOG10(D(I)) 
C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  CONUERT  THE  PORE  SIZE  DIAMETER 
C     TO  THE  LOG  OF  THE  DIAMETER 
300   CUSP(I-l)  =  CUSP (I) 

D(M+N+1)  =-3.0  -A 
C     DCM+N+1)  IS  THE  X  AXIS  UALUE  RELATIUE  TO  THE  ORDINATE  OF  THE  PLOT 

DCM+N+2)  =1.0 
C     D(M+N+2)  IS  THE  X  INCREMENT  FOR  EACH  INCH  OF  PLOT 

CUSPCM+N+1)  =-0.15 
C     CUSPCM+N+1)  IS  THE  Y  AXIS  UALUE  RELATIUE  TO  THE  ORDINATE  OF  THE 
C     PLOT. 

CUSPCM+N+2)  =0.15 
C     CUSPCM+N+2)  IS  THE  Y  INCREMENT  FOR  EACH  INCH  OF  PLOT 

CALL  SYMBOL  (A+3.G3.3.G4-B. .07.NL. 0.0.-1) 

CALL  SYMBOL  (A+3.85.3.G4-B, .035, 13HC0DE         W/C  =.0.0,13) 

CALL  SYMBOL  (A+3.05.3.50-B, .085.23HPERMEABILITY  =        CM/SEC. 
BO. 0,23) 

CALL  NUMBER  (A+4.25.3.B4-B, .085, Al, 0.0.2HA4) 

CALL  NUMBER  (A+4.5G.3.G4-B, .035, A2, 0.0.2HA2) 

CALL  NUMBER  (A+5.30,3.G4-B, . 085, W, 0. 0, 1) 

CALL  NUMBER  (A+4. 31,3. 50-B. .085, H, 0.0, 4HE7. 1) 
C 

C     THE  ABOUE  LISTS  THE  WATER  CONTENT  AND  PERMEABILITY  FOR  THE  PORE 
C     SIZE  CURUE 
C 

CALL  LINE  (D. CUSP, N+N, 1.1, NL) 
C 

C     THE  ABOUE  PLOTS  THE  ACTUAL  PORE  SIZE  CURUE 
C     STEPS  BELOW  PLOT  TOTAL  POROSITY  ON  THE  ORDINATE 
C 

XN=3.0*POROS/0.45 

CALL  SYMBOL  (A+l . 0, XN+1 . 0, 0. 105, NL, 0. 0, -1 ) 
C     B  IS  FOR  THE  PLACEMENT  OF  THE  PERMEABILITY  AND  WATER  CONTENT 
C     UALUES. 
400   B=B+0.335 
500   A  =  A  +  7.0 
C     A  IS  USED  FOR  THE  LOCATION  OF  THE  PORE  SIZE  PLOT 

CALL  PLOT  (0.0.0.0.339) 

STOP 

END 
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C 

C     THIS  PROGRAM  IS  ESSENTIALLY  THE  SAME  AS  THE  PRECEEDING  PROGRAM, 

C     EXCEPT  IT  PLOTS  A  SET  OF  DIFFERENTIAL  PORE  SIZE  DISTRIBUTIONS  ON 

C     .THE  UERSATEC  OR  CALCOMP  PLOTTER.   FOR  EXPLANATION  OF  THESE 

C     PLOTTERS.  SEE  PURDUE  DOCUMENT  J5  CALCOMP. 

C 

PROGRAM  PORE  ( INPUT. OUTPUT. TAPE5=INPUT. TAPE6=0UTPUT, PLOT) 

INTEGER  TITL 

INTEGER  TITLE 

DIMENSION  S(50),P(50).D(100).CU(100).SP(50),CUSP(100).DMN(50) 

CALL  PLOTS 

A  =  0.0 

READ.  NSET 
C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  PLOT  MORE  THAN  ONE  SET 
C     OF  CURUES.   NSET  IS  THE  NUMBER  OF  SETS  OF  CURUES  TO  BE  PLOTTED  AND 
C     IS  AN  INTEGER  UALUE  THAT  IS  TYPED  ON  THE  FIRST  DATA  CARD. 
C 

DO  500  K=1.NSET 

B  =  0.0 

CALL  FACTOR  (1.0) 

CALL  SYMBOL  ( A+l .85. .5. .7/6. . 33HLIMITING  PORE  DIAMETER  IN  MICRONS. 
A0.0.33) 

CALL  SYMEOL(A+l., 0.75, 0.7/6., 
F50H.01      0.1       1.0      10.0     100.0    1000.0.0,50) 
C     THE  ABDUE  PLOTS  AND  LABELS  THE  X  AXIS 

DO  15  1=2,5 

CI  =  I 

CALL  SYMBOL  (A+CI, 1.0. 0.1, 13, 0.0,-1) 

15  CONTINUE 

C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  X  AXIS 

CALL  SYMBOL  (A+0.5, 1 .5, .7/G. , 18HP0R0SITY  FREQUENCY. 30. , 19) 

CALL  SYMBOL  (A+O.G, 1 . 0, .7/G. ,3H. 00, 0.0,3) 

CALL  SYMBOL  ( A+O.G, 2. 0. .7/6. , 3H. 05, 0. 0, 3) 

CALL  SYMBOL  (A+O.G. 3.0, .7/6. . 3H. 10. 0. 0,3) 

CALL  SYMBOL  (A+. 6,3.3, .7/6. ,3H. 15, 0.0,3) 
C     THE  ABOUE  PLOTS  AND  LABELS  THE  Y  AXIS 

DO  16  1=1,4 

CI  =  I 

CALL  SYMBOL  (A+1.0, CI, 0. 1, 15, 0.0.-1) 

16  CONTINUE 

C     THE  ABOUE  PLACES  TICK  MARKS  ALONG  THE  Y  AXIS 

CALL  SYMBOL  (A+4.40,3.80, .085, 6HLEGEND. 0.0.6) 
C     THE  ABOUE  WRITES  LEGEND  ON  THE  PLOT 

CALL  PLOT  (A+1.0, 4. 0,3) 

CALL  PLOT  (A+6.0.4.0,2) 

CALL  PLOT  (A+G.0,2.5,2) 

CALL  PLOT  (A+4.0,3.0.3) 

CALL  PLOT  (A+1.0, 3. 0.2) 

CALL  PLOT  (A+1.0, 2. 0,3) 

CALL  PLOT  (A+G.0,2.0,2) 

CALL  PLOT  (A+G. 0,1.0,3) 

CALL  PLOT  (A+1.0, 1.0.2) 

CALL  PLOT  (A+1.0, 4. 0,2) 

CALL  PLOT  (A+2.0,4.0.3) 

CALL  PLOT  (A+2. 0.1.0.2) 

CALL  PLOT  (A+3.0.1.0,3) 

CALL  PLOT  (A+3.0,4.0,2) 


176 


CALL  PLOT  (A+4.0.4.0,3) 

CALL  PLOT  (A+4. 0.1.0,2) 

CALL  PLOT  (A+5. 0,1.0,2) 

CALL  PLOT  (A+5. 0,2.5, 2) 

CALL  PLOT  (A+G.0,2.5,3) 

CALL  PLOT  (A+6.0,1.0,2) 

CALL  PLOT  (A+l. 0,1.0,2) 

DO  18  I0U=1,5 

CIOU=IOU 

DO  17  10=2,9 

CIOA=IO 

CIO=ALOG10(CIOA) 

CALL  SYMBOL  (A+CIOU+CIO, 1.0. .05, 13, 0.0.-1) 

17  CONTINUE 

18  CONTINUE 

CALL  PLOT  (A+l. 0,1.0.3) 

DO  20  10=1.3 

CIO=IO 

DO  19  IOU=1.4 

CIOU=IOU 

CI0U=CI0U/5. 

CALL  SYMBOL  (A+l .0. CIO+CIOU, . 05, 15, 0. 0,-1 ) 

19  CONTINUE 

20  CONTINUE 

CALL  PLOT  (A+l. 0,4. 0.3) 

READ.  NDATA 
C 

C     NDATA  IS  THE  NUMBER  OF  CURUES  THAT  WILL  BE  PLOTTED  FOR  EACH  SET  OF 
C     CURUES.   NDATA  IS  AN  INTEGER  UALUE.  AND  IS  TYPED  ON  THE  FIRST  DATA 
C     CARD  FOR  THE  SET  OF  CURUES. 
C 
C 

C     DETERMINATION  OF  PORE  SIZE  DISTRIBUTION 
C 

DO  400  J=l. NDATA 

READ. DATE. SN 
C 

C     DATE  AND  SAMPLE  NUMBER  (SN)  ARE  TYPED  ON  THE  FIRST  DATA  CARD  FOR 
C     EACH  CURUE. 
C 

READ.  US. DS. DM, UP. WSP. WSPM, PE. PF. SR 
C 

C     US=WEIGHT  OF  SAMPLE.  DS=SPECIFIC  GRAUITY  OF  SOLIDS.  DM=DENSITY  OF 
C     MERCURY,  UP=UOLUME  OF  PENETROMETER,  USP=UEIGHT  OF  SAMPLE  AND  PEN- 
C     ETROMETER.  W5PM=WEIGHT  OF  SAMPLE  ,  PENETROMETER  AND  MERCURY. 
C     PE  =EUACUATING  PRESSURE.  PF=FILLING  PRESSURE. IN  MM  OF  MERCURY, 
C     SR  =STEM  READING  AT  FILLING  PRESSURE.  THE  ABOUE  UALUES  ARE  ALL  ON 
C     THE  SAME  DATA  CARD 
C 


US=WS/DS 

UN=USPM  -  WSP 

UM=UM/DM 

USA=UP-UM* ( 1 . +PE/PF ) -SR 

UUO=USA-US 

UOIDR=UUO/US 

P0R0S=UU0/U3A 

PU=PE*UM 

UAI=PU/PF 
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READ,N,M,SI 
C 

C  N= INTEGER  UALUE  FOR  THE  NUMBER  OF  LOW  PRESSURE  READINGS  MADE.  NOT 

C  COUNTING  THE  FILLING  PRESSURE  READING.   M=INTEGER  UALUE  FOR  THE 

C  NUMBER  OF  HIGH  PRESSURE  READINGS  MADE,  NOT  COUNTING  THE  INITIAL 

C  POROSIMETER  READING.   SI=INTRUSION  READING  AT  INITIATION  PRESSURE 

C  IN  ML.   THE  THREE  UALUES  ARE  ALL  ON  ONE  DATA  CARD. 
C 

C  LOW  PRESSURE  CALCULATIONS 
C 

cum  =  o, 

CUSP(l)  =  o. 
C     D(I)  PORE  DIAMETER  CALCULATIONS  BASED  ON  A  SURFACE  TENSION  UALUE 
C     OF  484  DYNES  /  CM  AND  A  CONTACT  ANGLE  OF  147  DEGREES. 

D(1)=-4.*484.*COSC2.5G56)*0. 145/(0. 01334*PF) 

DO  100  1=2, N+l 

READ, P(I), SCI) 
C 

C     PCI)  AND  SCI)  ARE  THE  PRESSURE  AND  INTRUSION  UALUES  FOR  EACH 
C     READING,  P(I)  IS  GIUEN  IN  MM  OF  MERCURY  FOR  THE  LOW  PRESSURE 
C     INTRUSION  AND  PSI  FOR  THE  HIGH  PRESSURE.   S(I)  IS  IN  ML.   THE  TWO 
C     UALUES  ARE  TYPED  ON  ONE  DATA  CARD,  WITH  A  CARD  FOR  EACH  READING. 
C 

PS=PCI)». 01934 

UA=PU/P(I) 

CUA=UAI-UA 

CSR=S(I)-SR 

SN=CSR-CUA 

IF  (SN.LT.O.)  GO  TO  10 

SG  =SN/HS 

GO  TO  11 
10    SG=0. 

SN  =0. 

ii   cum  =  cuci-i)  +  sg 

SP(I)  a  SNAJSA 

cuspm  =  cusp(i-i)  +sp(d 

Dm  =  ((-4.)*484.*C0SC2.565G)*.145)/PS 
UAI=UA 

ioo  sR-sm 

C     HIGH  PRESSURE  CALCULATIONS 

HCI=11.*.0011*UM/15000. 

DO  200  I=N+2,M+N+1 

READ, PCI). SCI) 

PCI)=PCI)+11. 

HC=PCI)«.0011*UfV15000. 

CHC=HC-HCI 

CSR=S(I)-SI 

SN=CSR-CHC 

SPCI)  =  5N/USA 

CUSP(I)  =  CUSPCI-1)  +  SPCI) 

SGH=SN/WS 

CUCI)=CUCI-1)+SGH 

D(I)=((-4.)«484.*C0S(2.5G5G)*.145)/P(I) 

HCI=HC 
200   SI=SCI) 

READC5.201)  NL.A1,A2,W,H 
201  FORMAT  CI1. IX, A4, A2, 1X.F4.1. 1X.F9.8) 
C 
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C     ML  IS  AN  INTEGER  UALUE  WHICH  DESIGNATES  WHAT  SYMBOL  WILL  BE 

C     PLOTTED  ON  THE  CURUE.   Al  AMD  A2  ARE  FOR  SAMPLE  CODE,  W  IS  FOR 

C     WATER  CONTENT  AND  H  IS  FOR  PERMEABILITY  IN  CM/SEC. 

C 

C     THE  PURPOSE  OF  THIS  DO  LOOP  IS  TO  CONUERT  THE  PORE  SIZE  DIAMETER 

C     TO  THE  LOG  OF  THE  DIAMETER 

DO  300  I=2.N+M+1 

DMN(I-1)=(ALOG10(D(I))+ALOG10(D(I-1)))/2.0 

SP(I-1)=  SPCI) 
300  CUSP(I-1)=  CUSP(I) 

DMN(M+N+1)=  -3.0-A 

DMN(M+N+2)=  1.0 

SP(M+M+1)=  -0.05 

SP(M+M+2)=  0.05 

CALL  SYMBOL  (A+4.20,3.6-B, .07, NL, 0.0, -1) 

CALL  SYMBOL  (A+4.40.3.55-B, .085, 11HSAMPLE  CODE, 0.0,11) 

CALL  NUMBER  (A+5.30,3.55-B, .085, Al, 0.0, 2HA4) 

CALL  NUMBER  (A+5.G0.3.55-B, .085, A2, 0.0.2HA2) 
C 

C     THE  ABOUE  LISTS  THE  WATER  CONTENT  AMD  PERMEABILITY  FOR  THE  PORE 
C     SIZE  CURUE 
C 

CALL  LINE  (DMN,SP.N+M,1,1,NL) 
C 

C     THE  ABOUE  PLOTS  THE  ACTUAL  PORE  SIZE  CURUE 
C 

C     B  IS  FOR  THE  PLACEMENT  OF  THE  PERMEABILITY  AND  WATER  CONTENT 
C     UALUES. 
400  B=B+0.2 
500  A=A+7.0 
C     A  IS  USED  FOR  THE  LOCATION  OF  THE  PORE  SIZE  PLOT 

CALL  PLOT  (0.0,0.0,939) 

STOP 

END 
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APPENDIX  E 


Purdue  University  Negative  Numbers  for  Photographs 

Figure  Negative  Number 

1.  76449 

4.  75529-24 

5.  75529-34 

6.  76455-2 
7#  76455-10 

8.  76455-26 

9.  75529-25 
10.  75529-28 
61.  76481 
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